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Abstract 

It is possible, though technically difficult, to produce beams of free 
electrons that exhibit beats of a quantum-mechanical nature. One may readily 
think of the following two applications: 

1) The generation of electromagnetic radiation, e.g. light, based on the 
fact that the beats give rise to alternating charge and current 
densities. 

2) A frequency shifter, based on the fact that a beam with beats constitutes 
a moving grating. When such a grating is exposed to external radiation 
of suitable frequency and direction, the reflected radiation will be 
shifted in frequency, since the grating is moving. A twofold increase 
of the frequency is readily attainable. 

In this report we show that 

1) . It is impossible to generate radiation, because the alternating electro- 

magnetic fields that accompany the beats cannot reform themselves into 
freely propagating waves. 

2) The frequency shifter is useless as a practical device, because its 
reflectance is extremely low for realizable beams. 

Of course, when the work on this grant was started, one could not foresee 
that, at a later daxe, somebody (like this writer) would arrive at such 
disappointing conclusions. It simply takes a lot of time to go through the 


analysis. 
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The Main Fart ( with page numbers Mi, MS, etc.) of this report and its 
Appendix (with page numbers Al, AS, etc.) contain a complete record of how 
this writer arrived at his conclusions. The report is addressed to someone 
who has some application for electron beats in mind, and who wishes to find 
out whether or not his idea has any chance of success. 
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Are there any useful applications for electrons that 
are made to exhibit quantum-mechanical beats? 


by 

Roland M, Lichtenstein 
Professor of Physics, Emeritus 
Rensselaer Polytechnic Institute 
Troy, NY 12l8l 
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Preface and Guide to the Reader . 

This report concerns the use that one might make of free electrons 
with quantum mechanical beats. The Main Part explains what the beats and 
what the contemplated applications are. Two applications are investigated 
in detail, the generation of light and a selective moving mirror. The con- 
clusions are; 

1. It is impossible to generate light. 

2. The mirror device is useless, because its reflectance is 
essentially zero. 

The negative results in these two instances induce the writer to believe that 
there are no other possible applications and that further work should be 
suspended, until somebody else comes up with a fruit fill idea. 

The reader should read the Main Part first. Then he may consult the 
Appendix, where the physical theories and mathematical methods are explained. 

This writer'is well aware of the fact that this report is not easy to 
read. He apoligizes to the reader for making his task so difficult. But it 
seems to be inherent in the subject that so many strands of argument have 'go 
be employed. Perhaps after two or more revisions, a more readable report 
might be produced. But in view of the negative results, such revisions are 


not worth the effort 
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Main Part of the Report 




X. Introductj.cn * 

A few years ago, H. Schwarz 'in Phys, Rev. Letters 42, ll4l (1979) 
and 43, 238 (e) (1979) suggested that quantum-mechanical beats of electrons 
could be exploited in a novel form of a free-electron laser. In this intro- 
duction we explain the meaning of the word ’’beats”. Throughout this report, 
we shall describe electrons with the formalism of non-relativistic quantum 
mechanics, in the style of Schrffdinger. A relativistic theory, in the style 
of Dirac, is not needed, because the speeds of the electrons are much smaller 
than the speed of light, c, in the device proposed by Schwarz. 

A single electron is described by a complex Schrodinger wave function 
i|<(x,ct), where r is the position vector, and t is the time. We find it more 
convenient to express the time by means of the product ct, instead of by t 


itself, because ct s r ;d the position r corns in Ibe same unit, e.g, the meter. 


This usage facilitates dimension checks. For the same reason, we use the 
product cp, whenever we deal with the electric charge density p, because cp 
and the electric current density J come in the same unit, e.g. the amp m" 2 . 
The wave function determines the associated charge and current densities, 
which are given by 

cp(r ,ct) = (-ec)\j» , (l) 

J(r,ct) = (-ec) — -/g(-ii|> ^ + iiJiVij) ) , (2) 

me 

Here, i - >£T, (-e) is the electronic charge (e = + 1.6 x 10" 19 coul), m is 

/ df 

the mass of an electron, is Planck’s constant divided by 2 tt, — — = 

me 

3.86 x 10“ 11 cm is the Compton wavelength divided by 2 tt (which implies that 
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is the Compton wave number); the star denotes the complex conjugate, and 
the symbol V denotes the gradient. (Later on, the symbols V» and Vx denote 
the divergence and the curl.) The wave function is normalized by the 
normalization condition 

/At r iM * 1 . (3) 


where the integration is over the entire position space, and where Ax^ is 
the volume element in this r-space. 

An electron is called "free" when no electromagnetic forces act on it. 
In that case, the wave function \p satisfies the simolest form of the 
SchrHdinger equation 


i Ji . . % 4 . 

3ct r 


me 


A particular solution of (1+) is the plane wave 


\l> .. = const, exp (i ior - — J$tc 2 ct) , 

T plane * me * 


(4) 


(5) 


where ic is the vectorial wave number, and tc 2 = k • *c . The general solution of 
(4) is the wave packet 


, _ 

^(r,ct ) = /Ax K F(K)exp(iK-r - — ^ic 2 ct) , (6) 

where F(k) is a complex amplitude factor that depends on the wave number ic. 

The integration is over the wave number space. Ax = Ate Aic Ax is the 

k x y z 

volume element in this ic-space. The structure and development of a wave 

packet is completely specified by the function F(ic) in the ic-space. 

Let us introduce the spatial Fourier transform t|>(‘K,ct) (pronounced 

psi hat) and the temporal Fourier transform ijj(r, —) (pronounced psi tilde) 

of the wave function ^(r,ct). (Whenever we deal with a frequency w, we use 

the ratio — because — and ic come in the same unit, e.g. m” 1 )* These transforms 
c c 


are defined by 


t 




M«,ct) * jh t e“ lK * r i|)(r,ct) , 


*0 0 * p 4 0t , 
* <ir'y S 


(7) 


(The integration is over the position space, At^ ■ AxAyAz is the volume 
elemeny in this r-space. ) and 


-i — ct 

#(r, —) = / +a> dct e c 4»(r,ct) . 
ct*-®' 

The theory of Fourier transforms yields the inversion formulas 

1 r i<*r 


and 


'M?,ct) * J27P“ / At k* ?(*,ct) , 
<|>(r,ct) 


_1 

2tt 


„ i - ct 

r d 7 a C V(r. ~> • 


( 8 ) 


(9) 


( 10 ) 


A consequence of the Fourier transform theorems are the Parseval relations 


/ A V* * s T2^T sr * * 


and 


/ + “ dctip% = —• /** d ^ . 

ct=-» 


( 11 ) 


(12) 


— - OB 

c 


When we compare (6) and (9) we see that 


$(<,ct) = F(k)exp(- i tc 2 ct) . 


(13) 




Then the Parseval relation (ll) combined with (3) gives us the following 
normalization condition for the amplitude factor F(ic). 

/At F*(k.)F(k) = (2ir) 3 (lU) 

(Whereas a. wave packet can be normalized, the plane wave (5) cannot be 
normalized, because the latter extends with constant absolute value 
throughout the infinite r-space„) 

Now let us write the charge density co(r,ct) in terms of the temporal 
Fourier transforms ij) and ij) . For the \p in (l) we use (10). For the ip 
in (l) we use a similar formula, but with a renamed integration variable. 


t>| o 
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Thus 


/<?,= t) - £ _ r a f e' 1 c Ct Vff,f ) • 


(15) 


The reason for this renaming is fehat we wish to write the product of two 
integrals aa the double integral over the product of the two integrands. 
Here confusion would result, had we not renamed. Then (l), (10) and (15) 
yield 


c p (?,ct) « {ffyr // 4 f 4 f 




♦<r#V(r.£> . 


cx 8 

We introduce two new integration variables — and — defined by 

c c 

a _ to 
c "* c “ c 

6 = 
c 

and vice versa. 


(16) 


(IT) 


w _ a + 
c ” c c 


£ 

C 


( 18 ) 


to 


The Jacobian of this transformation is unity. Furthermore, as — and — 

c c 

range over the entire (— , — )-plane, 2L and 0 range over the entire (— . — )- 


plane. Thus, 


cp I 


(?, ct ) * // d f d ® e c $(r, 7 + 7) , 


. ct . 
i — ct 


c ' .c c ’c 

or, after a change of notation for the integration 01 * dummy variables, (back 


to — and — ) 

° C ( \ i — ct 

cp(r,ot) = / d — e 


i (19) 


b) 


2tt 




— - 09 

C 


If we compare (19) with (10), the latter equation being applied to cp instead 
of we see that the temporal Fourier transform cp(r, of the charge 


density cp(r,ct) i 3 given by 
w» (-ec ) i 

op(r ’ -> * -sr _ i 


° r P °°» Quaut'? 


I ‘ 4 f »(r. » ♦ *)**<?. *> , 


where, of course, 

i ~ ct 

/ » , » * 1 r+*° uj c — t)i » 

cp(r,ct) * — - d-e cp(r, -) . 


— GO 


It is helpful to express the function $(r t — ) , which occurs in (20 ), 

c 

in terms of the amplitude factor F(k). We start with ( 6 ) .and use spherical 
polar coordinates k * |ic|, 0, $ in the < -space. Then 


At k * d*citd 0 K,sin 0 d$ * dKK 2 dO , 


where 


sK = d 0 sine d* 


is the solid angle element. Then ( 6 ) becomes 

tKr,ct) = r g \ 3 f* die k exp(-i ~^K 2 ct) dficF^e****, < 

K =0 < 

where ^ denotes the surface integral over the sphere of radius (c . We 

K 

express icin terms of the frequency “ , i.e. by way of the equation 

c 

a. _ . i 

c me 


Then (22) becomes 

i — ct , 

. /- . . 1 to . u» c me 1 

* ’ “ 2* a / » T?7^ 


KP(K)el 


^ me co 

" V- 2 T? 7 


We compare (23) with (10) and conclude that 


*<?. P TP'# anKF ^> el,: ' r f « f i 0 


ifJ'r , oi 


for t — ' 

1C = 

J -df C 


= 0 for - > 0 
c — 
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As an example for an application of equation (25), we consider an 
approximately mono-ar.ergetic wave packet, for which T(k.) * 0 unless k. lies 
on a very thin spherical shell (with center at Ie * <5) in the k -space. Then 
(25) tells us that #(?* ^) * 0 unless 7 lies in a very narrow interval I 
that straddles the frequency which is related to the radius of the shell hy 
(23). As a consequence, the integral in (20) is zero, unless — lies in a 

w 

very narrow interval that straddle's the frequency zero. The reason is that 
one of the factors in the integrand is definitely zero when the other is non- 
zero, unless 7 * 0, so that the whole integral is zero. (— + — and — cannot 
c ° c c c 

both lie in the narrow interval I that we mentioned above, unless — * 0.) 

c 

Thus, for an approximately mono-energetic wave packet, the temporal Fourier 

transform cp(r, 7} of the charge density p(r,ct) vanishes unless — lies in a 
c c 

very narrow interval that straddles the frequency zero. We express this state 
of affairs briefly, by saying that cp(r,ct) is apprximately d.c. ("d.c." for 
’’direct current”). 

On the other hand, for a wave packet that is not mono-energetic, equation 

(25) permits $(r, ^j-) to be non-zero over an extended range of — . And then 

(20) permits cp(r, 7) to be non-zero over an extended domain of the frequency 

7 , beyond a very, narrow interval around — =» 0. We express this state of 
c c 

briefly by saying that ep(r,ct) has alternating components or beats. Specifi- 
cally, the beat at a given frequency — is the product of the temporal Fourier 

0 

transform cp(r, ^) , regarded as a function of the position r, and the time 
i w ct C 

factor e c . The first factor is the (complex) amplitude function (it 
depends on r) of the beat. The second factor produces the alternating or 
a, c. behavior ("a.c.” for ’’alternating current"). According to equation (21), 



* 
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the charge density cp(r,ct) is a superposition of beats. 

Everything that ve have said about the charge density cp(r,ct) and its 
beats ^(r, ^e 1 c ct can be extended - mutatis mutandis - to the ciuftftt 
density J(r,ct) and its beats J(r, ^)e c ct . The analogs of equations (20) 
and (21) are 

ft*. 7) • r ^ f |(-i ♦*<?. ! + t> * 


— 3C_<J5 

c 


♦*?(?, f + f)W # (5,f» . 


where, of course, 


(i) 

— a_m 

C 


i - Ct 

c - C * 


(26) 


(2T) 


U 

To summarize; Both the charge density cp(r,ct) and the current density 

i w ct 

ir(r,ct) are linear superpositions of the beats c?i(r, — )e c v and 
i{?, £)e* t ct , as shown by equations (21) and (27). The complex beat ampli- 
tudeB cp(?, and j(r£) are given by (20) and (26). The amplitude J(r. £), 
which occurs in the latter two equations can be calculated according to (25). 
For an approximately mono-energetic wave packet, there are only beats of very 
low frequency. In the limit of a strictly mono-energetic wave packet (it 
cannot be normalized, but is a useful construct.), the charge density and 
current density do not depend on the time ct; they are strictly d.c., For a 
wave packet that is not mono-energetic, the charge density and the current 
density may contain beats of appreciable frequency. The word "beat" refers 
to the fact that, as shown by (20) and (26) the beat amplitudes cP (r,-) and 
J(r, £L) ^ Q f second order in the temporal Fourier transforms of the wave 
function, and that the beat frequency £ is the difference (“ + - § of the 
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two frequencies tint occur in the integrands of (£0) and (26). 

In the device th-vt was proposed by H. Schwarz a swarm of electrons 
is employed, each exhibiting beats in the optical frequency range, Since 
the charge and current densities cp(r,ct) and J(r,ct) are the sources of 
the electromagnetic field, and since both densities contain beats in the 
optical frequency range, one might expect that light is produced in Schwarz's 
device. 

However, this expectation is illusory. Two authors, A. Peres (Phys. 

Hev. A, 2627 (1979) and M. Peshkin (ibid, page 2629), attempted to prove 
that free electrons, even though they were made to exhibit beats, cannot 
produce electromagnetic radiation. But these proofs are not very convincing 
to this writer, because they employ only verbal arguments. (One of these 
arguments, namely the use of the superposition principle, is wrong, because 
beats involve products of wave functions, whereas the invoked principle 
applies only to linear expressions.) We, therefore, supplement the work of 
Peres and Peshkin with a more detailed mathematical treatment, which - though 
a bit tedious - can demonstrate exactly where the attempt to generate light 
fails. This analysis, which is based on the principles derived in the 
Appendix, will be presented in Section 3. But there we shall have to make 
use of the decomposition of a vector field into its longitudinal and trans- 
verse part. This matter will be explained in Section 2. 

We conclude the present section with some remarks about the physical 
meaning of cp(r,ct) and j(r,ct) as given by equations (l) and (2). These 
quantities are really the quantum-mechanical expectation values of the charge 
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and current density. But one can show, as it is done in the Appendix, 
that these expectation values are the sources of the electromagnetic field 
if this field is treated classically, while the electrons are described 
quantum-mechanically. The use of quantum mechanics for the electrons and 
of classical physics for the electromagnetic field is called the semi- 
classical method. It is described in the Appendix. We use it throughout 
this report, because we deem it accurate enough for our purposes, and 
because it is more easily handled than a fully quantized theory, in which 
also the electromagnetic field would be described quantum-mechanically. 

II. The Decomposition of a Vector Field Into Its Longitudinal and 
Transverse Part . 

We start with a vector field, e.g. the current density J(r,ct) and 

A 

take its spatial Fourier transform J(ic,et) (pronounced "Jay hat") defined 
by 

j(ic,ct) = /Ax^e - ** ^j(r,ct) . (l) 

This definition is the vectorial analog of (1,7). By the Fourier integral 
theorem, we have the inversion formula 

J(r,ct) = /Ax^e 1 * ** j( ,ct) , (2) 

the analog of (l,9). Now we define the longitudinal part J L (r,ct) of 
J(r,ct) by way of its spatial Fourier transform J (ic,ct). The definition of 

A 

J L (ic,ct) is 

j L (S,ct) = ^ 7 - KK*j(<,ct ) , for if ? 0 . (3) 

A A 

This equation states that J^(ic,ct) is the projection of J(ic,ct) onto an axis 
parallel to the wave number vector ic. For 1c « 0, the definition (3) becomes 
indeterminate, so that we must use a separate definition for this case. We 



MU 


define 3^(?j,ct) by* 
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J T (0*ct) - j j(0,ct ) . 


1 L' v kV * u * " 3 w ' w » v ' u / • (4) 

Equation (4) reaults from (3) if ve first take the average of (3) as 5 

ranges over the surface of a small sphere in K-space (vith its center at 

the origin B)$ and then ve let the radius of this sphere go to zero. The 

1 — -* 
factor ^ arises from the fact that the average <kk> of the tensor kk is 

equal to -jk 2 !, where I * II 4 JJ + RR is the unit tensor or idemf actor. 

A 

The definition of J^(K,ct) permits us to obtain J^(r,ct) by taking the Inverse 
Fourier transform, in analogy to the general formula (2). We have 
J L (r,ct) “ fg^yy /At^e^'^Oc^ot) B 7 g yp ^ A \- eik * r WK »?(* ,ct ) . (?) 

We assume that the convergence of the integrals in the K-space is good 
enough so that we may take spatial derivatives of ?-(r,ct)j such as the 
divergence 9*3^ and the curl $ x J^» by taking the derivatives of the inte- 
grand. Then (V n denotes the gradient vith respect -to r) 




l “ ! -<V lK ‘ r) • 


1 — 


and 




1 — . - 


And since 


we obtain 


V e 

r 


dx *r 








T^* 1 ' j *•?(? ,ct) 

k 


( 6 ) 


(T) 


and 


V x J L » 0 


( 8 ) 


Equation (8) states that the longitudinal part of a vector field is irrota- 


tional. 


! 
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If we take the divergence of (2) we obtain °00ft 1$ 

* ^ “ j 7fAte k * r i* •?(*', ct). 

On comparing (T ) and (9)* we see that 


V«J, 


V«J . 


(9) 


( 10 ) 


Inspection of equation (?) and (9) shows that the Fourier transforms 
/N /N 

(V»j_ )(tc,ct) and (V*j)(K,ct) are given by 

1 A* A 


(V*J l )(k ,ct)«(w) (£,ct) = iK«J(K,ct) . 
Then equation (5) may be written in the form 

J L (r,ct) = - /4 V i " r if (M)(8,ct) 


( 11 ) 


or 


J L (r,ct) * - Viji , 


with 


/\ 


(12a) 


(12b) 


ij»(r,ct) = /Ax^e rtK r ^-r (^•j)(ic,ct) . 

(There is no connection between the \p of equations (12a, b) and the SehrBdinger 
wave functions of Section 1. ) 

1 1 

Now, the inverse Fourier transform of -y is as one Jan easily show by 
direct calculation. Thus, by (12b), the Fourier transform of i)>(r,ct) is the 
product of two Fourier transforms, that of -j— — and that of Hence, by the 

convolution theorem of Fourier transforms, i|>(r,t) is the convolution of -j~- and 
V»«T, Thus, (12) may be written in the form 


J L (r,ct) = - Vi|» , 


with 


lp(r.ct) = /Ar g • (V*«T) (s,ct ) . 

! s-r * 


(13a) 


(13b) 


Equations (13) may be found in many textbooks, e.g. Morse and Feshbach, Methods 
of Theoretical Physics, McGraw-Hill 1953, pages 53, 5^ or Panofsky, and 


L 
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Phillips, Classical Electricity and Magnetism, Addison-Wesley 19&2, page 2. 
However, we find it more convenient to define 3 L (r,ct) by way of its Fourier 
transform (3). 

We are finished with the discussion of the longitudinal part J^(r,ct). 
The transverse part J T (r,ct) is defined by the requirement that both parts 
add up to J(r,ct), i.e. that 

j(r,ct ) e Jj. (r ,ct ) + ^(r,et ) . (l4) 

Then, a similar equation, namely 

A A A 

?(Jc,ct) = j (£,ct) + J (ic.ct) (15) 

must hold for the Fourier transforms, so that 

S T <£,ct) * j(fc,Ct) - ;^KK-j(£,ct) = 

1 _* - 1 

= - £-jr < x. (< x j(£,ct ) , for <* t 0 , (l6) 

and 

J T (o,ct) =~j(o,ct) . (IT) 

A A 

Equation (l6) shows that J t (k ,ct) is the projection of J(k ,ct) onto a plane 
that is normal to 

The remainder of the discussion is pretty much a repetition of what was 
done for J^. So we can be brief. The analogs of equations (5), (T), (8) are 

j T (r,ct) = j-p- /AT K e iK -* r J T («,ct) = 

= ipr (-K x (k. x j(£,ct))) , (l8) 

$ x /A Ik x, j(ic,ct) , (19) 

v • J T = 0 . (20) 

Equation (20) states that .f T (r,ct) is solenoidal. The analogs of equation 
(10) and (ll) are 
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V x 


7 x j , 


/\ . ' .X .. 

(V x J T )(<,ct) * (9 x if)(k,ct) * lie x j(E,ct) . 
Instead of (12), we have 


( 21 ) 

(22) 


with 


J T (r,ct) * 7 x C , 


(23a) 


C(r,ct ) * J^TyT / AT K eiK *^ £7 (^ * J)(«,ct) . (23b) 

When we take the divergence of (23b) by differentiating the integrand and 
when we also use (22), we see that 


7*C = 0 , (24) 

which implies that C(r,ct) is purely transverse. (Because of (l3), 7*C = 0 
yields = 0.) 

Finally, the analog of ( 13 ) is 

J T (r,ct) = V x C , (25a) 


with 

C(r,ct ) = -r- /At — (7 x J)(5,ct) (25b) 

|s-r| 

This is a well-known relation. [See the references quoted after ( 13 ) . 1 

Let us consider two x*eal vector fields 5?(jr,ci) and G(r,ct), We take the 
longitudinal part of the first and the transverse part of the second field. 
Then we consider the scalar product • f? T . Now, by (13), F^ is the 
negative gradient of some scalar field t|/, and, by (25), G^ is the curl of 
some vector field C. Then 

= -(7t|»)-(7 x C) = - 7 -(^V x c) + i|i7»(7 x c) , 
or, since the divergence of a curl is zero, 

F l -G t = - 7- (^7 X C) 
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We Integrate this equation over the interior of some large sphere of radius R, 

and convert the integral of the right-hand side into a surface integral, by 

virtue of the theorem of Gauss. Thus 

/ A t x C 

sphere r u 

vhere Af is the vectorial surface element. We assume that V • F and V x G 
vanish outside a finite region of space, as it usually happens in practice. 
(Somewhat less stringent requirements will suffice for the argument that 
follows. For instance, V»F and ^ * 5 should tend to zero for |r| -*• * at a 
sufficiently fast rate.) Then, by ( 13 ) and (25), tjf and G are of order R" 1 
for R ■+■ “. Then 7 x. C is of order R” 2 , and the surface integral in (26) is 
of order R" 1 . Thus, when we go to the limit R ®, equation (26) becomes 

/ At A' b t * 0 > (27) 

which is an important orthogonality relation. 

Another demonstration of (27) is by way of the Parseval relation 
/Ar r F L (r,ct)*G T (r,ct) = /Ar k F L (-K,ct)*G T (K.,ct) . (28) 

But the definitions (3) and (l6) show that the integrand of the right-hand 
side of (28) is zero (except at the single point it. = 0; this exception does 
not affect the integral.). 

A consequence of (27) is: When a vector field F is zero everywhere, then 
its longitudinal and transverse part must be zero individually. For, with 
F - F ^ + F t , we have (if F = 0) 

/At F*F t + /A t F-'F- + 2/At F *F m = 0 
The third integral is zero, by (27). The remaining two integrals (their 
integrands are never negative) can add up to zero only if individually 
F_ = 0 and F-, = 0. In the same way one shows that a longitudinal field F T 
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and a transverse field G^ can never add up to zero. For we have 

/ iT r<VV-<VV * / 4t AA * / 4 t A ,5 T + 2 / 4 t A' 5 T • 

By (27), the third integral on the right is zero. The other two are 
definitely positive, unless F^ * 0 and « 0 everywhere. Therefore, if 
and only if these two equations are satisfied, can F^ + G^ he zero everywhere. 

We also remark that the longitudinal part of the curl of a vector field 

vanishes. For V»(v x F) * 0, and hence (V x F) r » 0, as seen from (13). Also 

h 

from (13) we see that a solenotdal vector field (i.e. a field whose divergence 
vanishes) is purely transverse because its longitudinal part vanishes. And 
(25) shows that an irrotational vector field (i.e. a field whose curl vanishes) 
is purely longitudinal because its transverse part vanishes. Vice versa: A 
purely longitudinal vector field is irrotational, since - by ( 13 ) - it is the 
negative gradient of a scalar and V x (v<|>) = 0. And a purely transverse 
vector field is solenoidal, since - by (25) - it is the curl of a vector C 
and V«(V x c) = 0. Thus, the properties ’’purely longitudinal” and ’’irrota- 
tional'’ are equivalent, as are the properties "purely transverse” and 
"solenoidal". A non-vanishing vector field F cannot be both purely longi- 
dinal (meaning F T * 0) and purely transverse (meaning F^ = 0) at the same 
time, since F * F^ + F^, yields F = 0 under these conditions. (Of course, 
the assumption is always that F tends to zero for Jrj 00 rapidly enough so 
that the Fourier transform theory, on which everything was based, is appli- 
cable.) Therefore, because of the equivalences of properties spelled out 
above, we conclude that a field that is both irrotational and purely trans- 
verse is zero everywhere. The same conclusion holds for a field that is both 
solenoidal and purely longitudinal. 
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We apply the results of this section to the Maxwell equations of the 
electromagnetic field. 


? ” E + jft c8 * 0 * 

(29a) 

? K cS - sit E ■ n ° 3 • 

(29b) 

V • e « n Q cp , 

(29c) 

$ • c5' a 0 , 

(29d) 


Here 


E * Electric field strength , 

B » Magnetic induction , 

n = Impedance of free space * 3TT ohm. 
o 

The quantity n Q is related to 

_ _ 7 volt sec 

V = Permeability of free space = 4tt 10”' ■ ■ — * 

o amp m 


c = Speed of light * 3 x 10 


f; JSL 

sec 


= Permittivity of free space = , 

o 


through 



I _ y Q c - E 
o o 


Of course, the electric charge density p and the electric current density J 


’must be coupled by the continuity equation 

5-1 + sit Cp = 0 • 


(30) 


We apply (29) in the following manner: cp(r,t) and J(r,ct) are given; 
E(r,ct) and cB(r,ct) are to be found. Equation (29a) tells us that cB is 
solenoidal, hence purely transverse. Thus 


cB = cB T • 


(31) 


Then, by (24) and (25b there exists a purely transverse vector field cA,^ 


such that. 
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cB * V * cA^j, . (32) 

The vector c/L, is called the "vector potential in the Coulomb gauge" . We 
decompose E and also J into their longitudinal and transverse parts, so that 

E * e l + e t , (33) 

J » J L + J T . (3*0 

Then the equations (29) become (Note that V * E^ * 0 and V*E T * 0) 


"'Vs ' 5 ! 1,0 (35a) 

? * (5 x e V * aft e l - “It ^ * Vl + Vt (35b) 

$»E t ■ n cp (35c) 

Li O 

?•(? x c 5 t ) = 0 (35d) 


Equation (35d) is redundant. Equation (35a) tells us that E^ + cA^, 
which is obviously purely transverse, is also solenoidal, hence purely 
longitudinal. Therefore E^ + cA^,, being both purely transverse and 

purely longitudinal, must be zero everywhere, so that 


®T = ~ act cS T * ^ 

We no longer need equations (35a) and (35d). Equation (35b) must hold 
separately for its transverse and longitudinal part. Thus, with (36), 
this equation yields 

V * (V X cJ,) + eJj - n 0 J T , (37a) 


and 

- — E_ - n J T . 

8ct L o L 

We also have [see (35c)] 

V»E = n cp . 

L O 


(3Tb) 

(37c) 
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Now, by (13), there exists a scalar field <fr, such that 

E l - - . (3«) 

The quantity $ is called the "scalar potential in the Coulomb gauge". Then 
(3Tb) and (37c) become 

? sit * • <39a) 

and 
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- v»v$ ■ n o cp . 


(39b) 


Equation (39b) states that at any instant ct, the scalar potential is equal 
to thn electrostatic potential associated wiuh the charge density at that 
instant. But one should keep in mind that $ depends on the time, because cp 
depends on the time. 

The procedure for obtaining the electromagnetic field from a given 
charge density cp(r,ct) and current density j(r,ct) is then as follows. 

First one determines the potentials <j>(r, ct) and cA^r.ct) from [See (37a) 
and (39b)]. 




and 


- q o c$ , 


9 2 


5 x (5 * °V * 15^7* ’ Vt 


Here, n may be obtained either by way of $, namely by [See (39a)]. 

V\ ■ V - ■> 


(40a) 


(40b) 


(kOc) 


or one may bypass the determination of $ and obtain from J itself through 
(l6) or (25). After 4> and cA^, have been determined, one obtains the fields 
E(r,ct) and cB(r,ct) from 

(Lla) 
( Uib ) 


E = E^ + E^ , 


E L = - V* , 
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3 

B t ■ " act cS r • (Ulc) 

cB * V x cS^ . (4ld) 

Equations ( UO ) and (4l) are, of course, well-known. We showed how 
they are derived, not only for the sake of completeness, hut also because 
we wanted to exhibit why the source function for cA^, i.e. the right-hand 
side of (40b) is not 

One of the favorite ways of dealing with equations (40) and (4l) is by 
way of their spatial Fourier transforms. What corresponds to the operation 


Vin the r-space is the multiplication by i»c in the ic-space. 

Thus the trans- 

forms of and $ x ($ x cA^,) are 


ix • iic<{> * - k 2 <£ 


and 


* A A 

iK X (iK x CA.J,) * K 2 cA T - wc.clj 


A 

But since V*cA^, * 0, we have itc-cA^ * 0, so that 


A A 

iK X (i< X cA,j,) ■ K 2 cA,p . 


Thus, equations (40) and (4l) become 


A A 

K 2 (j)(K,ct) = n Q Cp(K,Ct) , 

(42a) 

K z cA tp (K,ct) + cA^K.ct) * n o J T (K,ct) , 

(42b) 

n o ? T <K,Ct) » n o J(K,Ct) - iK-^-^(K,Ct) , 

(42c) 

and 


A A A 

E(k , ct ) = E^(k , ct ) + E T (K,ct) , 

(43a) 

A A 

Ej^(<,ct ) = - iK?(K,Ct) , 

(43b) 

A 0 A 

^(K.Ct) = - ■^CA <r (K,Ct) , 

(43c) 

A • A 

C§(K,Ct)*iK x cA^.Ct) . 

(43d) 
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Equation (42b) states that the Fourier transform c5(E,ct) behaves like the 
excursion of a driven harmonic oscillator, whose resonance frequency flt'ic) is 
given by 

(£)2<k) . . (Mt) 

A 

The "driving force" for such on oscillator is n^^^ct). We shall return 
to this remark in the next section. 


The energy density P energ y th* electromagnetic field is given by 


p • a x" " '* " (E*E + cB*cS) 

"energy 2n e 


or, since E * Ej. + E^ , 


P = (E T *E r + 2E'E rn 4- E»E_ + cB-cS) 

"energy 2n Q C L L L T T T 


(45) 


(46) 


When we integrate this expression over the entire r*space we obtain the stored 
electromagnetic energy U, According to the orthogonality relation (27), the 


integral of E^'E^ vanishes. Thus 


U=U..+U , , 

stat rad 


with 


(47a) 

(47b) 

(47c) 


U stat * 2^ , 

U rad " 2^ /‘W 5 ! + ’ 

o 

■ *h °v + x **!>•(*■ * 

o 

The energy U is therefore the sum of two parts, the electrostatic energy 

associated with the longitudinal field E_ and the energy U , of what may be 
the h raa 

calledyradiation field, i.e., the field specified by cB and the transverse 
part of E, The decomposition (47a) is a natural one. We shall be inter- 
ested, mainly in U rad » since we can regard U gtat as a purely configurational 
energy. For it depends, by way of <f>, only on the, instantaneous charge dis- 


tribution. 
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We can convert the integrals over r-space in (4?) into integrals over 
tc-space by means of the Parseval relation for Fourier transforms. This 
relation states that for two complex vector fields F(r,t) and 3(r,t) ve have 

L " 5 /» . % 1 t * /■> ■ . l" * . t ( 1^8 J 


/dr r (P(r,ct) • G(r,ct) ■ ^ g - yj /Ax K (P(x,ct ) *S(K,ct) 


We use 


•"'X. 


A # 


(V K CA T ) • (v X cA^) * (-iK * (cA^p) >(i< ?< cA t ) 


A * 

- \ " 


4 .# 


* ^{cAjp) • cA T -k“*(c'A t ) K'cAp * K 2 (cI T ) »cA - 0 . 

A 

(Note that rcA^ ■ 0, since V 'c^ = 0). Then (47) becomes 


U ■ U stat + U rad' 


with 


U 


'stat * 5Te /*V Z(f t‘’ ot ) * ( *> ot) ’ 

o 

U rad * 2ST W'slt «V 5 ' et)1 * sit c V“' ct) + 


(1.9a) 


(49b) 


+ K 2 (cA r (^ct)*.cA T (^ct)] . (4?c) 

Equation (49c) states that U rad is the sum (actually an integral) of the 
energies of the individual field oscillators, whose excursions are given by 

A 

cA f j,(’ic,ct). Seemingly, each oscillator has three degrees of freedom, since 

A 

cA^KjCt) is a vector. But since \his vector is constrained by the condition 
that it be normal to the wave number vector k, each field oscillator has 
only two degrees of freedom. 

We are now ready to show (in the next section) that free electrons, 
even though they may have beats, cannot radiate. 

Ill . Proof of the Statement that "Quantum-mechanical beats cannot enable 
free electrons to* radiate in free space* 1 . 

We start with the expression (2, 49c) for the energy U rad of the 

radiation field. We repeat it here for the sake of convenience. 

U rad (ot) * 5T? ^ T[( 3rt °V V t) *' alt ♦ 

OK 

+ K2(c5 T (K,ct) 4t .(cA T (K,Ct)] . 


( 1 ) 
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We assume that in the distant past, ct - *, there was no electromagnetic 

field present as yet. Then what we wish to show is that in the distant 

future, ct •+ + ®, the radiation energy is zero, i.e., we wish to prove that 

lim U .(ct) * 0 . (2) 

ct •*> +« * ra 

In the meantime, there may be a build-up of radiation energy. But this 
energy will be reabsorbed again (not radiated away) in the end. 

A 

• We determine cA^KjCt) by means of equation (2-42b), which we repeat 
here for the sake of convenience, 

' fact) ' 2 + tc 2 cAfp( K,ct ) = n 0 ? T (K,ct) . (3) 

As we said before, this is the equation for a driven harmonic oscillator. 

A _ 

We determine the driving ’’force” n o J^(ic,ct) by taking the spatial Fourier 
transform of (see (1-42)) 

J(r ,ct) = (-ec) (-ii|/*V^ + iW*) . (4) 

(Afterwards, we extract the transverse part.) 

In order to do this, we start with (1,6), but use a different dummy variable, , 
namely I. Thus • 

iKr,ct) * /AT ? F(I)exp(il-r - 1 ct ^ * (5) 

Similarly, this time with the dummy variable n, 

V (r,ct ) = /Ax^F*(n)exp(-in*r + i ct) . (6) 

On taking the gradient of (5) and (6), we obtain 

(ViJ>) (r ,ct ) = j—p- !h t/(il>F(i)exp(i?-r - i ~ ^ ct) , (T) 


me , 2 
S 1 

yyy n ^ 


(V\p )(r,ct) = U (2„) J /^ T n (-in) P (n)exp(-in-f *• i ~ ^2 ct}* (8) 
We insert the expressions (5) - (8) into (4). But we write the product of 
two integrals as the double integral of the products of the integrands. We 


obtain 
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J(r,ct ) » (-ec) ~ // Ax g AT n F(g)F*(n) 

exp(i(£-n)*r - i %U 2 - n 2 )ct) . (9) 

me 

Instead of 5 and n * we use two new dummy variables If and T defined by 

< * Z - n , \ (10) 

X * ^ J, 

Vice versa: 

I * + + 1 , 

I 

n * - Jjic + X . 

The Jacobian of the transformation (ll) is equal to unity, as one can see 
when one examines the three Jacobians for each of the three cartesian 
components; all three are equal to unity, and 1 ■ 1. Thus At^At^ * At^At^, 
and equation (9) becomes 

si 
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(11) 


J(?,ct) = (-ec) II A T< AT x e ik,r THX «)F* tf - h<) . 


exp(-i ~ <• Act ) . 
me 


( 12 ) 


From this equation one can read out tL<* spatial Fourier transform J(tc,ct). 

See (2,2).) It is given by 

J(ic,ct) - (-ec) ~ 7’ g ' ^y3 jAr^a + %k)F ( A -^x)exp(-i j~Ic»Act) . (13) 

A 

However, in the equation of motion (3), we need J T (ic,ct), the Fourier trans- 
form of the transverse part J T of J. The definitions (2,16) and (2,17) relate 

A A 

J,j, to J. Thus 

J T (<,ct) - (-ec) /At a (T - kk«X)F(J + J<?k)F*(X - k<) 


and 


^ _ 

exp(-i - — <•! ct), for k ^ 0 , 
me * * 


(lfc) 


J T (3:,ct) = |<-ec) /ATa^F(I)F*(A) . (15) 

Ufix) 2 

We now make the assumption that the energy — ^ — associated with each 

A 

Fourier component \|/(K,ct) of the wave function i|i(r,ct) is limited to a value 
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that la much smaller than the rest energy me a . This assumption la fulfilled 
In the device that Bolivars proposed, Furthermore , wo have to adopt this 


assumption already in view of the decision, made at the outset, to use a non* 
relativistic treatment* This assumption implies that F(v) he aero whenever if 
exceeds some maximum value xv*> which is chosen in such a way that the correa-* 

y* m 

ponding energy , reckoned non-r e 1st ivi a tic ally , is much smaller than 


mc a . for instance* we might choose w... to he equal to t ^4* * so that 

KP , ' ™ ■■ 

— —in — * ^ mc ? , which is in the order of '&> kllo^eleotronvolt* We summarise 


Pm 


these remarks by the statement that 


P(x) * 0, for w > w * (lb) 

One consequence of (lb) is that 

*T;j,(w,ct) a 0, for w >■ » (It) 

ftie reason is that f(t + yw)F (t - -\k) in (ltd is aero, unless X lies simul- 
taneously in the two spheres (drawn in X-spaoe) jX *» and |X * (~ 

(Their centers are at y5 and - Vw respectively*, and both have the radius x » ) 
But when |Vr| v k , tlie two spheres do not overlap* thus the integrand in (It) 
is aero for every X, and (it) results * 

A more important consequence of (lb) is that, for each i, the quantity 


d,„(K,ct) is band-limited to a low-pass frequency band* In order to see this, 


m » M n , A that occur in the time factor 


we examine the frequencies *p « * ~~f x 

Vj. lUv 

expt-4 x*X ct) of (It). Then, what we are looking for, are the extreme 
values of x*X, when t Is restricted to lie in the domain that is common to 


the two spheres IX - h kv and IX * I <x , A sketch of the two spheres 
will show that the extreme values of x*X occurs when X lies at the two points 
1 (<* Vx * x )n aw) * (The surface of each sphere intersects the axis that goes 
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through - kx and in tvo points, whose distances from the origin are ^ 4 % 
and ^ - k*' The two critical vectors J that we need lie at the inter- 
section points that are closest to the origin. And th&i gives us the two 

points 4 ( -Yk 4 k ~ k)]. Therefore, for a given wave vector tc, the fre- 
m k 

quencies — ■ k»a that occur in the excitation function n J m (K,ct) of (3) 

c me o x 

ore limited to the band 

- K £ ( K ->(k) + -W . (18) 

me m — c — me m 

A 

The excitation function n Q j(K,ct) is thus band-limited to the low-pass band (l8) 
We could substantiate this conclusion by a detailed calculation of the 

A 

temporal Fourier transform r\ ^{k,— ) (pronounced J ay-tee-hat-tilde) of the 

OX c 

A 

excitation function J^K.ct). This transform is defined by 
Z -i <“■ ct A 

* / " dot e c J-<Tc,at) . (19) 

ct=-»> 

The inversion formula of ( 19 ) i3 

i ^ ct 

J T (K,ct) 3 2 “ / d “ e c J t ( k,|) (20) 

(1) 

— 

C 

However we do not need to perform this straight- forward calculation (which is 
done most convincingly with (20) as the starting point)* The reason is that, 
as we shall see later, we need J (k,~) only for one particular value of the 

X C 

frequency — . And this value corresponds to the resonance frequency k of the 

G 

harmonic oscillator 1 type equation of motion (3). Now, k lies well outside 
the low-pass band (18), (For our numerical example K m e » ihe resonance 

frequency k lies outside the low-pass band (lO) by a margin that is better 


than 2 to 1.) Thus we arrive at the important equations. 

V*’** * 0 * \ 1 — 1 . 

3* f for |k| f 0 . 

J, r (fc,-K) = 0 ,) 


(21a) 
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Now we take up the roason that we are interested mainly in the temporal 

^ ^ (a) 

Fourier transform n ?,„{£,* ) at the single frequency value — » < , rather than 

O l c 

A ^ 

in the complete time function J T (w,ct). We assume that the electromagnetic 
field was in its quiescent state in the distant past, so that 


lim cA,_(r,ct) * 0, lim cA_(r,ct) * 0, lim U ,(ct) * 0 . * (22) 

ct-*-« 1 ct+-« dcrc ct-+— °° raa 

The first two equations of (22) imply that similar equations hold also for 

the spatial Fourier transforms, so that 

lira cL(*,ct) = 0, lim ~r cl (<,ct) = 0 . (23) 

ct-Hi 1 ct-*-® 1 

We wish to evaluate how much radiation energy will have been produced in the 

distant future, i.e., we wish to calculate lim U ,(ct). We start with the 

ct"H-°° a 

expression (l) for U ,(ct). We then need lim cA m C<,ct) and lim -r~r cA-CKjCt) 

ct*H-°» T ct-M-® 9ct ‘ ^ 

We could determine cA^Ci^et) directly from the equation of motion (3). However, 

A 

we find it more convenient to determine the two auxiliary functions F(T?,ct) 

A 

and GCic,ct) which are defined by 


F(;,ct) = (-jJJ- cA T (ij,ct) - iKcA^t K,ct))e irct , 


(24) 


A a A A ___ 

G(’K,ct) = cA T ('K,ct) + i»C3A T (K,ct))e' 


-i rot 


• w 

where, as always, « = | . Note that both auxiliary functions vanish in the 


distant past because of (23). For the complex conjugates, we have 


F (7:,ct) = (~r cA^(K,ct) + ik cA T ( K,ct)^“ lK “ t 

§ # (K,ct) = (:— cA^dc.ct) - itocA^CkjCt^e^'^ . 
From ( 2U ) and (25) we obtain 

&($ (lc,ct).F(K,ct) + G (Ic,ct)*G(x,ct))=s 


(25) 


3 ^ 9 A ^ a 

= ^3ct cy^.ct ) +k 2 cA t ( <,ct ) .cA t (k ,ct ) 


(26) 


One should not confuse the (vectorial and time-dependent) F in (24) with the 
the (scalar and time-independent) F in (5). 

^What makes it possible to stipulate these initial conditions is that, in the 
distant past, the beats have not yet developed and that the wave packet is of 
infinite size. 


CWWWWl Wttft 
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The right-hand side of (26) is the integrand of equation (l). Thus, (l) 
may be written in the form 

U rad^ ct - = T+tmT I 6 k [F*(<»ct )• f (E,ct ) + S*(<,ct )• G(ic,ct) ] . (27) 

o 

A A 

Now we calculate the two. auxiliary functions F^Jc,ct) and G(<,ct). We 
take the time derivative of (24) and use the equation of motion (3). Then 
(there are two terms that cancel) 


3 £ _ / 3* r ..9? \ ikct 5 ii<ct 

15t F = ( IW CA T + K cA T )e " n o J T e 


let! 

32 


( 28 ) 


3s / t . 9 4 x -ikct . % -ixet 

art G = { Ttet7F cA t + • c V e = Vt* 

A A 

The functions F and G can then be determined by simple quadratures. Because 

these functions vanish for ct •+•-«, we have 

F(k,ct) = / ct dcs e'* KCS n J_(k,cs) , 

cs=— « °. T (29) 

G(k,ct) = / ct dcs q J t (k,cs) . 

CS=— 00 

Now we take the limit ct -► +». In view of the definition (19) for the 
temporal Fourier transform, we have 

A A .... 

lim F(K,ct) = n J.(k,-0 , 

Ct-H-® . 

A A 

lim G(K,Ct) = n 0 ?m(K,K). . 

Ct -H-« 

We use these expressions in equation (27). Then the radiation energy in 
the distant future becomes 

I*™ u rad (ct )=■!“. / at k [(J t ( <,-<))*• j T (K,-K) + 

+ (J t (k,k) «J t (k,k)] . (31) 

But according to (21), all the Fourier transforms in (31) are zero. Thus 


(30) 


lim U (ct) = 0 . 
ct~w rad 


( 32 ) 


we have the final result 
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In words: A single free electron in free space, even though it may have 
quantum-mechanical heats, cannot produce any radiation energy in the end. 

Of course, there may he some temporary production of radiation energy. But, 
in the end, this energy is reabsorbed by the electron. What made it impos- 
sible to generate radiation is that, for a heat of wave number ic, the 
frequency w/c is well below the required value k = |<|, And the reason for 
that was the hand width limitation (l8). 

So far we have considered Just one electron. Now let us consider a 
swarm of electrons. We enumerate the electrons with some digit p * 1,2,3, ... 
Let p7(r,ct) he the current density that belongs to electron #p. When we 
wish to adapt our previous calculation to the swarm of electrons, then we 

merely have to replace 7 hy the sum \ 7 . In equation (30) and (31) there 

-r _ _ p P 

appear the sums £ 7(1c,-tc) and J, J(k,k). Since each term in the sums is zero, 
P P P P 

the sums themselves are zero. Thus, equation (32) holds also for. a swarm of 

free electrons in free space. Note that we did not make any assumption about 

the instants at which the electrons are emitted. The conclusion (32) holds 

not only when the electrons are emitted at random instants, hut also when the 

electrons are hunched hy some gating device. 

This finishes the proof for the absence of radiation. We conclude this 

section with some miscellaneous comments. 

Several times we have mentioned the spatial and temporal Fourier trans- 

form J^CJc, “) of the transverse current density J^,(r,ct). We never needed to 

'7C 'X 

calculate it, because only the special values J,j(k,.k) and J^Ck,-*) were 

required in our proof. And these values were seen to he zero. Nevertheless, 

for the sake of completeness, we exhibit an expression for the general value 
» (0 

J,p(*, — ). We start with equation (lU) for the spatial Fourier transform 
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Jfj,(K t ct). We repeat this equation for the sake of convenience. 

J T (i?,ct) ■ (-ec) ~~ |2 tt) 3 + &t)F*(X - 1&) 

/ 

exp(- i ~ ? »X ct) (lh, repeated) 

me 

In the X-space, over which the integration in (lU) is performed, we introduce 
cartesian coordinates 5, ti, r, referred to the three mutually orthogonal unit 
vectors ~ , a, b. The first of these depends on the given wave number vector 

Is* % 

k , which appears qn the left-hand side of (lU). Then we can always choose the 
remaining unit vectors a and b in such a way that a, b ore mutually ortho- 
gonal. Thus we can write 

X * £a 4 Tib 4 C ~ 

1C 

and 


* d^dqdj;. 

Equation (lH) then becomes 

'fi* 1 tt 


J„(i<*ct) - (-ec) 


T 


HI 00 dedndcUa + n«)F(Ea 4 nb 4 u + |) M 
5, n, 


P ( £a 4- qb 4- (c - —)£•) exp(-i ~ Cxct ) . 

c k me 


(33) 


0) 


Here we replace the integration variable c by the integration variable — » 

C 


such that 


Then 


co 

c 


jf i 
me 


kK . 


, „ me 1 , to 

dC = — t — d — , 

*il K C 

and ^ runs from 4® to -®. We can remove the minus sign that results from 

the equation for dc by making — run from -« to 4® again. Then equation (33) 

c 


becomes 
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^(K.ct) * (-ec) {“J2”yy ///*** d jfdedr<Ca+nE)F(Ca*itf> + (“£ “■ “ + |’) i |) 

5 ,n,to/c *-« 

F ( 5a+nE+(^ ^ - |)“) «xp(i ~ ct) . (34) 

A 

Next we note that the spatial Fourier transform Jg,(K,ct) is related to the 

'X _ " " 

spatial and temporal Fourier transform “) by 

J T («,ct) ■ ^ /*" d “ J t (k, ^)exp(i ~ ct). (35) 

to 

— s_co 

C 

On comparing (34) and (35) we obtain the final result 

5 t (k, f) 1 (-ec)^ j' g ^ 5 // + ” d£dn(£a+nb)F(£a*nb + ^ “ + f)* ) 

p't^db ♦ <f ia.|)|). ( 36 ) 

Here, as always in our calculations, the quantity tc stands for |k|. On the 

other hand, ~ takes on both positive and negative values. 

According to our assumption about the energies encountered in the wave 

packet, the functions F and F in (36) vanish unless their vectorial arguments 

are smaller in magnitude than some maximum value tc , for which we suggested, 

m 

by way of an example, the quantity ^ ^ . Thus the effective domain of the 
( 5 ,n) - integration is constrained by the two inequalities 

r 2.„2.f' Inc 1 w . k \2 2 

(3T) 

ii tc c 2 m 

Therefore, the effective domain of the integral in (36) is finite, in spite 
of the 4j» marked on the integration signs. Hence, there are no convergences 
problems . 

Again, let us take the example k = — ^ . Let us now calculate J,j(ic., for 

— = ± k. Then, it cannot ever happen that both inequalities (37) are satisfied, 

c 
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(For — = the first is definitely violated, and for — * -k, the second is 
c c 

definitely violated). Thus the effective domain for the integration in (36) 
has shrunk to zero, so that 

? t (k»k) = o, 


for ic i 0 . 


= 0 , 


(38) 


This result agrees with the equations (21a, h), which were established without 
a detailed calculation of the spatial and temporal Fourier transform. 

We may establish the validity of equations (38)* on which the impossibility 
of radiation was based, also directly from equation (9). We modify this 
equation so that it yields the transverse part J,p. The modified equation is 


f(£)f*(h) 


exf [i(l-n) »r - i ~ U 2 -n 2 )et] . ( 39 ) 

me 

We generalize this equation to a form which is valid also in the relativistic 

region of velocities. But we still adhere to the description in terms of 

scalar Schrddinger waves, i.e., we refrain from using the Dirac spinors. The 

relativistic formulation allows us to dispense with the introduction of some 

maximum wave number fc , as we did in connection with equations ( 16 ) and (IT). 

m 

In order to obtain the relativistic generalization we have to replace the 

quantities — and — kq 2 in the exponent of (l4) by A^) + 5 2 and 
__ me me v 

+ r > 2 * These replacements follow from the relativistic relation 

E 2 = (me 2 ) 2 + (cp ) 2 (Lo) 

between the energy E and the momentum p, as well as from thede Broglie relations 

E = die ~ , p =/ 6 f, p = din . (hi) 

c c, n 
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The relativistic generalization of (39) is then 

4C 

T V- AT J JQT.CT_ L“TT“ “ „ 

U-nMI-n) 


J_(r,ct) « (-ec) ~T Tp I T S ' //AtJ*t . I|i] 

T me (2 tt)6 5 H 2 ( .:),( -) 2 


FU)F (n) 

exp[i(c-n) • r - i(>/^') 2 + + n 2 )ct] . 


(42) 


This equation shows that each doublet (C,ri) of wave number vectors contributes 
a sinusoidal wave to if^. The wave number vector ic and the frequency ^ of such 
an elementary wave are given by 

k = S - n t 

? -JW^-JW 2 ^ 2 ) ■ 

(Different doublets (!,n) can produce a common ic and 


(43) 


U) 


But this lack of 


a one-to-one correspondence need not concern us.) Now let us examine the 
ratio of an elementary wave. We have 


U/cl _ 
IkI 


M )2 ** 2 


'5 - n 


But since 


we obtain 


l - n i 


V - n 


/F - /tv* 


kM < 


W ** 2 -V<f> 2 +n 2 | 
/F - *F 


Now 




/Pf) 2 ** 2 + 

/(f) 2 + n 2 


5 2 -n 2 

vF"- *F 







so that 



original page is M3h 
0F "COR quality 


U/cl < 


{jj? ± /r? 


*i 

The expression on the right-hand side is clearly less than unity. Thus we 
arrive at the result 

(W) 

1*1 

The inequality (44) tells us that every elementary wave in (42 ) has a fre- 
quency-t o-wave number ratio that is less than unity. But J^(lic,k) and J t (k,k) 
are the amplitude factors of elementary waves whose frequency-t o-wave number 
ratio is equal to unity. Since we have seen that elementary waves with a 
ratio equal to unity do not occur in (4 2), we conclude that 

^p(M) * 0 , ) ... , 

: ^ for (tc| 4 0 . 

J t (k»~'0 = 0 *) (45) 

And these equations are identical with (38). 

We conclude this section with a remark about the spirit of the calcula- 
tions we have performed. We started with a given free wave packet of 
Schrtfdinger waves. The word "free" indicates that the wave packet is not 
influenced by electromagnetic forces. Then we calculated the electromagnetic 
field that is generated by the electric current density associated with the 
wave packet. Now in principle, this field would influence the behavior of 
the SchrSdinger waves. However, we neglected to take into account this back- 
reaction that proceeds from the electro-magnetic field to the SchrUdinger 
waves. We only considered the forward action that proceeds from the 
SchriJdinger waves to the electromagnetic field. What we have done then is 
in accordance with the perturbation calculus, whose principles are explained 
in Section 10 of the Appendix. We had to resort to the perturbation calculua 
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■because an exact calculation 13 impracticable. It is believed that the 
accuracy of the perturbation calculus is good enough for our purposes. One 
could, in principle, adduce some evidence for this belief by carrying out 
the perturbation calculus to the next higher order of approximation. But 
this would entail a lot of effort, more than we can afford to devote to 
this task. As far as this writer knows, there are no general principles 
that could obviate the necessity of those cumbersome calculations. Even if 
the next order of approximation were to predict some radiation, it would be 
too weak to be of practical use. 

We have seen that no radiation can be expected from free electrons, even 
if they exhibit beats. Naturally, there comes up the following question* 

Can one use the beats in some other way? For instance, one may ask whether 
electrons with beats interact .with externally produced radiation in some 
peculiar, and perhaps useful, way. We explore this question in the next 
section. 

TV. A Selective Moving Mirror , 

H. Schwarz suggested a device that produces an electron beam which is 

a superposition of two plane waves of different wavelength. In mathematical 

language: The SchrBdinger wave that describes the electron beam is of the form 

tj/(z,ct) - C{exp(i'Kiz - i ~ -Jr ct ) + exp(iK 2 z - i ct )} . (l) 

Here C is an amplitude factor. The beam travels in the z-direction of some 

cartesian coordinate system. Each of the two partial waves is a plane wave. 

The first partial wave has the wave number and the kinetic energy ~ Kj . 

in df K 2 

The frequency Wj is therefore given by ^ -J- . Similarly, *c 2 and * 2 2 

c me d 2m ■ 
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are the wave number and kinetic energy of tfo* second partial wave. Both 
kinetic energies are of the order of 30 x 10 3 electrovolt. And their 
difference is in the optical range, i.e. we have^^ ( K j 2 - *i 2 ) * 3 electron- 
volt. The constants df, m, c have their usual meaning, namely ~ * Planck’s 
constant, electronic mass, and speed of light in vacuo. 

We give a brief description of Schwarz’s design. An electron beam is 
sent through a beam splitter that splits it into two divergent beams #1 and 
#2, Several inches down their paths, they are made to converge again into 
the working region of the apparatus (near z * 0), where they form a single 
beam again. This writer is not competent in the field of electron optics. 
Therefore, we will refrain from describing how this may be accomplished. 
Perhaps one may imagine that the beam splitter and, further down the line; 
the beam joiner consist of diffraction .gratings, presumably in the form of 
single crystals. What is important to keep in mind is that the beam splitter 
does not sort out different individual electrons. Instead, it splits the 
SchrBdinger wave of each electron into two divergent portions. Similar re- 
marks may be made about the beam Joiner. So far then, the lay-out of the 
beams is as follows. First there is a straight and narrow beam along the 
z-axis. The beam enters the beam splitter and becomes two beams, diverging 
off to the north and the south respectively. After a few inches down the 
line, these two beams are bent back to the z-axis. Where they meet, they are 
Joined by the beam Joiner and form again a straight beam along the z-axis. 
(The beam joiner Will waste some of the beam by sending off additional beams 
to the sides. ) Eventually this beam enters the working region of the 
apparatus near z = 0, 
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In the region between the two curving beans , a device is installed that 
contains a magnetic field in the vertical direction. This field is made 
to change with time. According to Faraday’s induction law, an electric 
vortex field arises that speeds up the northern partial beam and slows down 
the southern partial beam. The beam that emergies beyond the beam joiner 
is then a superposition of two beams with different energies , the energy 
difference being chosen to be in the optical range. Because the beam splitter 
and the beam Joiner act on the Schrtfdinger wave of each electron, the super- 
position is linear as described by equation (l). 

If, in a more realistic manner, we consider that the Schrbdinger wave of 
an individual electron would be a localized wave packet instead of a c.w. 
wave (c.w. = continuous wave) of infinite extent, then we conclude that an 
additional device is needed. After all, since the northern beam was sped up 
whereas the southern beam was slowed down, the southern wave packet would lag 
behind the northern wave packet beyond the beam Joiner, so that the superposi- 
tion described by equation (l) could ( not take place. Therefore, wo have to 
install a delay section in the northern partial beam to remove the lag. 

From now on, we shall assume that a pe 'feet design of the electron beam 
in the working region has been achieved. By this we mean that the Schrddinger 
wave <Hz,ct) is the linear superpostion of two partial waves i|»i(z,ct) and 
MZjCt), 

ij>(z,t) = t|>,(z,t) + \p 2 (z,t) 


9 


( 2 ) 
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and that at the center of the working region (i.e, at z * 0) the quantities 
^ 2 ( 0 , ct) and i|>j(o,ct) are related by 

i|i 2 (o,ct) = i^ 1 (o,ct)e" iFc *' * (3) 

where F is a constant. In terms of k, 2 and Ki we have 

F -^% ( k . 2 2 - ^2) . (It) 

Of course, -ficF is the energy difference for the two partial beams. We 
assume that (3) holds, no matter how i|Ji(o,ct) depends on the time t. As 
we implied by the notation >J>(r«,ct), we ignore any dependence on x and y. That 
is, we assume that the SchrSdinger waves are plane waves. Certainly the c.w. 
wave (l) conforms to the specification (2) - (4). The specification (3) 
permits the investigation nqt only of c.w. waves but also of wave packets. 

To begin with, we examine c.w. waves. We will discuss wave packets later 

on. We start with the example of equation (l). The density of electrons is 

# 

given by i|/ i|>. In our case, we obtain 

= 2C*CU + cos[(k. 2 -Ki)(z - £*|!Sl Jl St)]} (5) 

The first term in the curled bracket describes a constant background, whereas 
the second term exhibits the beats. The beats persist with constant strength 
for all values of z. The beats have the wave number K . 2 - K ‘i» and their 
velocity v is given by 


V _ k) +IC.9 

c “ me 2 


( 6 ) 


Now we come to the matter of wave packets . Let us first consider the 
partial wave #1, Instead of the c.w. wave 

( 7 ) 

we consider a superposition of plane waves, each with a different wave num- 
ber < . However the ic-values should be closely clustered around the nominal 


i|>l(z,ct) = r; exp[i<iz - i ^ki 2 — ct] , 
■ me 
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value K\, Accordingly we write 

^l(z,ct) * Mj / +# ° d* f (k ) expfiKz - i k* 2 “ ct] , 
**-» roc 


( 0 ) 


where ]f(w)j z is sharply peaked near k = jcj. The quantity Nj is a (complex) 
normalization factor, which we determine later. The conclusions that Wj wish 
to draw do not depend on the choice of the weighting function f(»t), as long 
as it is sharply peaked. We choose f(«) to be equal to unity when k lies in 
a narrow range centered at Kj, and we set f(*c) equal to zero when * lies 
outside this narrow range. We therefore write 

'h(Zft) = Nx / <1+0tl d* exp[i<z - i ^Sk 2 “ ct) , (9) 

We demand that 


oil « K2 - K1 , (10) 

where K‘x and *<2 (<2 * <i) are the nominal wave numbers of the first and 
secnd partial beam respectively. We replace the integration variable k by 
? = K-Ki. Then (9) becomes 

^j(z,ct) « Nx /**! d? exp[iU*x+?)z - I^kx^kxS+S 2 ) ct] . (10a) 

5=-<M mC 

So that we may evaluate the integral in a convenient closed form, we decide 

to keep track of the wave packet only for those times that satisfy the 

condition 

«J 2 ~ |ot| « 1 . (11) 

Then we may delete the term ? z in (10a) and obtain 

||» 

(|) 1 (z,ct) = Nx exp[ Ik x z - ik Kx Z ct] / fll d? exp[,i(z-Kx ct)?] 

ttil. m iUv 

?=-«x 


or 
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^ 2sin[aj(z-Ki ct)] 

»|» J ( z ,ct ) ■ Nj exp[iKjz - 1 h<\ 2 Z“ Ct] - 


( 12 ) 


We might express the approximation we made when we deleted the ? 2 -term in a 


different way: We replace the exact parabolic relation 

a = H** 

c me 


(13) 


to 


= %1 2 ~ * <1 (tC-Kl) * 


(14) 


between frequency — and wave number k by its tangent line at Kj , i.e. by 

w 

£ 

me * me 

After having made this replacement we can afford to abandon the inequality (11) 
Thus the substitution of the new "dispersion relation" (l4) for the old re- 
lation (13) makes the wave packet formula (l2) an exact one. 

The exponential in (12) describes the ripples of the wave packet. The 
last factor in (12) describes the envelope. The envelope consists of a 

central peak with its center at z^ - — ct. The height of the peak is 2cq 

2it 

and its width is — . On both sides of the central peak, the envelope factor 
^ 1 

decays, with ups and downs, to zero. The envelope pattern moves rigidly with 
the velocity v* given by 


= K 




(15) 


c * me 

The rigidity of the envelope pattern is a consequence of our approximation 
('(l4) instead of ( 13 ) ] • An exact calculation would show that the peak becomes 
wider and lower for times that do not satisify the inequality ( 13. ) . But this 
is of no concern to us. 

The expression (12) describes a wave packet whose center passes the 
origin z = 0 at the time zero. If the center passes z = 0 at some other 
time s, then we merely need to replace ct by ct - cs. We are allowed to 
multiply by the constant phase factor exp[-i Jgkj 2 “ cs], or - expressed 
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y 

differently - to absorb exp[+ i k'Ki 2 cs] into the normalization factor Nj. 

IIlw 


Thus, for a wave packet whose center passes z * 0 at time s, we have 

^ 2sin[ai(z-Ki (ct-cs))] 
<|»l(z,ct) « Ni exptiKjz-i H<\ z —7 ct] ,!lc 


me 


z-Ki — (ct-cs) 
me 


( 16 ) 


The concomitant wave packet in the partial beam ff 2 is described by an 
analogous expression, namely 

SsinTa-jfz-Ko — (ct-cs))] 

(IT) 


, 2sin[a 2 (z-K 2 (ct-cs))] 
i|» 2 (z,ct) w W 2 exp[iic 2 z-i 1 *k 2 2 ct] ' 

ITIC 


z-k* 2 ZT (ct-cs) 
me 


We relate the constants a 2 , N 2 in ( 17 ) to the'eonstants a 1( Nj in (16) with 

the aid of the design conditions (3), (b). To this end, we examine the ratio 

* if 

sin[a 2 K 2 (ct-cs)] 




sintanKi — (et-cs)] 
* me 


( 18 ) 


The factor that follows the exponential should come out to be unity because 
of the design condition ( 3 ). This requirement yields 

a 2 K 2 B otiKi (19) 


$ 2 . 

Nl 


r 2, , or with (19), W 2 a 2 = N^aj . 


(20) 


The heights of the central peaks of the two wave packets are respectively 
| N iot 1 j and |N 2 a 2 |. They ore equal to each oth<°v, by virtue of ( 20 ) . 

Next, we determine the normalization factors Nj and N 2 . As we are 
dealing with a single electron, we have the condition 


/ 

z--« 

And, since <|i = < l'rH J 2 » ve obtain 


♦ 61 1 


(21) 


/ •fas W »+» TT »+oo w r +oo w 

dz\|»i 1 + j dz<|) 2 ^2 41 / ^2 14 / 8,zi|) 2 * 1 

z=_oo z--* Z=-« Z = -°° 


r+ m 


( 22 ) 
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These four integrals are evaluated most easily by way of the Fourier trans- 

* 

A At 

forms i|>i(K,ct) and $ 2 V K »ct). For, according to the Parseval relation, we 
have 

/ +# “ dzi^% 2 * o“ / + * * 2 3) 

Z*-» 


and similarly for the other three integrals. We can read out the Fourier 

transform ij!(K,ct) from the definition ( 9 ). But we must remember that we 

K, S'* , , \ 

linearized the dispersion relation, i.e., we replaced by 

Also we replaced ct by ct - cs, and we multiplied by the constant phase 

factor exp[-iK^Ki 2 ~ cs]. Thus, instead of ( 9 ), we now have 
* me 


Kl+Oj . . „u' 

i|>l(z,ct) * Nj / dK exp[iKz-i ) (ct-cs) - i & K i z — cs] 

We compare (24) with the inverse Fourier transform formula 


(24) 


♦ l(z,et) = 27 / “ dK * ill* ><&) 


(25) 


and read out 

iJJj^.ct) = 2irNiexp[-i ^(h K \ 2 + K i(*- K l ))(ct-cs ) -i fyej 2 cs ] 

.for I <_ ai 

$l(*c,ct) * 0 for |k-<i| > <*1 

A similar formula holds for $ 2 .( KiCt). We merely have to replace the subscript 
1 by 2. Because of the inequality (lO) and the related inequality o 2 « K2” k 1» 
the intervals in which $i(*,ct) and (|) 2 (ic.,ct) are non-zero do not overlap. Con- 
sequently, the integrand on the right-hand side of (23) is always zero. Thus 



» 

dzi|>j <i>2 = 0* and similarly 


/ fOD F 

dzd»2 

z=-°° 



(27) 


we obtain 
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Therefore, the functions ^(zjCt) and \|»2(z,ct) are orthogonal to each other. 
Since 


$ 1 *( * t ct)ij>i( K,ct ) » (2ir) 2 Ni N 1 for I’k-“<iI 1 ®1 

= 0 for | »c— ic 1 1 < ctj 

we obtain [from an equation similar to (23)] 

r+« 


/ dz i|/j * 2 ttNj Ni2ai. 
z*-“ 


Similarly , 


*>oo tt tt 

/ dz.^2 ^2 a 2 ttN 2 N 2 2ot2 • 


»V00 
Z=-“ 


The normalization condition (22) then gives us 

tt tt 

lnr(aiNi Ni + c*2 N 2 N 2 ) = 1 • 

When we combine this result with (20), we obtain 

‘ V*f«i+i 2 T ft’J 
N 2* N2 " Wa^Y 

Let us divide (30) by (29). With the aid of (32) and (19) we obtain 

rz s T , 

/:: ‘ az ' ^ ' vi 


(28) 


(29) 


(30) 


(31) 


(32) 


(33) 


(See (15) and the related equation for V2). This equation shows that the 

< _tc 

electron has a slightly higher probability (k, 2 > (q, ~ : x~ ' << to ^ eam 

#2 than in beam # 1, if the design conditions (3), (4) are fulfilled. 

tt 

Now we shall investigate the probability density xp t|». We calculate = 

(|)l + \p 2 from (16) and (17‘). In order to simplify the notation, we introduce 
the abbreviation 

j/ yf 

C = aiKi — (cs-ct) = a 2 <2 ~ (cs-ct) , (3*0 

* 1 1 me * me 
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where we have used <*i*x * a 2 »c 2 [see (19)]. We also use ajNj * a 2 N 2 [see (20)]. 
Then 


* ' a “lM Sl ;!niz 1Z ^ exptl *l z - i h*i 2 ct] 




(35) 


so that 


A . * Hgfff * 

♦ 2 siggzl cost (K 2 , Kl)z . f) g =t^ . 


The factors containing t; are envelope factors. The cosine term describes 
the pattern of beats. We see that - apart from the envelope factors - the 
beat factor is universal. This means that the peaks and valleys of the beat 
pattern of any particular electron coincides with those of its sister electrons. 
The SchrSdinger wave packets of different electrons are incoherent, because the 
normalization factors Nj for each electron contain uncontrollable phase factors. 

Nevertheless, as we have seen, the beat patterns are coherent , because they 

* 

depend only on the product Nj Nj in which the uncontrollable phase factors 
cancel. This coherence of the beat patterns is important in the application 
that we wish to examine later on. 

So far we have considered a single electron. Now we are going to evalu- 
ate the electron density for an electron beam, which consists of a great many 
electrons. We will ignore fluctuations (the shot effect). Thus we consider 
the electrons as evenly spaced, which means that the time interval ds between 
the instants when two consecutive electrons pass the origin z = 0 is given by 
ds = ^ , where R is the rate at which electrons pass the origin. Let us 
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express this statement in a different way, We are going to assign integer 
order numbers n to the electrons , in the order of their zero— passage time s> 
Thus, as the order number advances from n to n + dn, the zero-passage time 
advances from s to * + ds , and we have 

dn = Rds . (37) 

When 'we combine this with (34), we obtain (keeping z and ct fixed) 


, R 
dn = — 


dc . 


c -h 

me 

The (linear) electron density p(z,ct) (number of electrons per unit length 


along the beam) is given by 


p(z,ct) * /*" i|> 'I'dn , 


(39) 


n*-° 


( 1 * 0 ) 


or with ( 38 ) , 

p(z * ct) * t a^ ' , ’ ' ji ,C2 

me 

The integrand ^ 'I 1 is given by ( 36 ). 

We then have to evaluate integrals of the type 

<•+“ sin(c+c«i z) sin(c+« 9 z) / !.i\ 

1 12 “ 5 ’ 1 ' 

We use the Parsevnl theorem. It states that for two function f{ t) and g(c' 

we have 

r aewe) ]*«<?) - £ r «ktf(k)) # i(k) . iw 

£=-«> =-00 

where f (k) and g(k) are the Fourier transforms of f(c) and. g (0* In our 
case we have 


Since 


r(r) „ sin( ?;■*•« i.z_). , } „ 

sin(c+cnz) s JL [+ 1 dk e ika^ f 

2rr J 


5 +ajZ 


(1*3) 


( 1 * 1 *) 


k=-l 



mU6 
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we see that 


f(k) « ire ik0,lZ for |k| < 1 , 

* 0 for |k| > 1 . 

A 

Similarly for g(k). Thus 

^ rdk[?(k)]*5(k) . § r x , 

lc*“® k**l 


0*5) 


so that 

II . 2 " " (a 2 -ai)z 

The other two integrals we need are Ij j and I 2 ^ 2 . Their values follow from 


_ sin(a-j-ch )z 


( 1 * 6 ) 


(1*6) when we make a 2 oj and aj ■> a 2 . Thus 

Il t l * 12,2 * 71 • (**?) 

We insert (46) and (47) into (40), where the integrand was given hy (36). 

We obtain 


p(z,ct) * - 


R 

c 04*1 xfi 
me 


1 4(a 1 ) 2 N. 1 *N 1 2ir • 


ft 

Finally we use the expression (32) for Ni and also the relation (19). 
The end result is 

pU ’ ct) 1 t a fa;!ai ' )^ Z c °«CUa-«i)^ i 2 ) 

We see that there is a constant background, attributable to the 1 in the 
curled bracket, and a beat pattern, attributable to the cosine term. This 


(48) 


beat pattern fades with increasing distance from the origin, as described 

w -m,- 3in(a?-ctt )z . 

by the (a 2 -«i)z terJn ‘ 

We shall now discuss this fading term. According to equation (9), the 


quantity 04 is of the order of the quantum-mechanical wave number uncertainty 
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of the ^-packet. Thus is of the order of the momentum uncertainty. 
Similarly, ■•fia 2 * & 7 ^ c*i 3 'fiai is of the order of the quantum-mechanical 
momentum uncertainty of the if; 2 -packet. Nothing is known about these uncer- 
tainties. In the absence of any hard information about this matter, we use 
the simplest assumption, namely that we can make aj and a 2 i n (W) go to zero. 
Then the fading factor tends to unity, and (48) simplifies to 

p(z’,’ct) * 2. (~ •■ ^ --2-)“ 1 (l'»cos(( K 2 - K i)z- k(*2*- K \ 2 ) ct]} . (49) 

According to (49), the beats persist for all values of z. There is no fading. 

We feel fairly confident that we are permitted to replace (48) by (49)* 
because it is very likely that any fading of the quantum-mechanical nature 
is overridden by a fading of thermal origin, which we shall describe now. 

The electrons are emitted by a themionic cathode. So they come off with an 
energy uncertainty <5E, which is of the order of kT, where k is the Boltzmann 
constant and T is the cathode temperature. Since the cathode is hot, 6 E is 
of the order of ^ electronvolt . The same energy uncertainty 6 E obtains after 
acceleration and hence in the working region (near z= 0 ) of the apparatus. 

Now let us examine the i^-packet in the partial beam ft 1 . Since 

v 2 

- a"! • 


we have 


<5Et _ _ 5jej_ 

" a K, » 


El "1 

so that the thermal wave number uncertainty 6 *i is given by 

«*i - "i %• . (50) 

Now there is no uncertainty in the difference *c 2 2 - tcq 2 because of the design 
condition (4), namely 
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F B ~^^(^ 2 2 -Ki 2 ) , (4, repeated) 

where F is a precise constant. If we take the variation of (1*) and use 
$F * 0, we obtain 

*2^2 ” * 0 » 

so that 

$*•2 * ^ 1 . (51) 

*2 

Then we get [from ( 50 ) and (51)) 

% 

S(K2-»c 1 ) * - ^ (k2-*i) “ - ^(k2-ki) • (52) 

Therefore the quantity ( K 2 “ K 1 ) * n ^e cos i ne ”term of (49) has an uncertainty 
which is given by (52), whereas the quantity K 2 2 ~<\ Z is still exact. We take 
care of this uncertainty by taking a suitable average of the formula (1*9). 

The easiest average to take is the Gaussian average, as then the resulting 
integrals can be evaluated in closed form. We can afford to treat the quantity 
in the prefactor of (49) as exact, the resulting error is not important. 


Thus, if we denote averages by carets < >, we have 
<p(z,ct)> * ^(^-- ■ K ^ )~ 1 {1 + < cos[(k 2 ^<i)z - ^(k 2 2 -«1 2 ) “ ct] >}. 
Now 


(53) 


<cos[...]> = %<exp ![,,.]> + c.c, , (54) 

where c.c. stands for "complex conjugate". And, since we decided to take a 
Gaussian average. 


<exp i > = <exp[i( K 2 - K l)z - •§'( K * 2 2 ' K l 2 ) “T et]> = 

d * me 


= exp[- •^■( k . 2 2 - k i 2 ) — ■ ct] <exp(i( (C , 2 - K: i)z)> = 


/ + °° dn exp(i( K , 2 - K i+n)z ~(a 2 /2)n 2 ) 


* exp[- |( K - 2 2 -ici 2 ) “"ct] ^ 


/ dn exp(-( a 2 /2 )n 2 ) 
n--°° 


(55) 


MU9 
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where the quantity a is of the order of the reciprocal of 6(k2-kj) p as given 

by (52). Thus we nay write 

a , J! El {56) 

*2-*! «Ei * 7 

The integrals N and V in the numerator and denominator of (55) have the values 

J2* 


and 


H * — exp(i(k 2 -^l)z) exp(- jjgjO , 




Thus, (55) yields 

<exp[i(K2“ K l) z - 1 %( K .2 2 -^l 2 ) ^ ct]> * 

= exp(- ^y) exp[ i(* 2 -Ki)z - i ^(*2 2 ”*1 2 ) ^ ct l • 

Then (5*0 becomes 

<COs[(K.2- K .lh - H(k 2 Z -K\ 1 ) “ Ct]> = 

me 

= exp(- g^rjcosl (^2'-^ i )z - 2 Z ~* i Z ) ”” ct] • 

When we insert this into (53), we obtain the final answer 
<p(z,ct)> = f(^ ft )" 1 (ltexp(- ~j)cos[(^.2-^.i )z - H(k 2 2 -< \ 2 ) ~1), 

or, when written in a slightly different form, 

<p(z,ct)> * - - 1 {i+ exp ( - ~y)cos [ (k. 2 -**i ) ( z - K 1 ?' 2 ZT ct ) ] } . 


(57) 


2 me 

,2. 


(58) 


(59) 


( 60 ) 


c 'me 2 ^ 2a 

So now we again have a fading factor, namely exp(- g^y) . Let us calcu- 
late the fading length a. From (56) we have 


- li x __ 

a " 2 K.^.ic ^ 6Ex 


Ei 


( 61 ) 


Since the energy difference AE - E2 - Ej of the two partial beams is given by 


AE * ~ (< 2 2jc l 2 ) » 


we may write (6l) as 
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& , 

or, with Ei * tjj- «ci z , 


1 Ei 


I f% I n ± 

2 2m AE Wi * 
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a « k £l £l . (62) 

a V 2 AE <5E 

Since «2” k 1 k< '*i and E2-E1 « Ej # ve need no longer distinguish between kti 

and *2 or "between Ei and E2. Thus we can drop the subscripts 1 and 2 in (62) 

and write this formula simply as 

4 J5 J3 

a * < AE 6E * 


Finally we use E * k 2 or — * — 

dm k 


me 



a * 


~ , and obtain 

h /me 2 "' E E 
me sj 2E AE 6E ! 


(63) 


To repeat: AE is the energy difference between the two partial beams, it is 
typically of the order of 3 electronvolt . And 5E is the thermal energy un- 
certainty, which is of the order of y electronvolt. With E about 30 kilo- 


elect ronvolt , we have ^ 10^ and (1.5) * 10 s * Also me 2 a 500 kilo- 

* I * l X n 1 /inC ' /15o\ A AN. - 1, >. N .. N N ’ u . aQ ^ 



„ ~ 3. Thus a 'v 4 * 3 x 1.5 * 10 9 
u me 


y mC ^ 

— jgjjT 

V »/ a 

2 % 10 10 . Now — ■ is 3— timesthe Compton wavelength, or “= (3.86) * 


10” 11 cm. Thus finally, 


a * 1 cm . 


(64) 


So the best we can hope for is that the beats persist over a length of about 
1 cm. (Actually 2 cm, as there is 1 cm on either side of the origin. But 
we are only making rough estimates.) Of course we could stretch this length 
by installing a narrow-band energy filter in front of the beam splitter, thus 
reducing 6E. But then we would reduce the beam intensity. How the electron 
beam with beats will perform in its comtemplated use, which we shall discuss 
presently, will depend not only on the fading length but also on the beam 
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intensity. So, most probably, ve would be no better off in the end, had we 
installed the energy filter. 

We shall now describe the contemplated use of an electron beam with 
beats. The (linear) electron density in this beam is given by equation (60). 
Apart from a constant background (attributable to the 1 in the curled brakcet) 
there exists a sinusoidal density pattern, the beats, (attributable to the 
remaining term in the curled bracket). This beat pattern moves in the positive 
z-direction with a speed v given by 


x 

c 


2 me 


(65) 


The wave number of this beat pattern is *2“ Let us send a plane electro- 

magnetic wave in the direction of the negative z-axis, i.e, against the electron 
beam. We assume that this wave Is polarized with the electric field parallel 
to the x-aXis, i.e. normal to the z-axis. This electric field acts on the 
electrons in the beam and makes them oscillate in the x-direction. This 
oscillatory motion is, of course, superimposed on the steady motion, which is 
in the z-direction« The oscillating electrons give rise to a new wave, the 
reflected wave which propagates in the positive z-direction, i.e. against the 
incident electromagnetic wave. (The oscillating electrons produce also a wave 
in the negative z-directlon, i.e. in the direction of the incident wave. But 
this wave is of no interest to us.) Since the electrons in the beam exhibit 
a beat pattern, the process of reflection is similar to the reflection from a 
grating, as in x-ray crystallography. The reflected wave will be strongest 
when the Bragg condition is fulfilled, i.e. when the wave number of the 

Incident wave is related to the wave number <2“ K 1 °f the beats by 


K ln = l $( K 2“ K l)( 1 ” v / c ) * 


( 66 ) 
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The factor (l-v/c) comes in because the heat pattern is moving. We shall 
derive ( 66 ) after the discussion of the Doppler-shift formula { 67 ), Since 
the beats move with a speed given by ( 65 ), we are dealing with a moving 
grating. Because of this , the wave number k r of the reflected beam will be 
Doppler-shifted upward, according to 

K r * V iStI ‘ WV'lHW'M <«T> 

Now we are ready to derive equation (66). The Bragg condition is 

*r " ( -*in J ’ Vat * "a'" 1 * (67a) 

(Later on, when we develop the detailed theory of the reflection, we shall 

meet an independent proof of the Bragg condition for a moving grating. ) And 
indeed, with ( 66 ) and ( 67 ), the Bragg condition (67a) is satisfied. We thus 
have the possibility of constructing a highly selective frequency shifter. 

The selectivity comes from the fact that, for good reflection, the Bragg 
condition (66) must be satisfied. Expressed differently: One may . construct a 
selective moving mirror. The speed of the mirror is a good fraction of the 
speed of light. For electrons with an energy of 30 kilo electronvolt , equa- 
tion ( 65 ) gives a v/c-ratio of the order of 1/3. The frequency shift is 
appreciable. For v/c = 1/3, equation ( 67 ) shows that the wave number, and 
hence also the frequency, of the reflected wave is twice that of the incident 
wave. 

All of this sounds rather attractive. Nevertheless, the device is useless, 
because the reflection coefficient, i.e. the ratio of the reflected power 
density to the incident power density is exceedingly small. The reason is that 
there Just are not enough electrons available for a strong reflection to occur. 
Let us make a rough estimate to show that this is the case. A realistic 
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estimate for the electric current density J in the electron beam is J * 

1 milliamp per millimeter squared * 10” 3 amp * 10” 2 cm 2 * 10” 1 amp/cm 2 , 
about what may be found in a television tube. The electron current density 
is then — J, where e is the electronic charge * 1.6 x 10 19 coul. Thus ^ J 
is of the order of 10 18 electrons/sec cm 2 . The spatial electron density is 

then — J — . where v is the speed of the electrons. With v = 10 10 cm/sec 

e v 

(for v/c = 1/3), we get — J — - 10 8 electrons /cm 3 . But only a length of 
the order of the fading length a « 1 cm contributes to the reflection by the 
beam. The areal density a (number of electrons per unit cross sectional area 
of the incident electromagnetic wave) is than a = ^ J ^ - 10 s electrons /cm 2 . 

This is much less (by a factor of 10~ 8 to 10” 7 ) than the areal density of 
•optically active electrons in a monomolecular layer of a solid. Now we 
certainly do not get much reflection from a monomolecular layer, let alone 
a layer that is sparser by a factor of 10“ 8 to 10 -7 . 

The reader whom the preceding estimate convinces that the moving-mirror 
device is useless can stop reading right 'here. We continue a formal demon- 
stration to satisfy the more skeptical reader. Another reason is that the 
derivations that follow are instructive. They may be helpful for the feasi- 
bility analysis of other devices that one may wish to consider. 

We start with the equations (256) - (262) of the Appendix. These equations 
are the result of the perturbation calculus. We repeat these equations for 


the sake of easier reference. 


(i^ + 


( 68 ) 

(69) 





I 


ruo * nn A( !T y 


(Mt^l + »Xi) p = - (^o)p “ (-e)cA Q , 

+ H ^’XlL = " «X0)„ * 7 » 


P c 


TlhtT 2 ' cl o + v * C V ■ 0 * 


/ a ; - !ys CA X + V * (V X cAj) = -i- <P E (-e)J 2 '(xo%o + XO'J'O*),, » (73) 

' 3et; o T, low electrons p 

#P 

3 2 

( CA 2 + V x (V x cA 2 ) 

■ ^ I (-e)^(xo%i + xi%o + xo^i* + xi'l , o) T3 • (7*0 

o T,low electrons p 

#P 

Here in and (-e) are the electronic mass and charge. The index p refers to 
electron #p. The quantities with the subscript 0 are the zero-order approxima- 
tions. Those with subscripts 1 and 2 are the first-order and second-order 
corrections. 

We are not interested in equation (73), which tells us how to calculate 

the first-order correction of the electromagnetic field. To be sure, this 

correction does not vanish in general. However, as we showed in Section 3, 

there is no radiation that is attributable to these fields. It is for this 

reason that we are not interested in cAj. Furthermore, in our example in 

which the zero-order SchrSdinger waves of the electrons are plane waves in 

z-direction, the sum in (73) is purely longitudinal in character. Thus, when 

the projection operator J* acts on this sum, the result is zero, since OP 

T,low _ T,low 

retains only the transverse part, which is zero. Therefore, the right-hand 

side of (73) is zero in our example, so that (73) is satisfied by cAi = 0. 

As a consequence, the first-order correction cB* = V F cAi to the magnetic 


ORIGINAL PAGE IS 
OF POOR QUALITY 


M55 


field vanishes, and the first-order correction Ej B - simply 

Ej ■ - 741, where 4>, in "the electrostatic potential associated with the zero- 

order approximation to the charge density. A straight forward calculation of 

electrostatics, which we will not reproduce here, shows that the alternating 

part of Ei, which is attributable to the beats, is only of the order of 10“ 4 

volt /cm, much too small to be of any interest. (This magnitude of an electric 

field occurs in a laser beam with a power density of less than 10” 10 watt/cm 2 . ) 

Having disposed of the first-order correction cXi as something that is of 

no interest, we concentrate our attention on the second-order correction cA 2 , 

which we can calculate from equation (7*0. Even though c£ 2 will turn out to 

be a genuine radiation field, it will be too weak to be of any use. Let us 

substantiate this prediction by means of an explicit calculation. First we 

have to evaluate the quantities (tjh) and (xi) on the right-hand side of (74). 

P P 

The remaining quantities (^o) and (xo) are given by the properties of the 

P P 

electron beam. First we note that equations (70 ) and ( 71 ) are satisfied by 

('J'l)p = 0 » (T5) 

™(Xi) p = -,Uo) p ~ (-e)cA Q . . (76) 

— iK — 

The reason is first that xo Vtpn (from (68)) has only a non-vanishing 

ra 

z-component, since t|'o <ioes not depend on x and y. Furthermore , by assumption, 

i 

cAo has only a non- vanishing x-component, so that the dot product on the right- 
hand side of (71) vanishes. In the second place, on taking the divergence of 
(70), we obtain 

^(V.V^i) + m(V *xi ) = - ^(-e){( 7(jj 0 ) *cA 0 + (<|»o) V*cAo) . 

But the right-hand side of this equation vanishes, since ^•c'Kq * 0 (because of 
the Coulomb gauge) and the dot product vanishes (VifiQ has only a non- vanishing 
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z-component and cXg has only a non-vanishing x-component ) . This equation 
then becomes 

) + m(7*xi) = 0 . (TT) 

r Jr 

When ve combine (77) with (71) and use the fact that the right-hand side of 
(7l) vanishes, we obtain 

+ —■ V.fyj) = 0 . (78) 


+ S 5,5 '-p 


This equation is satisfied by (’i'i)p = 0, so that (75) holds. Equation (76) 
then follows from (70) and (75). 


Equation (7M then simplifies to 
— j cA 2 + V K (7 x cA 2 ) = - " -- 


( 3ct ) 


some' 


to 

l (^0 ^oL 

electrons F 

U 


cAq 


(T9) 


But the expression in the square bracket is the spatial electron density. 

It is related to the linear electron density <p(z,ct)> of equation (60). 

We merely have to divide by the cross section Q of the electron beam. Thus, 
with the abbreviation (65), 

l (^0*^0 )_ " exp(- 'r| 2 ')cos[(V 2 -K' 1 )(z - J ct)!} . (80) 

electrons P * 

fa 

TD 

We note that ^ is the electron current density of the beam. The electrical 

R 

current density is then (-e) ^ . 

Now we combine equation (79) and (80). In order to simplify the writing. 


we introduce the classical electron radius 


(S ’ 82 > • 10 ’ 13 > 


and the symbol 


j -5. 

J e Q 


(81) 


( 82 ) 


for the. electron current density of the beam. Thus 


We can delete the projection operator because the right-hand side 

of (79) is purely transverse and has only low-modes components. 


M57 


^SSy; 


(3ct)2 


cA 2 + $ * * CA 2 ) 


- ~ ( 7)” 1 r o J e (l+exp(- ■^•)cos[(h: 2 -k.i)(z- ^ ct)]}cAo 


(83) 


Finally, we use the fact that, by assumption, the incident electromagnetic 
field is polarized in the x-direction, Then the desired field c5 2 will also be 
polarized in the x-direction. Then, with the Cartesian unit vectors I, J, R, 
we can write 

cAo(z,ct) * icAo(z,t), 
c5 2 (z,ct) * IcA 2 (z,et) v 

where we are stipulating that both fields do not depend on x and y. We also 
have 


(8U> 


V * cX 2 * S * I ~ cA 2 * 3 cA$ , 


and 


■? r 


r 3‘ 


V x (7 x cA 2 ) = k * ] 7 cA 2 * - i jjpr cA 2 

Then, when we delete the unit vector I, which is common to all terms, equation 
(83) becomes 


2 g2 

— 7 cA 2 - — cA 2 


T8ctF VAZ 3z2 

- - — (— ) r J {l+exp(- •r| 7 )cos[(K 2 -ic 1 )(z- ^ ct)]}cA 0 . ( 85 ) 

c c o e d &*■ c 

We then have to find the solution *J>(z ,ct ) of a partial differential 
equation of the type 

f(z>ct) • (36> 
Here both z and ct run from - 00 to + 00 , and f(z,ct) vanishes in the distant 

past, i.e. for ct -*■ The explicit solut?lon of this differential equation is 

. z+( ct— c k) 

i|i(z,ct) = \ j c dcir / ^ f(?,cr) 

CT — — 00 C=Z-(ct-Ct) 


( 87 ) 
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There is no connection between the dummy variable C of this equation and 
the C of equation ( 3 ^). 

This writer could find the explicit formula (87) in only one book, 

namely " Mathematical Physics" by Eugene Butkov, Addison-Wesley , 1968, pg. 607. 

As this book is not readily available, we establish the validity of ( 87 ) 

here. We will not go through the process of creating this equation. This 

process is based on physical intuition, which comes from the examination of 

vibrating strings? officially it goes by the name of Green function techniques. 

Instead, we will check that the differential equation ( 86 ) is satisfied by 

the solution ( 87 ). We calculate the required partial derivatives. We start 

with 1 p. In ( 87 ), the ct occurs in 3 places. On taking the derivative 

with respect to each place in turn we obtain 

z+(ct-ct) 

/ at f(c,cx) + 

C=z-(ct-cr) at ct=ct 

+ % / Ct dcT - f(z+(ct-ct) ,ct) + 
cn=-°° 

+ % / ct dcTf (z-(ct-c t) , et) • 
cr=-°° 



The first term vanishes, because the domain of integration has shrunk to 
zero. In each of the other terms, the ct occurs in two places. If we denote 
the first partial derivative of f(z,ct) with respect to the first argument 
by g(z,ct),we obtain 

3 2 _x, 

; 3 Ct T2 ♦ = ^(f(z,ct) + / dex g(z+(ct-c.r}cT) + 

v c T=— 00 

+ f(z,ct) - / ct dex g(zi-(ct-cx) ,c t) J , 

Ct=-® 


or 
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7T-TV7 'll * fU.ct) + \ / ct dcx [g(z+(ct-cx) ,cx)-g( z-(ct-cx) ,cx)] . 

' a ' C X *- 00 

Now we come to the z-derivatives . In ( 87 ), the z occurs in 2 places. 
Proceeding, as we did before, we obtain 


( 88 ) 


U * h f Ct dcx [f(z+(ct-cx)-f(z-( ct-cx), cx)] , 
Cx“-« 


and 


I = h / ct dcx [g(z+(ct-cx) - g(z-(ct-cx),cx)3 . 


(89) 


CT -— 00 


When we subtract ( 80 ) from (88), we see that (-86) is satisfied indeed. 

We should require not only that the differential equation (86) is 
satisfied by the solution ( 87 ), but also that the solution contains only 
outgoing waves. Let us check that the second requirement is fulfilled as 
well. We assume that the excitation f(z,ct) is confined to a finite domain > - 
of the z-axis, say to the domain |z| < b, where b is some constant length. 

We therefore assume that 

f(z,ct) * 0 for jz| > b . 

Now let us examine the time derivative , rather than t|> itself. We already 
calculated this derivative and found 

\ J ct dcx f(z+(ct-cx),cx) + 
ex*- 00 

+ % / Ct dcx f (z-(ct-cx) ,cx) . 
cx=-°° 

Let us pick a point z to the right of the excitation domain; thus z > b. Since 
ct-cx > 0, we have z + (ct-cx) > b, so that the integrand in the first term 
vanishes. Therefore only the second term survives. But in this term, z and ct 
occur only in the combination z-ct, which signifies waves that travel in the 
positive z-direction, i.e. outgoing waves. Similar statements apply to a 
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point z which is to the left of the excitation domain. Thus, (87) contains 
only outgoing waves. 

Let us now apply the solution ( 87 ) to the problem of equation ( 85 ). 

We obtain 


cA 2 ( 2 ,ct ) * - H ~ (~) r Q J e . 

, z+(ct-cr) _2 v 

•} dcr / dc(l+exp(- “s-)cos[tc 2 “ ,c l)(C- 7 cr) ]JcA At.cx) . ( 90) 

cr=-» C“z-(ct-ct) 


We observe that this result is a linear functional of the incident electro- 
magnetic vector potential cAq. Let us assume that cAo is a monochromatic 
wave, whose wave number is matched to the wave number k 2 ~* 1 of the beats, so 
that we can obtain the strongest possible response. Accordingly, we write 
(with (66)) 

cAqU,ct) = cAo 'cos [ (l- J) (c+cr)3 , (91 ) 

where cAg is a constant amplitude factor. The plus sign in the (c+cr)-terro 
indicates that the incident wave travels in the negative z-direction. 

A precise interpretation of the formula ( 87 ) would require that we make 
the right-hand side of (85) tend to zero in the distant past. We can achieve 
this by first decomposing (9l) into two exponential terms, i.e. by writing 
cAgU,c ) = '-5c A 0 ' exp [ ( 1- ^)(i -M ( ?+cT ) ] 

+ %ca 0 '«*pKi- 1 L )(;+= T )] . (92) 

K -4C K i- jc ^ ^IC 

And then we should replace the wave numbers i 1 ^ and -i by i + 8 

K “1C 

and i + 8 respectively, where 8 is a positive real quantity. We then 

would replace ( 92 ) by 

cAq(c,ct ) = ^cAo'exp[(l- ~)(i + 8 ) ( C+cx ) ] 

+ ’$cAo’exp[(l- J)(-i - 2 J? + S)(c+c T )] . 


(93) 
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The presence of the quantity B in (93) would ensure that cAg(c,ccr) tends to zero 
for ex with sufficient rapidity. We would then calculate the respence 

CA 2 separately for each exponential term in (93) and add the results, which we 
are permitted to do because of the linearity of the functional (90). In the 
end w6 would make 6 approach zero. But this concession to logical precision 
is not worth the effort. Therefore, we shall evaluate (90) on the basis of 
expression (91). 

Furthermore, we shall ignore the constant background of the electron beam, 
when we calculate the response cA 2 » as being of no great interest. Accordingly 
we delete the terra "l” in the curled bracket of (90). Thus, with (91), we 
obtain 


cA 2 (z,ct) = - % ~ (J) r J cA 0 * . 
beats c c 0 e 


. z+(ct-cr) 2 

f dq exp(- •=^2 -)cos[(k2-<i)(?- — ct)] 


• r *ca / 

CT=-® C=Z-(ct-CT) 


2a 2 

cob[( 1 - J) ( 5 +ct)l . 

c d 


For the sake of greater convenience, we calculate, not cA 2 itself, but its 

g 

time derivative repeat the differentiation process that we used 

when we checked the validity of (87) and obtain 

3cA2 (z.ct) - - H (J) 1 Ve cA °' • 


Set 
beats 

/ ct dc T {exp(- ^^|^7^-^)cos[(K 2 -Ki)(z+ct-(l+ ^)cx] 
ct=-» 

cos[(l-J) Zz—L (z+ct)] ♦ 

C d 


+ exp(- ^ • Z "~ 2 ^ F' T , ^ ~) cos ( ( k 1 ) ( Z“Ct+ ( 1- 7 )cT)] 


cos[(l- — ) (z-Ct+CT)]) . 

C d 


( 95 ) 
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We assume that the field point z is far beyond the region of beats in the 
direction of the positive z-axis. Then z > 0, and z » a. Since also ct-cx 
is positive, ye see that the exponential fading factor in the first term of 
(95) is essentially zero. We can therefore delete the first term. In the 

second term we rewrite the product of cosines according to the formula 

> 

COS^COSH * JgCOs(5+Tl)+%COs(5-n) , 


Then 


1 frt (z > ct) * - * T ( ?" 1 r -. J » cA °' 


o e 


beats 

/ ot dot exp( illzigpsX. ) . 


CT =— 00 


• (cob[ — (2+(l- “)(z-ct) + 2 ( < 2~ k 1 ) ( 1" “) ct 3 + 

+ cos [" 1 ?■"* 1 ( 1+ 7) ( z-et ) ] ) . ( 96 ) 

d c 

Now we perform the integration over ex. Let us look at the first cosine 
term in (9^). This is a rapidly varying function of cx with average value 0. 
The rapid variation is attributable to the term 2 (k z -k i)(l- ^)cx in the argu- 
ment of the cosine. On the other hand the exponential fading factor is a com- 
paratively smooth function of cx. It is comparatively smooth, because 
(ic 2 _*c 1 ) (l - ”) is much larger than the inverse fading length — . As a result, 
the contribution of the first cosine term to the integral is essentially zero. 
There remains the second cosine term, which does not depend on the dummy 
variable ct and can therefore be factored out. The only thing that needs to 
be done is the intergral over the exponential fading factor. With the substi- 
tution £ = cx + cz-ct, we obtain 

ret , / (z-(ct-cx:) ) 2 \ _ tz „ t g 2 \ 

J def exp( ^7 ) = J <U expl- • 

or =-°° £*-• 
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Since z >> a, we may replace the upper limit by +•, without making an 
appreciable error. Thus 


/ ct dc t exp(-^" Z **^°|~° ^ ^ ) ■ a . 
c *-*» 

On combining (96) and (97) we. obtain 


3cA^ 

3ct 
beats 


z,ct) * -a /2 iT — (— ) r J cAo* cos[k (z-et)] 

c c o e r 


(97) 


(98) 


where we have used the expression (67) for the wave number tc of the reflected 

r 

wave. This wave propagates in the positive z-direction, as indicated by the 
minus sign in the cosine term, i.e. in the direction opposite to that of 
the incident wave. 

The transverse electric field E 2 X of the reflected wave is given by 

E 2 * - 3^ oAj. (99) 

(There is no contribution - from the scalar potential ♦2- beats • 

the first place, the domain of non-zero values of $2-b ea t s is conf i ned to the 
region of the beats, since < t*2- bea ^. s is a solution of the electrostatic Poisson 
equation 


^ * ^ 2 beats eoc cp) beats * 

which does not yield any propagating waves. On the other hand, the region 
for which we calculated the reflected wave is beyond the region of the beats. 
In the second place, where the scalar potent! a], does not vanish, its gradient 
is in the z-direction, thus purely longitudinal). We compare the x-components 
of the reflected and the incident waves. For the incident wave, we have, 
from (9l) and (66), 

Eg(z,ct) = - cAg * K in cAq' sin [^.(z+ct)] , (100) 
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whereas, for the reflected wave, we have, from (98) and (99), 
E 2 b eats^ z,ct ) * c (~)~ l i* 0 J e cAo * cos [k ^ ( z<-ct ) ] 


( 101 ) 


The ratio 11 of the reflected to incident power density is the square of the 


ratio of the amplitude factors in (101) and (100), so that 

R * ( ? _1 Ve > 2 • 


( 102 ) 


If we introduce the wavelength X. * — — ■ of the incident radiation, equation 

1 K in 

(102) becomes 


R ■ (a X lnJ?7 r oV 2 • 


(103) 


Let us insert numbers. For the fading length a we take the optimistic value 
of 1 cm. Suppose we adjust the device to an optical wavelength A in * 5 x 10" 5 cm. 
We also have £ ir . Then for J *v 10 18 cm" 2 sec, corresponding to an elec- 
trical current density of the order of 10" 1 amp/cm 2 , we obtain, with (8l), 


R a. 10" 18 . 


(10U) 


This exceedingly small value of the reflection coefficient R shows that the 
device is useless. This verdict stands, even if we boost the. electron beaju 
current by many orders of magnitude. Our detailed calculation substantiates 
our previous estimate, which was based on the sparseness of scattering electrons. 
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Appendix: The Semiclassical Method. 

The problem that we wish to address in this appendix is that of the inter- 
action between an atom (or molecule) and electromagnetic radiation. We confine 
ourselves to the semiclassical method in which the atom is treated quantum-mechani- 
cally, i.e. by way of the SchrBdinger equation, whereas the radiation is treated 
classically, i.e. by way of the Maxwell equations. Of course, quantum mechanics 
could be used also for the electromagnetic field. But for the kind of applica- 
tions we have in mind this is not necessary. 

Most textbooks of quantum mechanics discuss how a given electromagnetic 
field influences an atom. But what is hardly ever explained is how, in return, 
the atom reacts back on the field. This aspect is, of course, important in 
laser physics. Only a rule is sometimes used, but hardly ever explicitly stated 
nor derived. According to this rule, the atomic sources of the electron-magnetic 
field are the quantum-mechanical averages of the electric charge densities and 
the electric current densities. Here we wish to derive this rule. 

The derivation is based on the following consideration. We are dealing 
with a two-way process: Forward from field to atom and, in reverse, from atom 

to field. We know the equations of motion for the forward process, and we wish 
to infer the equations of motion for the reverse process. In order to make this 
inference one employs the action principle of analytical dynamics. One finds, 
as we shall see, that with the application of this principle the equations of 
motic'i for the forward process will determine those of the reverse process. 

In order to instill confidence in this procedure we shall employ it first in the 
case where both the atom and the radiation are treated classically. In this ca-e 
the outcome is, of course, known: The Maxwell equations will be the result. This 



case, then, provides a good confirmation of the procedure. Later on we shall 
apply the procedure to the case where the atom is treated quantum-mechanic ally. 

In order to state and use the action principle one has to define a number of 
terms. We will give these definitions, not in abstract generality, but for the 
concrete example at hand, an atom interacting with electromagnetic radiation. The 
terms that we wish to define are, in this order, the instantaneous configuration, 
the history, the instantaneous' dynamical state, and the driver fields. 

1. The Instantaneous Configuration. 

As far as the atom is concerned, its configuration is specified by the 
locations of each of the atomic” particles, if we neglect the spins of the par- 
ticles. These locations may be specified by a set of n cartesian coordinates, 
three for each particle. Thus n, the number of degrees of freedom of the atom, 
is three times the number of particles. In order to avoid cumbersome language 
and notation, we will say that the configuration is specified by a point q in the 
n-dimensional "configuration space" of the atom. If we so desire, we may specify 
the point q by the set {q 1 , q 2 , ... q 11 }, abbreviated by {q 1 }, of n "configura- 
tional coordinates" chosen In any convenient way. ({...} means "set of". Latin 
indices run from 1 to n. ) 

Each point of the configuration space determines an electric charge pattern 
in the three-dimensional space in which all of us live, our "home-space". 

Because of the print-like nature of the atomic particles, this pattern is, 
technically speaking, a distribution. The electric charge is concentrated in a 
number of discrete points (one for each particle); everywhere else the charge 
density is zero. So as to avoid the technicalities of the distribution theory 
we will envision the atomic particles as being of finite, though small, extent. 
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We also imagine that the charge density of each particle gradually decreases 
from a high value at the center to zero at the surface of the particle. We 
can then regard the charge density p of the atom as s vhole as being a smooth 
function of the home-space position vector r and of the n configurational coor- 
dinates q*. When we wish to indicate what variables the charge density p depends 
on we denote it by p(q*,r) or sometimes by p(q,r). - As the text of this para- 

graph shows, two spaces are involved in the theory, the home-sapce and the atomic 
configuration space. Vectors in the home-space will be denoted by an upper bar, 
as we have already done it for the position vector r. 

We have completed the discussion of the atomic configuration. We can now 
go ahead and define what is meant by the instantaneous configuration of the elec- 
tromagnetic field. Here we shall use the familiar notation and terminology of 
engineering electromagnetics, although there exists a more felicitous notation, 
namely that of the exterior calculus, which conforms more closely to the geometric 
imagination. As most readers may not be familiar with the latter notation, we use 
the former, thus trading ease of visualization for familiarity. All quantities 
and formulas will be expressed in the SI system of units (Systfeme International). 

One's first inclination might be to describe the instantaneous configuration 
of the electromagnetic field by means of the following two fields: the electric 
field strength E and the magnetic induction B. But then one does not account 
for the fact that E and B are subject to the constraint expressed by the first 
pair of the Maxwell equations, namely 

7x1+ ~r cB = 0 and 7 • cB = 0 „ 


(1) 


Ah 
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(We have inserted the factor c (c * speed of light in vacuo) in various places. 

The advantage is that E and cB are expressed in the same unit, namely volt ra” 1 
and that the differential operators V and glaring in the same unit, namely m” 1 ). 
In order to enforce this constraint, one adopts the usual procedure of expressing 
and cB in terms of the vector potential A and the scalar potential ♦, as shown 

E = - (-g—r cA + V<t>) and ci = V x cA (2) 

(Both cA and <f> come in the same unit, namely volt.) The fields of c A and $ are 
free of constraints. We use* them to specify the instan- 

taneous configuration of the electromagnetic field. 

At this point we must say something about the boundary conditions for the 
electromagnetic field. We imagine that the atom is located inside a large room 
(singly connected and bounded by a single surface, like a room in a house) with 
perfectly reflecting walls. The boundary conditions for E and cB are then that 
E is normal to the wall and that cB is tangential to the wall, i.e. that 

Af x E = 0, and Af • cB = 0 , (3) 

where Af is the vectorial surface element, directed outward. We can enforce (3) 
by stipulating the following boundary conditions for cA and <f>. 

Af x cA = 0, and $ = 0 at the wall. (4) 

The instantaneous configuration of the entire system comprising the atom and 
the electromagnetic field is given by the point q, (or the n configurational 
coordinates q*) of the atomic configuration space and the two field functions <J>(r) 

and cA(r) in the three-dimensional home-space. 

2. The History . 

The history of a system is determined, when we specify the instantaneous 
configuration as a function of the time t of its occurrence. In our case then 


mi 
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we specify the history by the functions q*{t), $(r,t), and cA(r,t). In general 
the history may be made to conform to any prearranged plan. But then the system 
will have to be driven by means of suitably chosen driving agencies, which will be 
discussed in section U. 

3. The Instantaneous Dynamical State . 

The instantaneous dynamical state of a system is determined, when we specify 
the time t of its occurrence, the configuration {q*,<J>(r) , cA(r)}at that time, and 
also the time derivatives of the configurational parameters at that time, i.e. 

i 3 g — i 

the set {v , rr- <$, — cA}, where we have used the customary abbreviation v 
d i 3 . 

for -rr q. . The partial derivatives tt indicate that the time derivatives are to 
ox 3« 

be taken while the position r in the home-space is held fixed. Altogether then, 

f i i • 3 3 

the instantaneous dynamical state is given by the set *4t, q , v , $, cA, c 

As far as the fields $, cA, etc. in this set are concerned, they have to be speci- 
fied throughout the room in which the fields exist. The utility of the concept 
” instantaneous dynamical state” arises from the fact that the specification of 
this state at one single time, e.g. at t = 0, is enough to determine the entire ' 
history, provided that the system is ” free-running" , i.e. not driven. 

We have seen that the instantaneous configuration of the atom, i.e. the 
set {q*}, determines a charge density field p(q^,r). Similarly the instantaneous 
dynamical state of the atom, i.e. the set {t, q*, v^} , determines not only a 
charge density field pCq 1 ,?) at the time t, but also a current density field 
^(q^jV 1 ,?) at that time. The current density field will be linear in the v* 
so that we can write 

JfaSv 1 ,?) ■ v^J^ (q 1 ,?) , 


(5) 
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where we have used the summation convention: We will always sura over any repeated 
index, here j, which appears both as a superscript and a subscript. Latin indices, 
range from 1 to n, where n is the number of degrees of freedom of the atom. 

One should keep in mind that the functions pCq 1 ,?) and Jj(q*,r) clo not expli- 
citly depend on the time t. The time-dependence of the charge density field and 
of the current density field is implicit* it arises from the fact that, when the 
history is examined, the q* and v * are seen to be functions of the time. 

The continuity equation 

P + V • J = 0 (6) 

entails a relation between the functions pfq*,?) and J (q 1 ,?). Since p depends 

J 

cn t by way of the q*, we have 


af- p * (r^Ap = 


‘ dt 


Also, from (5), 


J 


V 


= V J V . 


V 


V 


since the v° do not depend on r. Thus the continuity equation yields 

v* [~L. P + v . j ] = o . 

3q J J 

This equation must hold for any arbitrary choice of the velocity components v^ . 
Thus we get 


8q- 


j p + V • = ° . 


(T) 


3 i 

One should note that the derivative — r attacks one of the q in the functional 

i - K - - i - 

form p(q ,r) whereas the divergence V-attacks the r in J^(q ,r). 


Equation (7) is not the only relation connecting the functions (q 1 ,?) and 

J .(qSr). There is another important relation, which we derive in Note # 1 . 

J 


It is 


I 


of'2'M 




( 8 ) 

Sq 1 J 3q J 1 0 1 J 

The origin of equation (8) is more subtle than that of the continuity equation (6), 
on which equation (7) was based. Equation (8) is derived from the principle of 
’'gene-identity”. This word, meaning ident if lability, was coined in Germany, per- 
haps around the turn of the century, probably because it is easier to pronounce 
than its synonym. We say that the atomic system obeys the principle of gene- 
identity, if each element of charge can be regarded as an identifiable object, 
i.e. if it can be tagged. This is certainly the case if each atomic particle 
is a point-like object. We shall make the assumption that the n + 1 functions 
Jj (<!*,?) are structured in such a way that gene-identity is fulfilled 
even for smeared-out, i.e. continuous, charge and current densities. Equation (8) 
is the necessary and sufficient condition for gene-identity. 

Equation (8) is an important link between the history of a system and the 
driver fields. These fields will be discussed in the next section. - 
U, The Driver Fields . 

, In Section 2 we introduced the notion of driving agencies. These agencies 
are fields in the home-space, namely a force field with force density F d (r,t), 
an electric charge field with charge density p d (r,t) and an electric current 
field with current density J d (r,t). These three driver fields have to be 
impressed from the outside, i.e. from sources that are external to the system. 
Figuratively speaking, the driver fields are impressed "by hand” . 

We have affixed the superscript d (”d” for driver), mainly in order to 
distinguish the driver fields p and J from the atomic fields p and J (without 
a superscript), which were discussed in the previous sections. The driver fields 
p^ and have to be superimposed on the atomic fields p and J, so that the total 
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fields and ("t" for total) are given by * p + p d and = J + J d . In 

d. "”3 

actual practice the driver fields P and J Q are spatially veil separated from 
the atomic fields p and J, so that the superposition, when viewed as a technical 
task, is easily accomplished. 

We have to make a comment about the forces whose density we denoted by F d . 
Here only those forces are included which are impressed directly on the atomic 
constituents by the "man” who steers the system to make its history conform to 
a prearranged plan* Of course, this "man'would also have to impress forces on 
the electrical drivers, whose densities we denoted by p d and J d , in order to over- 
come inertia effects and to counteract the electromagnetic forces on the drivers. 
(The density of these electromagnetic forces is given by the usual expression 
p 1 ^ + J d x b). Perhaps the word "interface” used by computer engineers may help 
to explain why the density should pertain only to those forces that are im- 
pressed directly on the atomic constituents. The interface between the system 
under study (namely the atom and the radiation field) and the external world 
(namely the "man 1 ’ who steers the system and his actions) consists of nothing but 
the three driver fields r , p , and J . Thus the forces that the "man" impresses 
on the two electrical drivers (whose densities are p d and J d ) are entirely within 

the external world; they do not penetrate the interface. 

# 

When the history of the system is made to conform to some prearranged plan, 
then the system has to be "steered" by means of the three driver fields F^, p d , 
and J d . The required driver fields can be calculated by way of the action 
principle, which is the subject of the next section. The case that is usually 
treated in the literature is that of a "free-running" system, i.e. a system 
whose history unfolds according to spontaneous evolution. A free-running system 
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then need not be driven, so that the condition that a system be free-running 
is that all three driver fields be zero. This condition yields the veil-known 
Lagrangian equations of motion. 

5. The Action Principle . 

With each history (q*(t), $(r,t), cA(r,t)> and a set of two times t = a, 
t * b > a there is associated a quantity A . , the action accumulated between 
these two times. The units in which A is expressed are those of an action, 

& 9 D 

for instance Joule sec or Joule hertz” 1 . It is given by 


*a,b * 4 a dt L ’ 


(9) 


where L, the Lagrangian, is a function of the instantaneous dynamical state. 
The units in which L is expressed are those of an energy, for instance Joule. 
The nature of the function L is obtained essentially through revelation, of 
course not divine revelation, but revelation stemming from the long tradition 
of science. But one ought to check the revealed result against other things 
we know. And that we shall do to some extent later on. For the system con- 
sisting of an atom and an electromagnetic field, the Lagrangian L is the sum 
of three terms, the atomic Lagrangian L . the electromagnetic Lagrangian 

6L uOBl ) 


L , and the interaction Lagrangian L. . . 
q .m, int 

L - L . + L + L. . . 

atom e.m. int 

We now proceed to write down these three terms. In the non-relativistic 


( 10 ) 


approximation, to which we shall adhere, the atomic Lagrangian is given by 


L = — v^M* v'^ + u'v^ - V* ; v* = 
eitom 2 M iJ V U i * T dt * 


( 11 ) 


where the quantities M^, U’, V 1 depend on the atomic configuration {q }, but 
not explicitly on the time t. (in general these quantities could be permitted 
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to depend explicitly on t, but such a t-dependence is not needed for our 
system.) We have affixed primes to M^ t V, because ve wish to reserve 
the unprimed symbols for related quantities that are of more immediate practical 
importance. They will appear later on in our development. The first term in 
(il) is the kinetic energy, the last term is the potential energy. And the U’ 
are the components of a vector potential in the configuration space. For par- 
ticles without 3pin all the U' are zero. We shall restrict our treatment to 
spin-less particles. 

The electromagnetic Lagrangian L is the volume integral of the 

e.m. 


Lagrange density L of the electromagnetic field. 

e.m. 

L = /// AtL 
e.m. e.m. 


( 12 ) 


where At is the volume element. The integral is taken over the whole room in 


which the atom is situated. The Lagrange density L is given by 

e.m. 

*•.». = ! |E.E- %cB.cB , 


(13) 


where £ is the permittivity of vacuum. The two terms in (13) are, respectively, 

the electric energy density, and the magnetic energy density. Equations (2) 

enable us to express L more directly in terms of the configurational para- 

e.m. 


meters cA and <p of the electromagnetic field: 

3 .7 


+ + W --f <’ * eA) 

The interaction Lagrangian L^^ is a similar volume integral, this time 

of the interaction Lagrange density L. . . 

mt 

L int = M ArL int * 

The integrand is given by 

L 


(Ik) 


(15) 


’int 


-• (cA»J - <f>cp). . 
c 


(16) 


All 
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Note that 5 and cp are expressed in the same units, for instance amp nT 2 . Accord- 
ing to (5), equation (16) may be written in the more directly useful form 

L int * c ( c ^ # ^i vl - ^P) • (IT) 

Perhaps it is advisable to repeat (17), this time with all the independent 

variables indicated: 

I' int (<l i »v i # r,t) = — [c5(r ,t ) i Jj.Uf.rb- 1 - $(r,t )cp(q^ ,r) ] . (l8) 

The action principle is concerned with the variation 6A rt . of the action, 
when the actual history (q i (t), (r ,t ) , cA(r,t)].is replaced by a slightly varied 
history {q*(t) + 6q*(t ) , <|>(r,t) + 6<f>(r,t), cA(r,t) + 6c5(r;t)}. Here the varia- 
tions 6q i , 64», 6eA of the configurational parameters are very small in some sense 
(i.e. almost infinitesimals), but otherwise arbitrary* One then evaluates the 
variation of the action 6A ft ^ = A & ^ (varied history) - A ft ^ (actual history) to 
first order in the 6q*, 6$, 6cJ with the aid of the rules of the calculus of 

variations'. Although this is a very common calculation, we will repeat it here, 

« 

because it uses and generates some interesting ideas. 

But before we start with this calculation, we should adopt a policy of pru- 
dence and eliminate the somewhat vague concepts of "very small in some sense" and 
"almost infinitesimal". The procedure that we shall follow is also common. In- 
stead of examining just two histories, namely the actual one and the varied one, 
we examine a whole family of histories, ordered and labeled by some "variational 
parameter" n. The case of n = 0 corresponds to the actual history; and the more D 
deviates from zero, the more does the history labeled by n deviate from the history 
labeled by 0. Accordingly, the configurational parameters q*, 4>, and cA will 
depend on the label ru Thus the history labeled by n is given by the set 
(q i (n,t), <t>(n,r,t), cAvn,r,t)l Since q* depends not only on the 

# — , 

When we vary the action A we keep the end points a and b of the surveillance- 
time interval fixed* However , in some investigations, one finds it useful to 
vary also the end points by the amounts 6a and 6b. But we can get along without 
the introduction of these more general variations. 
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tine t, but on the variational parameter n as well, the configurational velocity 
component v 1 is given by the partial derivative ^ . Then we simply def in e 


fiq 1 , fi$, ficA, and fiA b by 

fiq 1 5 -4- «n, «♦ = ficX 5 ^ fin, «A a>b = 


3A 


a.b 

3n 


fin 


. i . 3V 1 . 
fiv =-^«n 


3 2 


3n3t 


fin 


(19) 


And the fin that appears in these equations could be of any magnitude; it need 
not be "small" of "infinitesimal". In fact, all we need to calculate is the 
partial derivative A^. In a way, the fin is mere decoration. In the end, 
if we wish to do so, we may multiply the resulting equations by fin and use ( 19 ) 
to recover the fiq 1 , 6*. fieZ, fiA^. During the calculation Ve will make free use 
of the fact that second partial derivatives do not depend on the order in which 

they are taken. Thus we may replace ~ ^ gf ^ by * 3q * etC ‘ In 

order to make it easier to read what is being done, we perform the calculation 
for each of the three portions of the action separately. These three portions are 


Aa,b 

atom 


r dt L x , A . 

't*a atom* a,b 
e.m. 

dt L 


r dt L , 

1 t*a e.m. 


( 20 ) 


't*a 


int * 


a. The Portion 6 A 


a,b * 
atom 


We have 


_i A 

3n a,b 
atom 


/t=a ^ 3n L atom 


In order to write the integrand in a useful form, we use the fact that 

L (t, q 1 , v 1 ) depends on n by way of the q 1 and v 1 . Thus 

atom i 

I ~ 


T 

3n natom 


z 3 \ Sv 1 . / 3 T \ M_ 

*__i ^atom 3n j^i atom 3n 


3v 
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i L atom * P i 
&v 

The v[ are the configurational momentum components. In our case, we have 

p t *' M’ v J + Ui , 

1 Li 1 

where we have used the symmetry condition M» * , which does not entail any 

loss of generality. This expression for p«, involving the primed quantities 

M’ and U!, shows why we have affixed a prime to p^,. We use 

id 1 4 i i 

3v _ __9_ 3q m 3 3q 

“ = an 3t at an ■ 


Thus 


J, 

3n atom 


. 3 ZgL 4 . /JL T ) s 

Pi ^t an + L atom 3n 

, , a^\ iai / jl l ) 2ai . 

(p i an } “ aT an + \ i atom ; an 


j* 

at 


Next we integrate over the time t. Afterwards we multiply by 6n (Note that 
6n does not depend on t.) , thus obtaining the final result 


atom 

b. The Portion 6 A 
We have 




op ^ 

5A a,b “ (p i 6ql) | b “ ^t-a dt( “3t ~ TT L atora 
t=*a 


W 


a,b* 

6«Hli 


-1 A = /5 dt L 
an a,b J t=a an e. 


m. 


e.m. 


and 


s| L ,.m. ■ /// 4T S? L 


an e.m. 


( 21 ) 


Equations (lU) yields 
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if L e. a . - - ^ • E o g - <» • E o® - <» - ^ • V s • 

where we have used equations (2) to shorten the notation. We rewrite each of 
the three terms in a way that is preliminary to an integration by parts. 

3 ,3cA _ «i . 3cA 3 


a L 

3n e«m. 


f .^Q-A , jcA o *, 

3ct 3rt 6 o 3n 3ct V 

3n o 3n o 


- ? • X e cB) - — • (5 X £ cB) . 

9n o sn o 

We integrate this equation over the volume of the room that contains the atoms 
and the radiation field. As for the first term on the right-hand side, we may 
interchange the space-integration and the time-differentiation, obtaining 


3cA 


- T"T {/// it . ,_1> . 


3cA 

3n 


By virtue of the theorem of Gauss, the contributions made by the third and the 
fifth term gets transformed into the surface integrals. 

- # « • e o E ft • 


and 


- ^ Af • x c^cB} « -<# (A? x 


3cA < 

3n 


e cS . 
o 


Now, the derivatives 


M li 


3n 3n 

that were given by equations (h), namely 

A? x * 0, and * 0 at the wall. 


and r 1 - fulfill the same kind of boundary conditions 


( 22 ) 


(Simply differentiate (h) with respect to n.) As a consequence, the two surface 
integrals vanish. Thus we get 


3 T 

3n e.m.“ “ 3ct ‘ JJJ 3q 


3 </// At 3oA 


e o E) + /// ix |i ? . e 0 E - 


- HI 4t Tn * t 5 x V s - lit • 

Next we integrate over the time t. Afterwards we multiply by fin (Note that fin 


does not depend on r and t), and obtain the final result 
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^ $X& 


* ' { c M At(Sc * ) • e o 5) 


e.m. 


t*a 


* 4a dt{ /J7 4T < «»> 9 ’ E 0 ®> ' 


- 4a 4t( /// 4 t(«=A) • [? « e 0 cB - ^ e o i]) . 


(23) 


c. The Portion 6 A 


a,b * 
int 

Here our development will become a bit more perspicuous, if we vary 
the atomic configurational parameters q* and the electromagnetic configurational 
parameters <fr and cA separately. We are not forced to do it that way, but it 
eases the task of writing and reading the equations. Accordingly we examine 
a family of histories that depend on two variational parameters n and 5. The 
history labeled by the pair {n,C) is given by the set {q*(n»t), t(c,r,t), 
cA(?,r,t)}. Then we define 6 A by 

SL| D 

int 

* A a,b * A a,b )5 " + ( 95 A a,b )5? • 


int 


int 


int 


and 


First we calculate t— A . . We have 

ol^ A j D 

int 

i Aa,b ' 4a « l£ L i„t • 

int 


3n L int = M AT 3n L int * 


We obtain t— L. . from equation (l 6). But there only J and p depend on the 
on int 


variational parameter n. So we have 

3 , 1 , ? 


3J 


± _ 


3p 
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The next task Is to express and in terms of the . Equation (5) yields 


V, * i 3J. 

3J - 3v i 

3n * J i 3n 3n v i 


j -2 iSL. + —L y 1 
3t 3n 9n 


We rewrite the first term in a way that is preliminary to an integration by 
parts. And in the second term we replace the dummy index i by J. Thus 


i 3J. 


J2U-! (J ifll) . ialZl + V J — . 

3n 3t w i 3n ; 3n 3t 3n 

Since the and J. depend on t and n by way of the q k , we get (Notes = y*^) 


&..»( j sal, + iai yJ ( !fi . lij . 

3n 3t ' J i 3n 9n V * i 


3q W 

Now we use equation (8), which followed from the principle of gene-identity. 

At the same time we can bring the and v^ behind the curl fx, since neither 
of these quantities depends on r. Thus we obtain 

We use (5) to replace v^ by J. Also, since the combination will 

occur frequently later on, we rewrite the first term in such a way that this 
combination will appear. Thus 

f* if tP ( i j t *£>].♦ ** x5] • 

Next we calculate in a similar manner. We have 

3n 


(25) 


3fi.. 


3n 


or, with (7), 


3q 


i 3n * 


iai 

3n 

depend on r. Thus 


ifi. = _ iSL v • j . 

3n 3n i 


The — may be brought behind the divergence V • , since this quantity does not 
3ti 


an 


r . / 


1 - 


Bn 


i 
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Equations (25) and. (26) may "be written in a shorter way. We multiply 
these equations by dn» The left-hand sides then become 

dJ s dri and dp = |^- dn . (27) 

On the right-hand 3ides we may bring the dn behind the differentiation operators 

3 - - - 

— , V x, v, since dn does not depend on r and t. Then there will appear the 

combination 


p. 1 3t) p i 

for which the customary abbreviation is dr. 

will be explained presently. Then, with 


9 

The physical significance of 


dr 


d? = i j^q 1 , 

equations (25) end (26) become 


( 26 ) 


dj = g|- (pdr) + V X (6? X j) , (29) 

dp = - V • (pdr) . (30) 


Another derivation of equations (29) and (30), which are used in our 
subsequent discussions, may be found in Note #2. There the derivation is . 
based on the vector calculus. 

We interrupt the process of calculating dA in order to discuss the 

^ D 

meaning of dr, called the field of virtual displacements. According to the 

principle of gene-identity, each element of charge may be tagged, i.e. it may 

be considered as an identifiable object. Let us consider one such element, 

and let a be its position in the home-space. Now a depends explicitly only on 

the atomic configurational parameters q*. When the q^ depend on the time t, 
do. 

the velocity of this charge element is given by 
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But the velocity ought to he related to the local charge density p and the local 
current density J by p * J. With the aid of (5), this equation yields 


da 

dt 


K vl 


da i 

We compare the two expressions for and use the fact that the set {v > of the 


configurational velocity components may be chosen arbitrarily. Then we find 

(31) 


3a _ 1 - 

1 ‘ 


Now let the q also depend on the variational parameters n, as before. Then 
we have 

_ 3a . 3a 3q* . _ 3a . i 
Sa = ^ an ^ in = ' 


With the aid of (3l), this equation becomes 

^ J 1 dq i . 


(32) 


An equation of this type holds not only for this particular charge element, but 
for every charge element. Of course, the p and the in (32) have to be evalua- 
ted at the position at which the charge element is situated. We may therefore 
regard (32) as a generic relation. The customary symbol for a generic 6a is 
6r, although 6a would be preferable on strictly logical grounds. With this 
change of notation, (32) becomes 

6r = |j i 6q i (28), 

i aa 5- repeated. 

According to (28), each set of {5q x = 6n) engenders a field of 

virtual displacements 6r. Additional material, especially as it concerns 

equation (31), may be found in Note #1. 

g 

We may resume the calculation of r- A . . At this point we have become so 

9 9. 9 u 

accustomed to the use of quantities which axe prefixed by a 6 that we can skip 

g 

a few steps and calculate (r— A , )6n directly. To repeat: No infinitesimals 

0 T1 SL j D 

int 
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are involved. A quantity, say Q, that is prefixed by a 5 simply stands for 
the n-derivative multiplied by fin. And fin is of finite magnitude. Further- 
more, fin does not depend on r and tj it is a constant. We have 

(tT- Jfin * /? ^ 4 . » 

fin a,b ; t*a. int 

int 


and 


Equation (24) yields 




/// Ar 5 L 


int 


SL * ^(cA • SJ - <J>cfip) . 
int c 

With the aid of (29) and (30), this equation becomes 

5l. * - {cA (pfir) + cA • [V x (fir x j)] + <f>c7 ' • (pfir)} 

int c ou 

We rewrite each term in a way that is preliminary to an integration by parts, 
and obtain 


fiL 


int 


alt (cl ’ o6?) “ * 

+ 7 V • [(fir X j) x CA] + i<V X C A) • (fir x j) + 
c c 

+ — V • [<|>cp6r] - — (V<|>) • cpfir , 
c c 

or, with the aid of equation (2), 

fiL . = (cA • pfir) + — V • [(fir x J) x cA + <J>cpfir] 
mt dec c 

+ (cpi + J x c b) * fir . 

We integrate this equation over the volume of the room that contains the atom 
and the radiation field. As for the first term, we may interchange space- 
integration and time-differentiation. The second term gets transformed into 
a surface integral by virtue of the theorem of Gauss. But this surface integral 
vanishes because of the boundary conditions (4)*, it is treated in a way that 
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1» very similar to vhat was done in the preceding section, We then obtain 

3 rfffi T _ . 1 


6L 


int 


* at {///AtA • pfir> + i ///dt(cpE + J x cB) • 6r 


Finally we integrate over the time t and obtain 


(•— A„ v ) 6 n - (///AtA • p6r) ' 4- r dt {///A.t( E + J * B )•«?}. 

3T1 f t=a t=a 

int 


(33) 


Next we calculate (-r- A -) 65 . Again we work with quantities that have 

at, a, d 


a 6 as a prefix. Here, of c&Srse,the 6 indicates the 5-derivative multiplied 
by $5. Equation (l6) yields 

The only thing we have to do here is to integrate this equation over space and 


.. . * — (J • 6cA - cp6<|>) . 
int c 


time. We obtain 


(■jj“ A a b )fiC * dt{///Ai ^<tT • 6c A - cp6<|))> . 


8? a,b' 
int 


(34) 


Finally we add equations (33) and (34). Then, according to the first 


equation of this subsection, we obtain the desired quantity 6A_ . . We arrive 

cl 9 D 

int 


m 




-•lib 


6A a,b = * / Ifo# * 


int 


t=a 


+ 1° dt{///Ax[(pE + J x B) • 6r + ^<J. • 6cA - cp6$)]} , 
t=a C 


(35) 


which is the last equation of this subsection. 


Now that we have finished expressing each of the three portions of 6 A 


a,b 


in a form that is useful for our purposes, we obtain 6 A itself by adding 

SI 9 D 


equations (21), (23), and (35). We change the order of the various terms so 
that related terms appear together. Since in the final expression the varia- 
tional parameters h and 5 no longer occur, we can replace the partial 
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3p i dp i 

derivative in (21) by the total time-derivative — jr* . (The partial deriva- 
3p ^ 

tive —st - appeared at an earlier stage of our discussion, because then the p! 

OU i 

were considered to depend not only on the time t but also on the variational 
parameter From now on we can ignore the n-dependence, since in essence it 
was only a procedural artifact in pur derivations.) We obtain 

i 


6A a,b * p i 6(1 + /// At P a * - ///A t£ q E 

dp 


m 


b 

t=a 


- P - ///ixtpE + J «5) • Jr + 

t=a 3q 


+ x r 5 -if o 


e E - J) • «A - (V • e i - p )6<p] > , 


(36) 


where we have introduced the permeability u of the vacuum, given by u = - \ . 

o o sc* 

o 

Later on, when we will state the action principle, it will be helpful if 
we replace those terms in (36) that contain the 6q* by terms in <Sr. To this 


-i 


end we introduce n vector fields s in the home-space. (Note that n was the 
number of atomic degrees of freedom) which are "reciprocal" to the n vector 


fields J^, i.e. which fulfil the relation 


///A tS 1 • Jj = . 


(37) 


Here 6^ is the Kronecker-delta: * 1 for i = J , 6* - 0 for i ^ J . The s 1 are 

not uniquely determined by (37). But our further considerations do not depend 
on the feature of uniqueness. At any rate, later on the s* will occur only in 
the products ps* ; and these are determined uniquely, because the atomic configura- 
tions are not subject to the kind of constraints that we find in engineering 
mechanics, in the form of mechanical linkages. 

We can express the Sq 1 in terms of the 6r-field if we multiply equation (28) 
by ps^ and then integrate over the volume. Then, by virtue of (37), we obtain 
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page , s 

- ///4tpl J • 5? . 0R ^ITy (38) 

We use this result to rewrite equation (36). Here the pi and the expression 


dpj 

dt 


(• 


aq. 


i L atora^ may moved bellind ' fche integration signs ///, since these 


quantities do not depend, on r. We obtain 

i 

6 A - M*"*-/-i-’ 


b “ ^//Atp(s p t + A) • fir - ///Axe o B • fiA} 


-/^ a dt(///M(%--^L aton , 


t*a 


Jpi 1 - (pE + J x B)] • fir 


+ ///At[(V x ~ B - ~ e Q E - J) * 6A - (V • ti - p 
*o 


(39) 


At this point we are ready to state the action principle. But before 
we do so, we briefly review what we have done so far. We started with an 
expression for the action A & as a functional of the history. This step was 
a part of Physics. Then there followed a rather lengthy derivation resulting 
in the expression (39) for the variation of the action. This derivation was 
mostly of a purely formal nature, except for the use of the principle of 
gene-identity [as in equation (8)], which is a part of Physics. 

The action principle, a part of Physics, states that the driver fields 
f , p , and J a are related to the variation of the action by the equation 

b 

fir - ///fi iTD • 5A} 


SA . = {///dtp 
a,b 4ii mom 


t=a 


- /+_ {dtC/J/Axi^ • fir + / // AT [ • 5A - p d 6<{>]} , 

U— fit 


(4o) 


where p is the momentum density of the atomic constituents, and where D is 
mom 

the electric displacement. 

When we subtract equation (39) from equation (40) we obtain 

r-i 


b 

t=a 


{///Ar[p mom - p(sp| + A)] • fir - ///At[D - e Q E] • fiA} 

dp * 

- /*. a dtt///4T[F a + ( s * j x BD , «? - 

- //Jat[(v x — 5 - ~ e^E - J - J^) • fiA - (V • e o E - p - p^)fi<|>]} . (4l) 

W o 
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In order to extract useful information from equation (4l), we remember 
that the atom possesses enough degrees of freedom so that there are no 
restrictions on the displacement variations 6r. (Where there are no atomic 
constituents the 6r is completely arbitrary by default.) Furthermore there 
are no restrictions on fiS and 6$, except that, at the walls of the room, these 
variations are subject to the boundary conditions 

Af x 6A = 0, and 6<J> * 0 at the wall. 

These conditions are merely a- different version of (22). In addition, the end 
points t = a and t * b of the time interval of surveillance are arbitrary. We 
observe that the left-hand side of (4l) is contributed only by the end points 
of this time interval, whereas the right-hand side is contributed only by the 
interior points. Each side is therefore a linear functional of a different 
set of quantities. Thus they can be equal to each other only if they are zero. 
Because of the arbitrariness of 6r, 6A, and 6<J> each side can vanish only if the 
factors of these quantities are zero. Thus we obtain 

P_ OTn = p(s i p' + A) , (42) 

mom I 

(43) 

(44) 

(45) 


F 4 - ( 


d p! 


D = e E , 
o * 


« - rr W ’ 51 - (pE + 3 * S) ■ 

j d a 7 X i 5 . e E - j , 
u 9t o * 


p = V • e E - p . 
o 


(46) 


r-d 


Equation (44) - (46) serve to determine the required driver fields F , 
J^, and p^. Additional information is furnished by equations (42) and (43). 
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But we will not use them. Equation (44) is of great interest. It states that 
in the absence of electromagnetic fields the required driving force density is 
given by the first term, whereas in the presence of electromagnetic fields the 
required force density is reduced by the amount pE + J KB, This means that the 
electromagnetic fields exert forces on the atomic constituents with a force 
density equal to pE + J x B. Since this agrees with the predictions of electro- 
magnetic theory, our choice of ^ irrfc » as given by equation (16), has been con- 
firmed. [it was the interaction part of the total action that produced the 
term pE + J x 8, as shown by equation (35).] It is hard to imagine - and probably 
also impossible - that any other choice for ^ would have produced the same 
result. But then equations (1*5) and (1*6) follow automatically. In terms of 
the total charge density p^ = p + p^ and the total current density * J + J d , 
these equations yield the second pair of the Maxwell equations 

V x jp- B - e Q E = J t and V • e Q E * p t , (47) 

which supplement the first pair (l). 

The most familiar application of the action principle concerns the behavior 

of a free-running system. In this case the three driver fields T* 1 , J d , and p d 

are zero. Before we set equal to zero in (44), we multiply this equation 

by ^ Jj and then integrate over the volume. With the aid of (37), we obtain 
d Pj 9 

dt “ TT L atom = - J , • (pE + J x B). (48) 

oq P j 

This equation tells us how the electromagnetic field influences the atom. 
Equations (45) and (46) yield directly 

’ * £ 5 - sf c o S ' J 
o 


and 
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These two equations tell us how the atom influences the electromagnetic field. 

Of course, in all three of these equations we have to express E and E by the 
potentials cA and $, as shown by (2). At this point we can appreciate the great 
utility of the action principle. If the action functional is formulated in such 
a way that it yields the expected equation of motion (1*8) for the forward process 
"radiation •+ atom", then it will automatically produce the equations of motion 
(1*9), (50) for the reverse process "atom -► radiation”. 

We have written down the equations of motion (1*8) - (50) for the free- 
running system in order to show that something we are familiar with results from 
the action principle. But it is more useful to work directly with the action 
principle, even in the case of a free-running system, and to leave the equations 
of motion as something that is implied by the principle. In this case the action 
principle may be stated in a rather concise form. Equation (4o) shows that, 
when the three driver fields F* , J , and p a are zero, the terms contributed by 
the interior of the time interval t * a to t = b vanish and that only the terms 
contributed by the end points t = a and t = b survive. Thus the action principle 
states « imply that 


5A 


a,b 


= E.T.O 


« » 


(51) 


where E.T.O, stands for "End Terms Only". In order to apply the principle (51), 
one first writes down the action functional A . . For the system consisting of 

8. 9 D 

an atomand a radiation field both treated according to classical physics, this 

functional is given by . 

, . dq 1 dq^ i 

A a,b = /t=a dt{ 2— M ij— +U i^t “ V + 

+ ///a T[~| (JL CA + 7$) . (—• cl + f*)-^(7x cA) • (? X CA) + 


3ct 


+ ~ (cA • J i - ij»cp)]} 




(52) 


A2 6 
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as we have seen. Then one works out the variation 6 A. "by purely mathematical 

8>|D 

manipulations and obtains (39). The next step is to set the ‘’’interior term* 
/^ -& {dt ...) equal to zero. Here one uses the fact that the variations 6r, 
and 6$ are arbitrary. The result is the equations of motion (48) - ( 50 ) . The 
end terms themselves are of no interest in the applications we have in mind. 

We shall follow the same procedure later on when we treat the atom according 
to the quantum mechanics. The only modification occurs. in the very first step^ 
namely in the formulation of the action functional. This step will be made 
easier, if we rewrite the various terms of (52) in a different oyder, namely 

A a,b * ^t=a dt( l^at M ij + [u i + //Mt • * [v ' + * 

* I!! A*t-§ =A + ?*) • cA + V$) - (7 x cA) • (? x cA)]) . (53) 


3ct 


We see that the vector potential component in the atomic configuration space 
is augmented by /// ArJ'^ • A, and that the potential energy V’ is augmented by 
///A rp$. We may regard ///AtJ^ • 5 as the ’’pull-back" of the home-space vector 
potential and /// Axp<> as the pull back of the scalar potential <|>. In the pull- 
back operation one starts with a quantity defined in the home-space and con- 
structs a related quantity in the configuration space. A third example of a 
pull-back is furnished by equation (48). In treatises an analytical dynamics 
this equation is usually written as 


21 

dt 


j , r — — ,e«nu 
8q .J atom " } 


t e.m. 


is the component of the generalized force (here of electromagnetic 


where F 

J 

origin as indicated by the superscript e.m. ). We see then that F 


e.m. 


is 


given by 
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p e.m. , jjj Ar I? . ( 5 + J x 5) (5k) 

J P J 

Thus F®* m * nay be regarded as the pull-back of the electromagnetic force density 
pE + J x B. A fourth example is furnished by equation (38), which describes the 

i — 

pull-back of a displacement field fir. A related equation would be the pull-back 
of a velocity field w from the home-space into configurational velocity components 
v , namely 

v^ * ///fitps^ • w . (55) 

The inverse of the pull-back is the ’'push-out" . In this operation we 
start with a quantity defined in the atomic configuration space and construct a 
related quantity in the home-space. We can always recover one of these two 
operations from its inverse by means of equation (37). For instance the push- 
out that is related to the pull-back (55) is 

v = ~ JjV 1 . (56) 

As we go from pull-back to push-out we use (37) to check the result ,• whereas we 
use (37) directly when we go from push-out to pull-back. Equations (28) and (38) 
constitute another pair of these two inverse operations. 

The action principle (51) in conjunction with (52) or (53) seems to be all 
that is required for the analysis of the free-running system. But when one 
comes to the equations of motion (48) - (50), one meets a difficult situation. 

For the sake of simplicity and physical reality one would like to treat the 

wm 

atomic particles as point-like objects. But E and B tend to infinity as we 
approach the object. As a consequence, the meaning of the right-hand side of 
(48) becomes obscure, as the following consideration will show. The factor ^ 
in the integrand presents no difficulties. As equation (56) shows, this factor 
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must be finite, so that w, which. here is the velocity of the particle in 
question, comes out to be finite for any arbitrary choice of the set {v 1 }. 
Integrals of the type ///A t(~ j^)p and ///a t(~ Jj)J present no difficulties, 
even though p and J become infinite at the location of the particle. The values 
of these integrals are & nd (~ ) Q v, where Q is the charge of the 

particle. But unfortunately the integrand in (1+8) contains also E and E, which 
tend to Infinity as the particle is approached. And it is this feature that 
makes the meaning of the integral obscure. As far as equations (1+9) and (50) 
are concerned, the situation is manageable, as the theory of the Li^nard-Wiechert 
potentials would show. But we will not pursue this matter, since the difficulty 
about equation (1+8) still remains. 

The discipline of quantum electrodynamics was created for the purpose of 
coping with this difficulty. But we shall not use this discipline here, because 
it is too complex. Furthermore it treats the particles relativistically, and 
the electromagnetic field quantum-mechanically, whereas we decided at the outset 
to treat the particles non-relativistically and to use the classical theory for 
the fields. Instead of quantum electrodynamics , we are going to employ a more 
elementary method. This method is based on two approximations, which will be 
explained in Section T. 

But before we deal with these approximations we will describe some simplifi- 
cations that result from the choice of a particular "gauge” for the electro- 
magnetic potentials cA and Choosing a gauge means to subject the divergence 
of cA to any requirement that will assist further developments of the analysis. 
Our ability to stipulate what 7 • cA should be is based on the following argument 
According to equations (2) the electromagnetic fields E and cB remain unchanged 
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if we subject cA and | to a gauge transformation, i'.e., if we replace cA and $ 

- - 3 

by cA + Vx and $ - — x respectively^where x is an arbitrary function of the 
position r and the time t. There is then some redundancy in cA and which we 
can exploit to impose conditions on ? 4 cA. The most common conditions are 

V • c J + <p * 0 , ( 57 ) 

which embodies the choice of the Lorentz gauge, and 

? • cA * 0 , (58) 


which embodies the choice of the Coulomb or radiation gauge. We choose the 
latter. The consequences of this choice will be described in the next section. 

Besides the constraint. ( 58 ) we will also impose a similar constraint on the 
variation 6cA, namely 

V • (6cA) = 0 (59) 

The statement that we are going to use the Coulomb gauge is not quite precise 
enough. What we are really going to do is to expand the vector field cA in terms 
of longitudinal and transverse cavity modes. And then we impose the constraint 
that cA be a linear combination of only the transverse modes. In other words: 

Those expansion coefficients that belong to the longitudinal modes are set equal 
to zero. The equations (58) and (59) still hold. But cA and <5cA when constrained 
in this manner will have additional properties. For instance, the boundary conditions 
Af x cA = 0 and Af * 6cA = 0 (60) 

will be fulfilled. [See (4) and (22)3. All of this will be explained in the 
next section. Instead of saying that we will subject the vector potential c5 
to the Coulomb gauge (58), we will say that we constrain cA to be "transverse”. 
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We shall base our discussions on the formulism that way he found in the 
hook "Microwave Kleetronies" by John G» Slater, (Van Noatrnml 1950), As this 
hook may he out of print we shall repeat the arguments that are presented there 
hut sometimes in modified form. We shall use the word "cavity" , customary in 
this context , for the room that contains the atom and the radiation field, 

. 'this cavity is assumed to have perfectly conducting walls, which implies that 
the electric field 1 and the magnetic induction cB satisfy 1 the boundary condi- 
tions 

hf k J? »* 0 , and A? • cB * 0 (a) repeated 

Slater ' a formalism pertains to the expansion of fields in terms of cavity modes, 
Before we get started we review three theorems of geometry, The first 
theorem states that the surface integral of the curl of an arbitrary smooth 
vector field V is aero, when this integral is taken over a closed surface. 

That is, we have 

'Sf./a 4? « ? * 0 ■ 0 , (dl) 

Xn order to prove (^l) we transform the surface integral into a volume integral 
over the enclosed space, hy virtue of the theorem of Gauss, Thus 

A * • F * V w /// At V ■ V * V . 

Gince the divergence of a curl is aero, equation (6l) follows. 

The second theorem states that, for a scalar function \J) which has the 
constant value G on some smooth surface F, we have 

hf * $\|> « 0, if d* * 0 on F, (6g) 

i.e. the gradient is normal to the surface. This should be clear, when the 
geometrical situation is visualised, A more formal proof goes as follows. We 
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can always extend the surface, so that it forms a part of a closed surface. Equa- 
tion (6?) is equivalent to the statement that 

^(Af x ty) • w * 0 

for any arbitrary smooth vector field w. And this statement results from the 
following calculation. 

•£/>(A f x 7i|>) • w * </-/Af • (Vt{i x w) s ff&£ • [7 x (||>w) - lf*V x w] . 

The integral of the first term is zero by virtue of (6l). The integral of the 
second term vanishes, because * 0 on the surface. Since 7i|/ is not a.ffected, 
when we add a constant to 1 ft, the statement ( 6 2) may be generalized to read 

Af x fi(i 5 0, if ip = const on F. (63) 

The third theorem states that, for a vector field E, which is normal td a 
smooth surface , i.e. for which Af x E * 0, the curl is tangential to the sur- 
face, so that 

Af • 7 X E - 0, if Af X E - 0 on F. (6k) 

In order to prove ( 6k ) , we again extend the surface, so that it forms a part of 
a closed surface.' 1 Equation (6^) is equivalent to the statement that 

4$ Af • (V x E)w = 0 

for any arbitrary smooth scalar field w. And this statement results from the 
following calculation. 

<#Af • (7 X E)w = £ #Af • [7 x (Ew) - (7w) x E] 

= ff>bf • 7 x (Ew) +c#Af x Ew . 

The first term vanishes by virtue of (6l). The second term vanishes because of 
the assumption Af x I = 0.* 

« 

We assume that i|» * 0 on the extended surface. 

+ _ — 

We assume that Af * E = 0 on the extended surface. 
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We now consider the following two eigenvalue-eigenfunction problems; 

7($ • F ) + » 0, with the boundary condition $ • F =0, ( 65 ) 

a a a a 

and 

) + m 2 G ■ 0, with the boundary condition Af • ff * 0. (66) 

a. a, a 3* 

Here the i 2 and m 2 are the (positive) eigenvalues, and the P and G are the 

El Gl a cl 

vectorial eigenfunctions. These eigenfunctions are irrotational (i.e. the curls 
vanish), as one can see when one takes the curl of ( 65 ) and (66). (Note that 
the curl of a gradient is zero.) The index "a" labels the eigenvalues and 
eigenfunctions jthere is an infinite number of them. We assign the indices "a" 
in such a way that the eigenvalues are ordered according to their magnitude, 30 
that 

i-! 2 1 l ?. 2 1 l 3 2 1 • • • » 
mi 2 <_ m 2 2 <. m 3 2 <. . . . . 

In order to make sure that the problems (6 5) and (66) conform to the 
classical paradigm of such problems, we have to prove that the operator 77»is 
self-adjoint. The self-ad Jointness conditions are: 

/// A:t [f. v( 7 • F' ) - F* • 7(7 • P) ] = 0 for two functions F and F* 
that satisfy the boundation conditions 7 • F = 0 and 7 • F* - 0, ( 67 ) 

and 

///At[G • V C V • G') - G’ 7(7 • G )] = 0 for two functions G and G' 
that satisfy the boundary conditions Af • G = 0 and Af • G 1 ~ 0. (68) 

To prove these relations we note that 

F • 7(7 • F») = 7 • [F7 • F' ] - (7 • F) (7 • F' ) . 

Similarly 

F’»7(7 • F) = 7 • [F*7 • F] - (7 • F’)(V • F) . 
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On subtracting, we obtain 

F * 7($ • ?’) - F’ • 7(7 • F) =7 • [FV * F' - F'V * F] . 

We integrate this equation over the cavity and transform the right-hand side 
into a surface integral with the aid of the theorem of Gauss. Thus 
///At[F *7(7 • F* ) - F* • ^(7 • F) ] = “/^Af • (F7 • F’ - F»7 • F] . 

But the surface integral vanishes because of the boundary conditions that are 
stated in (6j). Hence (67) is established. We prove (68) in the same way. The 
only change is that here the surface integral vanishes because of Af • G * 0 
and Af • G* = 0. 

Next we prove the orthogonality relations 

///AtF • F, * 0, for l 2 f l 2 , ( 69 ) 

fit D 8. D 

and 


/// AxG a • 3^ * 0, for t m* . (70) 

We take the scalar product of ( 65 ) with F^ and obtain 

a 2 ? • F = - F • V(V • F ) . 

a a b b a 


Similarly , 

• F = - F • 7(7 *> F, ) . 
bo a . a b 

We subtract and integrate over the cavity. Then, with the aid of the self- 
adjointness relation ( 67 ), we obtain 

*g)//Ji • F b - 0 . 

And, since we assumed Si 2 - l 2 ^ 0 t equation ( 69 ) follows. Equation (70) is 

a b 

proved in the same way. 

So far we have proved ( 69 ) under the assumption that ^ But when an 

eigenvalue is degenerate, i.e. when it admits several different eigenfunctions, 
we may always redefine these functions (by taking suitable linear combinations) 
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in such a way that (69) still holds, even when A 2 * A 2 . Similar comments apply 

fit D 

to equation (TO), Finally we normalize the eigenfunctions so that 



/// 4Tp a • p a * 1 ’ 

(71) 

and 

UM a ’ % * 1 • 

(72) 

On combining ( 69 ) 

and (71), we obtain the orthonormality relation 



/// ATf a - f b * S ab • 

(73) 

Similarly 

W ** a ' B b * S ab • 
* 

(7^) 

Here, 6 . is the 
at> 

Kronecker delta: 6^ = 0 for a ^ b, 5^ 38 1 for a = b. 

There 


is no orthonormality relation between an F and a G^. 


W 3 can get a better insight into the behavior of the vectorial eigenfunc- 
tions F and G , if we relate them to scalar functions \p and x defined by 
a a a a 

A J n * - V • F . (75) 

£1 E SI 


and 


m x = - V • G 
a a a 


(76) 


Here A and m are the positive square roots of the eigenvalues A 2 and m 2 . 

& & a a 

By direct substitution of these relations into the first terms of (6 5) and (66), 
we obtain the inverses of (75) and (76), namely 


and 


A F = Vtp , 
a a a 


m n G a = V X Q ‘ 
a a a 


(77) 


(78) 


Inspection of (75) and ( 78 ) in conjunction with the 

boundary conditions for the F and G stated in ( 65 ) and (66) yields the 

a a 

following boundary conditions for ip and x • 

a a 
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* 0 at the wall, 

E 


(79) 


Af • * 0 at the wall (8o) 

A 

On combining ( 6 2), (77) and (79) we f/ee that the f automatically satisfy also 

E 

the boundary condition 

Af x p = 0 , (8l) 

E 

i.e. the F are normal to the wall, 
a 

When we take the divergence of ( 65 ) and (66), divide by - A and - m 

E E 

respectively, use (75) and ( 76 ), and adjoin the boundary conditions (79) and 
( 80 ), we see that the and x could have been determined also directly from 
the following two eigenvalue-eigenfunction problems: 

^ • Vip + l z 4> =0, with the boundary condition i|i * 0, (82) 

E E E E 

and 

9 • Vx + m 2 x = 0, with the boundary condition ; • Vx * 0 . ( 83 ) 

E E E E 

Here the self-adjoint-ness relations for the operator V • V are: 

///At[(|<^ • T V • Vt|»] = 0 for two functions ^ and 4*' 

that satisfy the boundary conditions ip * 0 and ^' = 0 , (84) 

and 

///At[xV • Vx* - x’V • V X 1 * 0 for two functions x and x f 
that satisfy the boundary conditions Af • Vx = 0 and Af • Vx’ = 0 ( 85 ) 

These relations are established by the same technique that was used in the 
proof of ( 67 ) and (68). One starts with the formula 

tj/V • Vtp * - ij/’V • Vip = V • [ipVtp* - i|>’ViJ>] . 

Consequences of the orthonormality by relations (73) and (74) are the 
following orthonormality relations for the <J> and x J 

E E 
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fff^A * 5 ab’ 


and 


/// 4T Vb = S ab' 

In order to prove (86), we use (75) and arrive at 

= (V « Fj(V • F v ) * V • [F a ($ • Fj] - F a • V(V 


( 86 ) 


(87) 


Va^b " w * a' ' ’ ‘b 
or, with the aid of ( 65 ), 


V - 


= v • [F ( V • f )] + £*F . p , 

We integrate this equation ever the cavity. By virtue of the theorem of Gauss, 
the first term gets transformed into a surface integral, which vanishes because 
of the boundary condition V • F^ * 0. Then (86) fellows. Equation ( 87 ) is 
proved in the same way. 

According to established eigenvalue-eigenfunction theory the set } is 
complete, i.e. any function ^ that occurs in practice may be expanded in terms 


of the , so that we may write 
& 


* = I e * 


a=l 


a a 


( 88 ) 


We can determine the expansion coefficients 0 in the following manner. We 
multiply (88) by ^ and integrate over the cavity. It is permissible to use 
term-by-term integration on the right-hand side.* Then the orthonormal. ity 
relation (86) yields 

- ///A . (89) 

Similarly the set (X } is complete, if we adjoin. the (constant) function 

8 * 

x a „ .A- where Vol is the volume of the cavity. Here we have to supple- ... 

0 

ment ( 87 ) by the orthonormality relation 

J/M 0 x b ■ « ob , 


( 90 ) 


Term-by-term integration is permissible when the expansion (88) converges in 
the mean. See equations (107) - (ill), which pertain to the analogous case 
of vector fields. 
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which is easily proved. For b * 0, (90) is obvious. For b t 0, we have, 
with (83), 

Vb * - it? *o’ • \ • 

Then (x q * const), with the theorem of Gauss, 

///4tx 0 x b - - • V 

But the surface integral vanishes because of the boundary condition Af • $x fe * 0. 
The eigenvalue m Q 2 is, of course, zero. 

The convergence of the expansion series is most rapid, when the function 4* 
that is to be expanded satisfies the same boundary conditions that the eigen- 
functions satisfy. Thus we would prefer'to expand a function 4* with the boundary 

condition 4< a 0 (alternatively Af • V4> * 0) in terms of the 4* (alternatively x_)* 

a a 

Any irrotational vector field may be expanded in terms of the eigenfunctions 

F , or in terms of the eigenfunctions G . But we are not going to pursue this 
a a 

matter yet, because we are interested in the expansion of general vector fields, 

not Just of irrotational ones. For general fields, the sets {F } or {G } are 

a a 

not complete. Therefore we must construct • additional sets of vectorial eigen- 
functions. There will be two sets {E } and (K } of solenoidal eigenfunctions. 

a Ci 

('Solenoidal* means: the divergence is zero.) They arise from the eigenvalue- 
eigenfunction problem 

V k (V x 1 ) - k 2 E = 0, with the boundary condition A? * E * 0, (91) 

a a a a 

and 

V * (V x H ) - k 2 H = 0, with the boundary condition Af x (V x H ) =0. (92) 

a a a a 

The (positive) eigenvalues k 2 are the same in both problems, as we shall see 

a 

shortly. Again we assign the indices "a" in such a way that the eigenvalues are 
ordered by magnitude, so that 
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ki 2 <, k2 2 <. k3 2 < ... . 

That the E and H are solenoidal is seen easily when we take the divergence 
a a 

of (91) and (92) and observe that the divergence of a curl is zero*. Here the 
self-adjointness relations are 

///At[E • 9 * (7 k!') • E 1 • $ x (7 x E)] * 0 for two functions E and £' 

that satisfy the boundary conditions if x I * 0 and if x E' * 0 > (93) 

and 

• 7 * (7 x H’) - H' • 7 x (V x H)] » 0 for two functions fi 
and B' that satisfy the boundary conditions if x (7 x H) * 0 and 
Af x (7 x h* ) = 0. (94) 

In order to prove (93), we note that 

E • 7 x (7 x E») * 7 • [(7 x S») x E] + (7 x E* ) • (7 x E) . 

Similarly 

1* • 7 x (7 x E) = 7 • [(7 x E) x E'] + (7 x E) • (9 x £' ) . 

Subtraction, integration, and the theorem of Gauss yield 
// / At[E • 7 x (V x E*) - E’ • 7 x (7 x E)] 

= <^Af • [(V x S') x E - (V x E) x E' ] . 

But the surface integral vanishes because of the boundary conditions. The 
relation (94) is proved in the same way. 

Now we show that the eigenvalues k 2 are the same in the two eigenvalue- 

E 

eigenfunction problems (91) and (92). Instead of determining the H by way of 

(92), we can obtain them also through the definition 

ki = 7 x g (95) 

a a a 

from the E , Here k denotes the positive square root of k 2 . If we substitute 

E E E 

(95) in the first term of (91) we obtain the inverse of (95), namely 
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k E * V x R , 

61 ft fll 


(9«) 


Then, when we take the curl of (91) (and also divide by k )» and use (95) and 

Si 

(9 6), we see that (92) is a direct consequence of (91). 

Equation (95), in conjunction with the theorem (6L), tells us that the 3 

Cl> 

automatically satisfy the boundary condition 

Af • 3 s 0 . (97) 

a 

Just as we established the orthaiormallty by relations (73) for the f , we 

6L 

obtain the orthonormality relations 

///AtB a • \ - « ab • (98) 

All we need to do, is to substitute E & for E and E^ for E' in the self-adjoint- 
ness relation (93) and to use equ&tlon (91). The use of the formula 
k a B a • k^ * (7 x E a ). • (? xys? . (v *1^)3 + 

+ E a • V x (7 x ^ = 7 . [E a x (V X E^)] + ^ , 

integration, the theorem of Gauss, and the boundary condition Af *1 * 0, 
give us the orthonormality relations 

UI^K ■ \ - « ab <*9) 

as a direct consequence of (96). There is no orthonormality relation between 
an J$ a and an H^. 

Finally we prove the orthogonality relations 

///AtF a • - 0 , (100) 

///AtG a • ^ = 0 . (101) 

and 

/J/AT? a .« - 0 . (102) 

For we have 


f . • \\ ' \ ’ * ®b * V * * V + *b • ? * f a • 

S « ’ ’‘A * K ' ’ * ®b * ? ‘ (E b * V * ®b • ’ * 5 a • 


A^O 


\ ' k b H b ' P a * ’ * \ "■ ? • (E b * V * ®b • 5 * f a • 

But V x P * 0 and 9x3*0. Then* integration and the theorem of Gauss give 

ft Cl 

us 

k b ///ArP a • ^ * #A f • Hh x , 

\III^ 3 a • Bh * $Af • Eh x G a . 


and 

♦ 

V//AtV ®b s * (E b * V 

But the surface integrals vanish, because f and are parallel to A?. 

Before we continue with our discussions, we summarize the most important 
properties of the eigenfunctions in boxes, in order to facilitate the task 
of referring to an equation. 

Box #1. Eigenvalue-eigenfunction problems. (BC stands for boundary condition.) 


9(9 • P ) + SL 2 F * 0, BC: 9 • P * 0 , 

a a a a 

(65) • 

9*9i|> + i 2 •<!» * 0, BC: ill * 0 . 

a a a * a 

(82) 

9(9 • 3 ) + m 2 : G * 0, BC: A? • 3 * 0 . 
a' a a * a 

(66) 

9 • Vx + m 2 x * 0, BC:. Af • Vx * 0 . 
a a a a 

(83) 

9 x (9 x E- ) - k 2 ff = 0, BC: Af x E = 0 . 
a a a a 

(91) 

9 * (9 x H ) - k2 •£ = 0, BC: Af x (9 x B ) e O’ . 

& & 3> Q, 

(92) 


We note that there is a duality in the BC ? s for each pair. One BC involves the 
function itself, whereas the other BC involves a derivative. 



Relations between eigenfunctions 


a a a 

# P a _ 

a 


* m 3 

A a a a 


Box #3. Additional boundary conditions. 


if x f = 0 . 
a 

( 81 ) 

Af -x Vip * 0 . 
a 

(77 + 81) 

A? • a * 0 

a 

(97) 

Af • (7 x E ) = 0 

E 

(95 + 97) 


Box #4. Properties of derivatives, as obtained from Box n2. 


? X f s 0, 5 K 3 * 0, $ • H * 0, $ • E 
a * a * a a 







Box #5. Orthonormality relations. 
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/// ATf a • 'b * 4 ab • 

(73.) 

* 4 ab ‘ 

(86) 

///ATg a • 0 b - « ab 

(74) 

///A T X a X b - « ab 

(87) 

///Ati a • ^ > « ab 

(98) 

///ATH a • \ > 6 ab 

(99) 

J//AT? a • \ = 0 

(100) 

/// AtC a * ^ ” 0 

( 101 ) 

II M a • *b * 0 

(102) 


Ak2 


So far our treatment of the eigenfunctions was strictly formal,! But they 
allow a direct physical interpretation. For instance, a possible free-running 
(no charges, no currents, except in the walls) electromagnetic field in the 
cavity is given by 


E(r,t) = C ! (r)sin k ct , 
a a a 


(103) 


cB('r,t) “ C 3 (?) cos k ct , 
a a a 

where C is a constant. One easily checks that the Maxwell equations (1), 
a 

(1*7), (for = 0, * 0) and the boundary conditions (3) are fulfilled. In 

this case the fields oscillate in the single (normal) mode #a. Therefore the 

functions E (r) and H (r) are called the electric and magnetic part of the 
a a 

mode pattern #a. The other eigenfunctions, namely the pair (r), ? (r) and 
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the pair x (r), G (r), do not permit such an interpretation in terms of electro- 

fit €L 

magnetic oscillations. But they are associated in a similar vay with acoustic 
vibrations in an air-filled cavity. Thus, the pressure excursion p and the 
velocity v in a cavity with ’'hard 1 ' walls, which impose the boundary condition 
A? • v * 0 , are given by 

P ■ c a X a sin m ft ct , 

v = C — G cos' m ct' , (104) 

a sc a a * 

for mode Here, of course, c is the speed of sound, and s is the ambient 

mass density of the air. One easily checks that (lOU) is consistent with the 
first-order (for vibrations of small amplitude) equations of motion and 
continuity, namely 

S V = - V P , 

~ p = - SC 2 V • V , (105) 

where we have used the relation B = sc 2 for the bulk modulus B = s -j- . We 

may also contemplate a cavity with "soft" walls, which impose the boundary 

conditions p = 0, even though it is hard to envision how such walls might be 

constructed. (A liquid drop may be one way in which such a soft-walled ''cavity’' 

can be realized. ) In this case, the motion for a single mode #a is given by 

p = C ifi sin m ct 
* a a a 

v = C — F cos m Ct. (106) 

a sc a a 

We have related the modes with the patterns (E , H ) to electromagnetic 

si a. 

oscillations and the two classes of modes with the patterns (x_,G ) and 

a a 

(t|/ ,F ) to acoustic oscillations. Since plane electromagnetic waves are trans- 
a a 

verse, and since plane acoustic waves are longitudinal, one uses the term 
"transverse modes" for the (E ,R ) - modes and the term "longitudinal modes" 
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for the (x ,&) - modes and for the (tfc ,F ) - modes. This terminology is 
u & a a 

somewhat inaccurate, because the words "transverse" and "longitudinal" make 
sense only for plane waves, which are characterized by a definite direction. 

No such direction can be assigned to the complex patterns of the cavity modes. 
Nevertheless, these two words are commonly used, because they simplify the 
language . 

We accept, without giving a proof here, that the two sets {E } and {F > 

a a 

taken together are complete. This means that any vector field V in the cavity 

that occurs in practice may be expanded in terms of the E and F , so that we 

a a 

can write 

00 00 

v = l c-F + [ o s . ( 107 ) 

a=l a=l 

Here, the C’ and C are the expansion coefficients. The first series in (107) 
a a 

is called the longitudinal part of V and is denoted by V^, the second series is 
called the transverse part of V and is denoted by V T . The series in (107) are 
"convergent in the mean". Thus, if we approximate V by the first N terms, i.e. 
if we replace it by 

N N 

v . = y c'f + y ce , (io8) 

appr,N a a z '. a a * 

a=l a=l * 

has the property that the integral of the square of its 


the error V - V .. 

appr ,N 

magnitude tends to zero, as N tends to infinity, so that 


lim IH At(V - V appi . jH ) • (V - V_, M ) = 0 


N 


appr ,N 


(109) 


The expansion coefficients C ’ and 0 may be determined from V in terms of 

St St 

integrals, as shown by the following argument. We start with the inequality 
of Schwarz 

[// /AtU • W] 2 < [///drO • U][///AtW • W] , (110) 


Alt 5 


which holds for any pair U,W of vectorial functions. Let us make the choice 

U » V - V „ and W * P . Then, with the normalization condition 

appr ,N b 

///AT-? b . F b = 1 

we obtain 


mis 
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[// /at(v - V ) 
appr,N 


< //M(v-f J • (?-? 


,) 


V — 'apprjN' ' ’ 'appr,N' 

Wow let us take the limit for N + ®. According to (109), the right-hand side 

tends to zero. Thus we get 

/I/*’ %p r , W • F b = ///A ? . F b 


But as soon as the index N exceeds the index b, the integral on the left assumes 
the value C^ f , by the orthonormality relations. Thus the limit is also equal to 
C^ f , and we obtain 

<y - ///* tv • P b 

Similarly, 


C b * IfjL tv • Eh . (Ill) 

The result (ill) agrees with the following procedure: Take the scalar product 

of (107) with 5^ (or 1^), integrate term-by-term, and use the orthonormality 
relations. Now that we have established the legitimacy of term-by-term inte- 
gration, we can use it freely in our future work. 

Let us assume that the vector field V is of sufficient smoothness so that 
the derivatives V * V and V x V exist and that these derivatives can be expanded 
in terms of eigenfunctions. We choose the types of eigenfunctions in the expan- 
sion for V • V and 7 * V in such a way that the new expansion coefficients can 

be easily expressed in terms of the coefficiet .1 C f and C , which occurred in 

a a 

the expansion (107) for ?. Accordingly we write 


? • V 


l K L\ 

a 


and 


At(6 
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Indeed, we have 


X V « J D' 5 + J D H 
L a a ** a a 
a a 


(113) 


K» * / / /Atip V • V - ///Vt[? • (ifi 7) - V • ] 

a '''a .a a 

= #4? • * a V - « a ///4rV • F a . 

The surface integral vanishes because of the boundary condition ip * 0. Then, 

& 

with (ill), we obtain 

K 1 * - t C 1 , (ll4) 

a a a 

so that 

VTV^-Tc’tip . (115) 

a a a 
a 

We would have obtained the same result (115), if we had taken the divergence 

of (107) term-by-term, since 

V • F = - l ill and V • E . - 0 . 
a a T a a 

Next we calculate the, expansion coefficients D. ’ and D in (113). We have 

a a 

D’ = ///AtG • V x v = /// Ar[V • (V x. G ) + V • V x 5] * 

fi> Q. El ® 

a *$Af • v X G + ///At' ? 1 f x 5 a 

= -f/Af • V X (J + o , 

since V x * 0. Similarly 

D = ///AtR • V X v = /.jAf • V X f{ + ///AtV • V x H 

= £fAf • V x H + k C , 

0 0 0 

where we have used the relation V x H * k E and the second formula (ill). On 

a a a 

inserting these results for D and D* into (113) we obtain 

0 0 

V X V a HffA f • V x S ]3 + J[^A? • V X R + k C ]H . 

0 0 0 0 0 0 


(116) 


PAGE is a 
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Had we taken the curl of (113) term-by-term, we would have obtained a different 
result. For we would have lost the surface integrals that appear in the expan- 
sion coefficients of (116). However, if V happens to satisfy the boundary 
condition 

i? x V * 0 , (117) 

then these surface integrals vanish. In this case then, we are permitted to 
evaluate the curl term-by-term. 

If the series (107), (112), (113) for the fields V, V • and V x V were 
uniformly convergent, then we could deduce (by the usual continuity arguments) 
that these fields satisfy the same boundary conditions that the eigenfunctions 
satisfy. Thus, for uniform convergence, we would have 

Af x y = o, V • V = 0, Af • (? x I) = 0, at the wall. (118) 

Unfortunately, the proofs for uniform convergence, at any rate those that are 
most easily carried out, start with conditions that are not satisfied in the 
application we wish to make. (What foils these proofs is the point-like 
nature of the atomic particles . ) However , later on we are going to make an 
approximation in the Lagrangian. As a result of this approximation, the 
infinite series will be replaced by finite series, so that all the convergence 
problems will disappear. And then the boundary conditions (118) are certainly 
satisfied. 

So far we have considered a general vector field V and its derivatives 
V • V and V x v. Now we apply what we have seen to a specific case, the vector 
potential cA. In analogy to (107) we write 

cA * I Q' P + I Q E 
a a a a 

a a 


ORIGINAL 

OF POOR 


9 


(119) 


ORIGINAL PAGE IS 
0F P0 °R QUALITY 

where the Q f and Q are the expansion coefficients. And in analogy to (118), 

EL Cl 

we have the boundary conditions 

A? x CA = 0, V • cA = 0, A? • {$ x cl) = 0, at the wall. (120) 

The third of these conditions is a consequence of the first, as we have seen 
in the proof for (61*). Again we remind ourselves of the circumstance that, 
in the end, the series (119) will have only a finite number of terms, so that 
we will not be bothered by convergence problems. [See the remarks following 
(118)3. 

Now we are ready to formulate the gauge condition for the vector poten- 
tial cA. We require that all the primed expansion coefficients Q* be zero. 

8 . 

Then the longitudinal part of cA, i.e. the. first series in (119) vanishes, and 
only the transverse part, i.e., the second series, remains. Our gauge condition 
is then that the vector potential cl is purely transverse. As a consequence 
[see (115)3 the divergence of cA vanishes everywhere. Thus we have 

V t cA = 0 , (121) 

But, in addition, the boundary conditions (120) are satisfied. The second of 
these is already implied by (121). Thus 

Af « cA = 0, Af • V x cA = 0, at the wall. (122) 

Similar remarks apply to the variation 6cA: Merely insert the symbol 6 as a 
prefix to cA, o* , and Q . 

& £L 

We summarize the discussion of this long section by the statement that we 
choose the "transverse gauge" for the vector potential. This term is more 
restrictive than the term "Coulomb gauge". The latter implies only (l2l). But 
the former implies that cA can be expanded as a series in the eigenfunctions E 

EL 

alone, the F are not needed. A particular consequence is that the boundary 
8 
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(122) are satisfied. 

In the next section we will introduce the approximations we have already* 

* 

alluded to. They will make it possible to replace infinite series by finite 
ones. 

7. Two Approximations 

We start with the expression (53) for the action functional* which we 
write down again for the sake of easier reference. 

Mij SaL + [ V + ///irjj . I]4ai . [ V . + ///atp*] 


+ ///a t[-| (■— r cA + V4) • (•“■ cA + $4) 1 (V X CA) • (V x cA)]} . 

(53 repeated) 

The functions that are to be obtained from the action principle are q*(t), 
4(7, t), and cA(r,t). Since we agreed to use the transverse gauge for the 
vector potential cA, we can simplify (53). For, when we multiply out the 
scalar product that contains V4, the following sum of three integrals appears 


in (53). 

///AT • V* ♦ ///At <3^ eS) Oj) ♦ ///At ^ 

The last of these vanishes because of our choice of the transverse gauge for cA. 


cA . 


For we have 
///AT£ o V* 


9 

9ct 


cA = e o ///Ar[V • (4 7 ^ cA) - 4 


9ct 


cA] , 


where we have interchanged and V in the last term. But since V • cA 8 0, 


we obtain 


///Ate 0 ^4 3ct — - •grr-- t 3ct 


cA , 


^ cA - e/^Af • 4 

when we convert the remaining volume integral into a surface integral with 


the aid of the theorem of Gauss. But this surface integral vanishes because 
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of the boundary condition $ * 0. As a result, the expression (53) is 


A50 


simplified to 


A * /” dt{| ^ M’ ^ + + UU^i - 5] - [V + ///ATP*] 


tip" 2 dt « 

e 

+ ///A t[-~ • v<p + 


p (T~ CA) • (~r c5) - 4 (J X cX) • ($ X CA)]} . (123) 


2 TY Y ' 2 3ct 3ct w# _ 2 

We see that the scalar potential $ and the vector potential cA have become 

"decoupled", since terms that are mixed in <p and cA no longer occur in (123). 

* 

One of the consequences of the action principle is equation (1*6), namely 

d - 


7 • e Q E - p . 


(46) repeated 


(it results, when we examine <5A . under the assumption that only <J> is varied.) 

ft f D 

For a free-running system, to which we confine our attention here, the driver 

i 

field p d is zero. Thus (1*6) yields 

^ • e Q E - p . 


And with 


E = - V* - cA, V • cA = ° , 

we obtain 

e Q 7 • 7<J> = - p. (124) 

This equation shows that <f> is equal to the electrostatic potential associated 
with the atomic charge density p. Since p depends explicitly only on the con- 
figuration {q*} and the position r, the scalar potential <j> will depend explicitly 

only on these variables. Then the two integrals 

e 

J7/Atp 4> and — ^ ./// At 7<J> • 7<f> , 

which appear in (123), will depend only on {q*}, since the r-dependence of 
and p was wiped out by the integration process. These two integrals are 
closely related. Indeed with (12U), we have 
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/ / /Ai:p* ■ - c o ///At* 5 • ?* * - e Q ///AT[V • (*V*) - V* • V*] , 

The integral over the first term is converted into the surface integral 
<$)Af • *?*, which vanishes because of the boundary condition * * 0. Thus 

///At P * .^///At?* • 9* . ( 125 ) 

We use this result to simplify (123). With the abbreviation 

V(q i ) = V’(q A ) + -f ///ATfy • V* , (126) 

equation (123) becomes 

\ b * ^ dt(~ M 1 + fn » + ffUr 3 • A] - V + 
a,D t*a l 2 dt ij dt lU i JJJ * d i AJ dt 

+ !IM - 1 CA) • (^- cl) - \ (? « CA) • (7 « CA)]} . (I2T) 

We note that the quantity V depends explicitly only on the atomic configuration 
{q*}, and that the scalar potential * has disappeared from the action functional. 

The procedure that enabled us to eliminate the scalar potential * from the 
action functional may be applied also to the vector potential cA. In order to 
get ready for it, we use the eigenfunction expansion (119). Since we are using 
the transverse gauge, in which all the coefficients Q f are zero, the expansion 
simplifies to 


ca = y q e . 

« n o 

(128) 

cL a, 

a 


Similarly, 


3 r.lff = 

3ct c dt a * 

a 

(129) 

and, because term-by-term differentiation is permissible , 


V x cA = [k Q H 
u a a a 

(130) 


Next we express the three volume integrals in (127) that contain the vector 
potential in terms of the expansion coefficients Q . We have 


4 
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e * . c dQ dQ. 

///* T -i (ilt •*> • 'ilt '*> ■ UMl it V • <z tK> 

c dQ d& a e dQ D 


(131) 


a b ' a 

where we have used the orthonormality relations (98). (See Box #5.) Similarly 

///At ^ (V x cX) • (V x cA) « % I (k ft Q a )2 . 


In the remaining integral we expand the functions in terms of the and E a , 
i.e. we write 

J 1 * l J i, + l J i,a E a • (l32 > 

a a 

The expansion coefficients J’ and J, depend on the atomic configuration 

1 )8* I jQi 

{q^}. The first series in (132) is the longitudinal part ^ of J^, the 

second series is the transverse part J. _. Then, with the aid of the ortho- 

1,1 

normality relations (98), (lOO) (See Box it 5) we obtain 

• 5 * | I J l,a°'a • 
a 

When we insert (131) - (133) into (127), we obtain 

^a b = ^ dt^ M* +[U'^y J 0] — ^ 

a * t=a Utt 2 dt , dt * lU i c L J i,aV dt 


(133) 


V + 


iJ. 


* b l hk> 2 - ■§ KW 2 ’ 

a a 


(134) 


(We have used the same symbol "a" for two different purposes. It indicates 

not only the lower limit in the time- integration, but also the mode-numbers. 

This double usage should not cause any confusion, since the place where the "a" 

is written makes its meaning unambiguous.) The functions that are to be determined 

from the action principle are the n functions q*(t) and the infinite set {Q (t)>. 

a 

In order to derive the equation of motion for one of the Q & (t), we examine 
the variation of the action under the assumption that only this particular Q 


is varied. With the aid of the calculation 
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. dq dQ dQ dQ , , dQ d 2 Q 

4 * Tt s h£> ■ 


we obtain 

c dQ 


d 2 Q 


- J ti *» u Iflj 

SA a,t> ■ SThft 4 V + J a dtt c J l,a "3t - ^ Itf - c o (k » ,2<5 a 1S «a ' 

t-a * a 

The action principle 


(135) 


SA 


a,b 


E.T.O. 


(51) repeated 


for a free-running system causes the square bracket in the integral of (135) 
to vanish. Thus 

d 2 Q 

-dtf ♦ - 0 


. i ’ 


c Q i» a dt • 

This equation shows that the expansion coefficient Q (t) behaves like the 

cl 


(136) 


- i 
*3^ 


excursion of a harmonic oscillator that is driven by the ’’force — - J. — jx . 

» i,a at 

0 


The resonance frequency ft of this oscillator is given by 

a 

ft * ck 
a a 


(137) 


Equation (136) suggests the kind of approximation we are going to make. 

We will first go through the procedure of this approximation. Only later will 

we discuss the motivation for it. At this point we merely mention that it 

serves to eliminate infinities that arise when the atomic particles are treated 

a3 £oint-like entities. We divide the modes into two groups: the ’’low” modes 

whose resonance frequencies lie below some critical frequency ft^^, which will 

be selected later, and the ’’high’-' modes whose resonance frequencies lie above 

ft ... For the low modes, we retain equation (136). But in the high modes, 
crit d 2 Q 

the "intertia term” is small when it is compared to the "stiffness term" 

(ck ) 2 Q . We may therefore neglect the inertia term. Accordingly, we re- 
a a 

place (136) by 
d 2 Q 


. i 

■vrr + (ck ) 2 Q = — J. , for the low modes, 

dt z a a e i,a dt 


(138a) 


A5fc 


and 


original rage is 
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(ck ft ) 2 Q a ■ ~ J A a , for the high modes, (130b) 

We decompose the sums over the modes in (l3M into two sums, one over the low 
modes, the other over the high modes, so that 

1 * 1 + l . ( 139 ) 

a a, low a, high 

The deletion of the inertia term in (138b) ic< accounted for by a change in 

dQ a 

the action functional ( 13 * 0 , where we would have to omit the term I (“^r-) 2 . 

a, high at 

Accordingly, we approximate A . by 

A = ft dt{~ ^2— M* + U ' 4 . i ? JO + 

a,b .{ azi 2 dt M iJ dt * u i dt + c r. . J i,a W 'a dt 


appr 


a, high 
dQ 

l 


t*a 

* ~.,L '« ‘ V * ..law ‘ dt 


(ito) 


We express the expansion coefficients Q for the high modes by means of (138b), 

a 

We therefore write 

(ck )* e i,a dt (ck e , 1 ,a dt 

Cl O cl O * 

Here we must be prepared to use two different summation indices i and J , in 
order to avoid confusion later on. The two terms in (l4o) that contain the 

a 

of the high modes are then 

« - J i. A ^ * r ^ ( . L„ Ickl? J i,a J J ,a> ^ ' 


a, high 


a, high a 


and 


l ^ k a Q a^ 2 = 2e~ ^dt ^ £ ~ (ck ' T f tT i a J .1 a^dt * 

a, high a a 2e o dt a, high ,a 
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We see that they resemble each other. We combine these two terms and use 
the abbreviation 

<l4l) 




Then (l4o) becomes 


„ - /” ai. + U {aa_.v] + | l J 0% 

,b t 2 dt ij dt i dt c , i,a a dt 


a 
appr 




+ _fa i (— ) 2 - — i <k o ) 2 > 

^a,L 4t 5 a, low A ' 


( 11 * 2 ) 


The functions that are to be determined from the action principle (51) are the 
n functions q*(t) and the, now finite, set {Q (t)} for the low modes. 

8l 

When we reverse some of the procedures that led from (127) to (134), we 
can write (l42) also in the form 


A * { [i ^3— M ^2^- + U' _ v] + / ( [trJ • A + 

a,b ' dt ij dt i dt VJ r JJJW A i 0 v 

appr a 


+ til" n ' S >lov • <4 ' 5 lov> - T <’ * ° S lov> • • 


where 


I «A • 


low “ "a a 
a, low 

Because of the orthonormality relations, we may also replace the complete 
current density J in (l43) by 1qv , the low-mode part of the transverse 
portion of J, which is given by 


(143) 

(144) 


T,low 


l ^dt J i,a E a • 


(145) 


a, low 

However, the best policy is to work with equation (l42). 

We are now in the position to discuss the motivation for the approximation 
we have made. The expression (123) contains the integral 
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which is equal to the electrostatic energy associated with the atomic particles. 
(Remember that ♦ was the electrostatic potential, as shown by (12U) . But for 
point-like particles, this energy is infinitely large, because |v$| increases 
to infinity with decreasing distance from the particle so rapidly that the 
energy integral diverges. The other terms represented by triple integrals in 
(123) offer similar difficulties. For instance* the term 

III At< 7 « CAM? x cl) , 


which is the magnetic energy, is infinitely large for moving point-like particles, 
because the immediate vicinity of a particle contributes a diverging amount to 
the integral. Thus for point-like atomic particles, equation (127) is really 
without meaning. 

On the other hand, the approximation described by (lU2) appears to be 
meaningful, as neither diverging triple integrals nor infinite series occur in 
this equation. Of course, this appearance is deceptive, because the second 
terms in the equations 

e 

V* V* + ~ /// Arty • , (125) repeated 


and 


(1U1> repeated 

are infinitely large for point-like particles. But it is conceivable that the 
infinities of the second terms in (125) and (lUl) are compensated by infinities 
of the opposite sign in the first terms, so that the left-hand sides V and M,, 
turn out to be finite. It is generally believed, although there are no 
rational foundations for this belief, that this actually happens. It is 
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■believed furthermore that the square bracket in (l^2), which looks like an 
atomic Lagrangian, is in fact equal to the conventional atomic Lagrangian, 
if the critical frequency which formed the dividing point between low and 

high modes, is chosen low enough. This means that, for + 0* the 

expression ^ is equal to the conventional kinetic energy (in the 

sense of elementary particles) of the atomic particles, that V is their conven- 
tional potential energy (again in the sense of elementary particles) and that 
the quantities U^’ are 7<ero. 

We have introduced two approximations. 

(1) The approximation that was made when we deleted the inertial 
tern for the high modes, as in (138b). 

(2) The approximation that was made when we replaced the square 
bracket in (l42) by the conventional Lagrangian of the atom. 

The quality of the first approximation is improved, when we raise the critical 
frequency ft And the quality of the second approximation is improved when 

we lower the critical frequency ft j,^. Therefore, we have to make a compromise 
when we choose We believe that the quality of the overall approximation 

is good enough for our purposes, when we choose ft ^ to be of the order of ten 
times the dominant frequencies, e.g. ten times the frequency of visible light 
for the example of the laser. Unfortunately we have no evidence for this 
belief at the moment. But to adhere to this belief is the best we can do. 

At any rate, the approximations we have made and the beliefs we have 
adopted furnish us with the definite action functional (l42). Since our future 
work is based on it, we write it down again in a box. We use the statement 


that all 
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the ' are zero, and we indicate the variables on which the various quantities 


depend. 




*7 l + 

a, low * 

+ o?T I CflT Qjt )] 2 - ~ I r. 


2,c* if dt a 
a, low 


^a,lov [ V'a< t)12) ' 


(1U6) 


The expression in the double square bracket is the conventional Lagrangian ci‘ 
the atom, when the atom is isolated from* the radiation field. 


This action functional is of the standard type that one always meets in classical 
dynamics. It pertains to a system with (n + N) degrees of freedom. Here n is 
the number of the degrees of freedom of the atom, and N is the number of trans- 
verse modes whose resonance frequencies are less than the critical frequency 8 
The configurational coordinates are the n quantities q^, and the N expansion 

coefficients Q . The action functional depends on the functions Q (t) in a 
a a 

particularly simple way, characteristic of harmonic oscillators. There is no 
direct coupling among these oscillators, as shown by the absence of mixed terms 
of the type etc. However, the "radiation field" oscillators are coupled 

to the atom by the term on the second line of (1^6). Hence they are coupled 
to each other in an indirect way. 

The other form, analogous to (1^3), of the action functional is 
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We not^ that low (r,q k ) depends on the time t in an indirect way through 
the functions q (t). 

There is another way of looking at (lL6). We retain the first and second 
lines as they are. But in the third line we delete the subscript ’'low" . 

Then equation (138a) holds for all transverse modes, not Just the low ones. 

But the excitation term, i.e. the right-hand side of (138a) vanishes for the 
high modes. Thus if we use the initial conditions 

Q = 0»7rQ = 0 at t = 0 for the high modes, (1^8) 

we find that the high modes are, never excited, so that the are auto- 

matically zero at all times. The preceding remarks may be interpreted in the 
following way. Since, in the spirit of our approximation, the interaction term 
(i.e. the second line of (l46)) is a sum that extends only over the low modes, 
the real atom has been converted by our approximation into an object that is 
completely inert to radiation whose frequencies exceed the critical frequency 
^crit* ^is wa y then our approximate treatment differs from reality. But 
this deviation is not harmful in the applications we have in mind. In these 
applications we are interested in radiation whose frequencies lie in or near 
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those of visible light. 

We close this section with a qualitative argument which is designed to 
show that the approximation in (138b), where we deleted the term Q & , is 
good only in the non- relativistic domain, i.e. when the speeds of the atomic 
particles satisfy the condition (speed) 2 << c 2 . Let us consider Just one 
atomic particle of charge e. We assume that it travels with constant velocity 
w along the x-axis of some cartesian coordinate system. Its position r is then 
given by r = wtl, where i is the unit vector in the x-dlrection. Here we need 
only one configurational coordinate q^, which we choose to be equal to x. From 
the expansion (132) we determine J. by integration. We have 

X y B* 

J l,a * U^ 3 1 • \ 

Since J 1 is a distribution that is concentrated in Just one point, the integra- 
tion is easy to carry out. The result is 

J, = ei • E (wtl) 

X |8> & f 

As the time t progresses, the value of J. will change appreciably when the 

1 }8< 

1 

particle has travelled a distance, which is of the order of “ . The elapsed 
! 

time is then — r— . The inverse W = wk is the dominant frequency of . Now 

WK & 0. 1 jS. 

Or 

we calculate Q from (138b). It is proportional to J. . Thus the dominant 
a l ,a 

frequency of Q is also wk . Now we ask: Were we permitted to delete the term 
a a 

• 2 j2 

-rrr Q in (138b) ? Since rrrr Q is of the order of (frequency) 2 Q , our question 
dx* a crt a a 

becomes: Is(wk ) 2 |q I much smaller than (ck ) 2 | Q 1 ? The answer is yes, if and 

8. & 8 El 

only if (— ) 2 << 1. And this is the condition for the non-relativistic velocity 
c 

range. 

We have seen that the approximation in (138b) is of poor quality, if the 
particle velocities lie in the relativistic range. The converse of this 
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statement is not necessarily true. Non-relativistic velocities do not neces- 
sarily guarantee that the approximation is of good quality. After all, in the 
foregoing estimate, we have considered only a particle that moves with constant 
velocity. Nevertheless, we stick to our approximate treatment, because it is 
presumably the best that can be done without the use of advanced quantum electro- 
dynamics. Furthermore, it yields a definite action functional, which permits 
further analysis. 

In the next section, we are going to make the transition to quantum mechanics. 
We shall apply quantum mechanics only to the atom, not to the radiation field. 
However, the procedure that we are going to use can be easily adapted also to 
the case in which quantum mechanics is used for the radiation field as well. 

8. The Transition to Quantum Mechanics . 

The usual method for making the transition to quantum mechanics (or !, quan- 
tizing") is to define the generalized momenta and the Hamiltonian function, and 
then to proceed to the Schr<3dinger equation. But this means that additional 
concepts must be used. Another drawback of this method is that it can tell us 
only how the atom responds to the radiation. It will not tell us how, in reverse, 
the atom influences the radiation field. However, both of these questions are 
answered, when the action principle is used. Therefore, we stay with the action 
principle. All we have to do, in order to make the transition to quantum mechanics, 
is to modify the action functional (ll*6) in the appropriate manner. 

In order to simplify the presentation of this approach, we shorten the 
notation in (lU6). We use the abbreviations 

« i i(t) = vl(t) 

(This one we have used before), and 


r 


( 11 * 9 ) 
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c l (t) - U (150) 

a, low * 

(This U. is without a prime. The primed quantities were all zero, as we had 
mentioned in the paragraph that followed equations (125), (l4l) (repeated) of 
the previous section.) We also, delete the variables on which the quantities 
of (146) depend, i.e. we delete the variables in the round parentheses. Then, 
if we put the term U^v* inside the double square bracket, the shortened version 
of (l46) appears as 


a,b 

appr 


/ b dUCtl A^v* + U,v J - V]] + 


t=a 


i 


* I - "I I tk a 9 a ]2) 

a, low a, low 


(151) 


In order to quantize the atomic subsystem, we introduce a time-dependent 
complex scalar field ']>(q m ,t) defined throughout the configuration space of the 
atom. And we also introduce a time- dependent complex vector field in this 
space , with components X^(q. m »t). These components are referred to axes that, 
locally, conform to the coordinate lines of the configurational coordinates q^ . 
Then we replace the double square bracket in (151) by an integral over the 
entire configuration space of the atom, i.e. we make the replacement 

[ [| v J M Jk v k + UjV 1 - V]] -*■ 

/ Aa[[~ ((4)* II* - 4» ||- + - x^ ~rj ^ )) + 


ecs 


?q 


sq 1 - 


+ ((| x* J M jk x k + |(x% + xV)Uj - 4»%V))]] . 


(152) 


where the stars indicate the complex conjugate quantities. The symbol "h denotes 
Planck's constant divided by 2ir. The "i" in front of stands for vQT. The 
symbol A o denotes the volume element in the configuration space. In order to 
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explain the meaning of Ac, we consider a volume element that is shaped like a 
small parallelepiped whose edges conform to the coordinate lines of the configura- 
tional coordinates q 1 ^. The extents of the edges are A q l , Aq 2 , ... Aq 11 . The 
volume element Ac is proportional to the product of these extents, so that we 
can write 

Aa = MAq 1 Aq 2 ... Aq 11 . ( 153 ) 

The proportionality factor M is equal to the positive square root of the deter- 
minate of an n-hy-n matrix whose matrix elements are the quantities in (152). 
Thus, in simplified notation, 

M = (DetM Jk )^ . (15^) 

The letters ecs under the integration sign stand for "entire configuration space." 
They will he omitted henceforth. 

Some insight into the physical meaning of the field quantities ip and x 1 ^ is 
gained, when they are related to probabilistic statements. In quantum mechanics, 
the instantaneous configuration of the atom is no longer given by a definite 
point {q*^} in the. configuration space. Instead, we deal with a probability dis- 
tribution, whose density s is given by 

s = ip*ip . (155) 

This means that the probability A(prob) of finding the configuration to be in 
some volume element Ao is given by 

A(prob) = sAo = ij> t^Ao * ip MIAq^q 2 ... Aq 11 . (156) 

(The last member of this string pertains to the case in which the volume 
element is shaped like a parallelepiped of the type we described earlier.) 

Since the probability of finding the configuration to lie somewhere in the 
configuration space is unity, the field quantity ip should satisfy the 


normalization condition 
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Associated with the probability density s is the probability current density 
(in the configuration space) whose components are given by 

+ xV) • (158) 

The probability density and the probability current density are linked by the 


continuity equation 


*w-J W k3) ’° 


(159) 


as we shall see later. 


When we make the replacement (152) in (151) we obtain the quantum mechanical 
action functional . 

M a,t 4< i n H -♦ * X*' 1 fy’ 1 


appr » quant 


+ ((%X* J M Jk X k + + xV)Uj - *\v) ) ] ] 

e_ dQ o e ? 

*2& £ [ "at ] l I fc a\ ] } • 

a, low a, low 


(160) 


The functions that are to be determined from the action principle are i|)(q ,t), 
the n functions x^(q*,t), and the N function Q (t). We repeat the definitions 
of n and N. The number n is equal to the number of the atomic degrees of 
freedom, and N is the number of transverse modes whose resonance frequencies 
lie below the critical frequency n crit * We wight also say that N is the number 
of electromagnetic degrees of freedom that we retained after the approximations 
of the preceding section. When we apply the action principle (for a free-running 
system), we may regard the variations 6 i|j, 6i|» , , $X 0 as independent variations. 


The Uj in (l6o) are given by (150). 
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even though ^ is tied to ^ and X is tied to x by the condition that the 
starred quantities are the complex conjugates of the unstarred ones. What 

i 

enables us to do this is that we will have to deal with quantities of the type 

# # 

F6\f> + F , which have to be equated to zero. Choosing 6 'I' to be purely real 

# ft 

yields F + F * 0. And choosing 6^ to be purely imaginary yields F - P * 0. 

* 

These two equations imply that F * 0 and F * 0. These are Just the results 

ft 

that we would have obtained, had we regarded 5\p and 6ip as independent varia- 
tions. 


Had we desired to quantize also the electromagnetic field, then we could 
have extended the procedure that led from (151) to (l60). In this case we would have 
used, instead of an integral over the n-dimensional configuration space of the 
atom, an integral over the (n + N) -dimensional configuration space of the total 
system. We need not elaborate any further upon this theme, since we decided 
at the outset to quantize only the atom and to treat the radiation field 
classically. 


Now let us apply the action principle 

SA = E.T.O. 

EL j D 


(5l) repeated 


for a free-running system. We use the action functional (l60). If we vary 

ift 

only the x » we obtain 

6A = / b dt{/Aq l Aq 2 ... Aq 11 M fix + ^Mj V X k + 'll) } 


a ’ b t-a 


2 9q J 


Equation (51) tells us that the factor of 6x J must vanish. Thus 


(161) 


We introduce the matrix with elements M J , which is the inverse to the matrix 
with the elements M^j. (Both are symmetric.) The condition that these two 


Of POO* 
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matrices are the inverses of each other is reflected by the following equations 
for the matrix elements 


V 


Jk 


* i M 1J M <V 
i «Jk 


• ( 162 ) 
(The summation convention remains in effect), where the 6^ are the Kronecker 
deltas. Then, if we premultiply (l6l) by M^, we obtain 


M 


* J (- 1* - M) 


Similarly, 


3q' 

* 




( 163 ) 


X* 1 * M* J (+ & u 

# , ^ « 

Now let us vary only v» . In 6 A . there appear the terms - rr- and 

a, D dt 

J 3 * 

X — - , which require some manipulation. We have 

jjqJ 

- ^ ~ 5\j/ * - ^ (WW*) + M H' * 

The first term, being a full time-derivative contributes only to the end terms 
E.T.O. The last term vanishes, since we assumed that the matrix elements M,. 

jK 

do not depend on the time t. Thus * (Det M does not depend on t, so that 

JK 

M 

3t 


* 0 . Similarly 


- Mx J « - -4 + a/ -4 <Mx J ) 


The first term is a full configurational derivative. When we integrate it over 

1 * 

the entire configuration space, the result is zero, because MX' vanishes in 

# 

the far reaches of the configuration space. (We could even require that 6^ 
vanishes except in a finite domain. ) 


Combined with the results of these manipulations, equation (l6o) yields 


6 *a,b * " ^2 **• A{ i n M^ # > 

t»a 
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* / b 4t{/iq'40.2 ... Aq n H(«»*)[[Hf §f + H ^ (M X J 0 * 


t*a 


2 M 3^ 


+ |«jX J - V*]]} . 


The action principle (51) tells us .that the expression in the double square 
bracket must vanish. Thus 

~iv J + v <" 

Similarly 


2 *V 

On combining (163) with (165), we obtain 


(165) 


(166) 




(167) 


or, in somewhat abbreviated notation, 


** f * w . ( - - u j )m3 * ( - ^ ^ - v* + 


(168) 


Equation (168) is the well-known Schrddinger equation for time- independent 
"kinetic coefficients" M.. , mJ\ m. We will refer to the pair of equations 

jK 

(l6l) and (165), from which (168) could be derived, as the "pair of SchrBdinger 
equations". The influence of the radiation field is contained in the quantities 
U , which - according to equation (150) - depend on the electromagnetic vari- 

V 

ableB Q (t). Similarly, when we consider the variations of tji, we obtain 
a 

- itf ||- = -jj (+ Mi'-— - U )MM Jk (+ + Vtf>\ (169) 

^ 3q 3q 

This equation is, of course, the complex conjugate of (l68). 


At this point we are ready to check the continuity equation (159). We have 

_3 

3t 


TT* 8 + | 


ij (Mk j ) - (<m») + [M i(/ x J + 'lor 5 )) 

# w + 1 * J ^ + »* in 

where c.c. stands for conjugate' complex. With (165) we obtain 





where, in the last step, ye have used the complex conjugate of (l6l). But 

» * 

this expression is purely imaginary, since the factors of 1 are purely real. 
Thus, when we add the complex conjugate "c.c.", we get zero. And this com- 


pletes the check. 


So far we have exploited the action principle only in part, in that we 

* 1 *1 

examined what happens when the quantum-mechanical quantities \ p 9 , x » X 

are varied. Now we subject the electromagnetic quantities Q to variations 
6Q . Here we must not forget that, according to equation (150), the Q a are 
contained in the U^. We therefore write down the action functional (160) 
again, but this time with the Q a exhibited wherever they occur. All the 
other terms are indicated by a sequence of dots. Then, with the use of (150 ), 
we obtain 


A = / b at{ l a [/iqV ... i4 n M §<x% + J, a U k >] + 

a * D t=a a, low 

e dQ ? e. 2 

+ ac? I [ “dt J ““2 l [k a Q a ] } * (1T0) 


a, low 


a, low 
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<!Q a 2 


In the variation $A ft b there appear the terms » vhich require some mani* 

pulation. We have 


r d l,2 

dt 


dO. dQ 
dQ 


dQ . 

2 dt dt 
d 2 Q 


’ 2 dt [ "3t 5, « ] * 2 "at* 4 V 

The first term, being a full time-derivative, contributes only to the end terms 
E.T.O. We then obtain 

g dQ 

4A a ,b * • • • + ^ i [_ dt 6 Q a ] 

a > low i fa 

e d 2 Q 

♦ / b dt l 5Q o { - -££7 - «JkJ*Q.+ 


o a 


t«a a, low 

+ / Aq^q 2 ... Aq^|<x% + X’V)^ Jj #a (q k ) } • 

» 

The action principle tells us that the expression in the curled brackets must 

vanish . Thus 
g d 2 Q 

c? “at* + e o (k a )2 °'a = c ^ q ' 4q2 — ^ + xV)J J>a 

These N equations (one for each of the low transverse modes #a) describe how the 
atom influences the electromagnetic field. Each of these equations is constructed 
like the equation of motion for a driven harmonic oscillator. With the aid of 
(153) and ( 158 ), we can write (lTl) in more compact form, namely as 


(171) 



after we have divided it by g q . This equation is the main result of this section. 
It is the quantum-mechanical analog of the classical equation (138a). 

In order to explain the meaning of (172), we write it down again in con- 
junction with (138a), also divided by e . Only this time we indicate the vari- 

d i 

ables on which the various quantities depend. In (138a) we replace by its 
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0fr po 0r 

Classically! 

o 

Quantim-mechanically : 

Tdiry + M ’ T7 Ao K J (q k ,tu (q k ) . <m> 

o u w 

Of course, the words "Classically" and "Quantum-raechanically" refer only to the 
treatment of the atomic subsystem. The electromagnetic subsystem is always 
treated classically in our presentation. 

In order to make the transition from (173) to (17^) more evident, we rewrite 
the right-hand side (R.H.S.) of ( 173 ) in a cumbersome, but instructive, way. We 
could regard the precisely defined configuration (q^Ct)} as the reflection of a 
highly concentrated probability distribution in the atomic configuration space. 
Let s(q^,t) be the density of this distribution. It is sharply peaked near the 

j. f 

point {q (t)} and zero elsewhere. When we integrate s over the entire configura- 
tion space, the result must be unity. Of course only the immediate vicinity of 
the point q^(t) contributes to the integral. We can therefore write the R.H.S. 
of (173) in the form 

R.H.S. = -i— {/Aa s)v^ (t)j Cqk'tt ).) , 
e o c -3 

Here s in the integral is evaluated at the generic or "running" configuration of 
the integration process, whereas M. is evaluated at the precisely defined 

J 

k k 

configuration { q. ( t ) } . But, because s is so sharply peaked near {q (t)} we 

may pull J, into the integrand and evaluate it at the 1 unning configuration, 

J ,a 

without changing the value of the R.H.S. We may pull in the v^(t) as well, 
because this quantity figures as a constant in the integration process. Thus 


p *ee is 

Walit? 
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But sv^ is the component of the probability current density IC^ in the atomic 
configuration space. Like s, the as well are sharply peaked near the point 
{q^(t)}. We then obtain 

R.H.S. * /doK^J, a . (175) 

e 0 c j » a 

And now the R.H.S. of (173) looks like the R.H.S. of (17^). This shows the 
intimate connection between the classical and quantum-mechanical equations (173) 
and (17^). 

We may recapitulate the considerations of the preceding paragraph in the 
following way. The R.H.S. of the classical equation (173) is deterministic in 
form and in substance, since it refers to a precisely defined configuration 
(q k (t)} and velocity with components (v^(t)}. But we may contrive to write it 
in a probabilistic form, as in (175)* In substance, it is still deterministic, 
because the probability current density components are such highly concentrated 
distributions. When we go over to the quantum-mechanical equation (17^) » the 
same probabilistic form is retained. Only now the R.H.S. is probabilistic not. 
only in form but also in substance, because the quantum-mechanical are no 
longer highly concentrated distributions. We may therefore regard (17^) as the 
direct probabilistic analog of (173). 

Now that we have elucidated the meaning of the quantum-mechanical equation 

(17^), we are permitted to make the transition from (173) to (17^) in a rather 

mechanical manner, i.e. by a purely typographical procedure. All we have to do 

is to replace the symbol v^ in (173) by the symbol /AoK*^ and to pull in the 

symbol J. under the integral. We do this, while the independent variables 
j » a 

are deleted. Afterwards we restore the appropriate independent variables. Thus: 
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(q*(t)) + v J J Jia * /AoK^Jj^ * /4oK J (q k ,t)J J((l (<l k ) . (176) 

We ought to make 3ure that the ultimate independent variable i3 the same 

on both sides of (173) and also of (.174). The left-hand sides depend only on 

the time t. And in fact, the right-hand sides have only t as the ultimate inde- 

lc 

pendent variable. In (173) thi3 is assured, because the q depend on t. And 

k. 

in (174) this is assured, because the q get "integrated out". 

The equations (173) and (17 1 *) tell us how the atom influences the electro- 


magnetic field. In the form we have written them, they describe how the atom 

affects the individual mode excursions Q Q (t). We may also write down equivalent 

equation for the vector potential field cA(r,t). To this end, we multiply each 

of the N equations (173) and (l7 l 0 (one for each mode H a) by the corresponding 

transverse mode pattern E (r). Then we sum over the low modes #a. We use 

a 

equation (9l) (See Box of Section 6), namely 

(k ) 2 E = V x (V x E ) , (91) repeated 

a q a 


and equation (l44), namely 

cA loy (r,t) = l Q a (t)E a (r) * (144) repeated 

a, low 

Then we obtain 
Classic ally: 


w «*!»«.*) - 5 11 (5 «<=A l0W (5.tn 

Quantum mechanically: 

+ ’ * (5 x cS lov (? * t>) 



= 7^/Aa>C>( q k ,t)J 3ia (q k )E a (?). 
o 


(177) 


( 178 ) 


Next we introduce the abbreviation 


I 7 

a, low 


J i a 


(q k )E a (r) 


,T J,T,low (qk ’ 3 


! 


(179) 
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which is suggested by equation (l45). Then equations ( 1 T 7 ) and (178) can be 
written in more compact form: 


Classically: 

TOFF cZ Iov (?,t) + * * [? * 0? lov (? * t)1 " r? • 


(180) 


Quantum mechanically: 

TiSF oS 1 ow (f ' t) - 9 X TCV x cA l0W (?,t)] = ^ /AaK J ( q k > t)j J>T(low (, k .?) . (181) 

We may also reverse the procedure that led from (173) or (17*0 to (180) or 
(l8l). All we have to do is multiply (l8o) or (l8l) by the transverse mode 
pattern E (r), and then to integrate over the cavity. As far as the details 

6L 

of the calculation are converned, one uses the series expansions (l44) and 
(179), and also the orthonormality relations ( 98 ). (See Box 05 of Section 6 .) 

One also uses term-by-term integration and differentiation of the expansion 
series. This is certainly permissible, since all the series are finite. 

We Just went through a procedure that carried us from equations ( 173 ) and 
(174), the equations of motion for the mode excursions, to equations ( 180 ) and 
(l 8 l), the equations of motion for the vector potential field cA(r,t) as a whole. 
We could have gone through this procedure also in the action functional. In the 
case of classical dynamics for the atom, we have done it already, as we went 
from (l46) to (l47). When the atom is treated quantum-mechanically , the field- 
analog of ( 160 ) would have been 
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l Rfb - /* « .{/Aq*Aq* ... A^AfCC^ </ H - * ^ ♦ X* J -2f - X* 1 **J) * 

appr, quant. W 3q 3q ‘ 

♦ ((| x #J M Jk x k + §(>c4 + x¥)Uj - Av))3) + 

+ "fW 4Tt »St cI iov (?tt > ] ‘ [ alt c W*’ t)J - 
- ///at[? X cl l 0 v (?,t)] • [7 » cA low (?,t)]> ( 182 ) 


For Uj(q ,t) in this expression, we would prefer, instead of (150), the equivalent 


formula 


UjUV) ///at; J J(T> low (a k ,?) • ci low (?,t) . 


(183) 


All the equations of motion for a free-running system follow from the action 


principle 


E.T.O. 


(51) repeated 


The transition from (173) and (174) to (l8o) and (181) in the equations of 
motion for the electromagnetic field does not produce any new results. It is 
merely a change in the mathematical idiom. However, it affords us some insight 
into the physical significance of the excitation terms, namely the right hand 
sides. In the classical equation (100 ) the excitation term (except for the 


factor — -~) is the low-modes portion of the transverse current density 

C C 

T 4- \ - . J / + VT / h / . \ ~ \ 




(18k) 


In the quantum-mechanical equation (l0l), this quantity ^ ow (r,t) replaced 
by its quantum-mechanical expectation value. 

So far we have considered only one atom in the cavity. The next question 
that arises naturally is: How does one proceed when several or many atoms are 
present? We take up this question in the next section. We omit the classical 
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treatment of the atoms and proceed right away to the more interesting case 

in which the atoms are treated quantum-mechanically. 

9. The Quantum-Mechanical Problem of Many Atoms. 

In order to explain the principle of the method that we shall employ, we 

first discuss the case of Just two atoms. We assume that the two atoms are 

far enough apart so that they do not influence each other directly. But they 

are permitted to interact with each other in an indirect way, namely through 

the low-modes portion of the vector potential field cA(r,t). 

The first atom has n' degrees of freedom. Its configuration is described 

k * * 

hy the n’ coordinates q , (k* = 1,2, «.. n*). The second atom has n" degrees 

k" 

of freedom. Its configuration is described by the n" coordinates q , (k" * 1, 


2, ... n” ). The correct procedure would be to treat the two atoms as one super- 
atom with n = n’ + n M degrees of freedom. One would then start with the action 
functional, e.g. of the form (182). But this procedure is much too unwieldy, 
especially so, when many more than two atoms are involved. Therefore, we shall 
look for some assumptions that will simplify the treatment. In order to spell 
out these assumptions, we have to devise a consistent scheme of notation. 
Single-primed symbols and symbols with single-primed indices refer to the first 
atom, while double-primed symbols and symbols with double-primed indices refer 
to the second atom. Unprimed symbols and symbols with unprimed indices refer to 
the super-atom. Single-primed, double-primed, and unprimed indices run respec- 
tively from 1 to n' , from 1 to n M , and from 1 to n = n' + n 1 '. We relabel the 
index set 1,2, ... n' , n* + 1, n' + 2, ... n 1 + n" as 

1*, 2’, ... n', 1", 2" ... n" . 
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One consequence of this relabeling is that the product of differentials Aq l Aq 2 . . . 
Aq n that occurs in (182) can be rewritten as (Aq 1 Aq 2 ... Aq n HAq 1 Aq 2 ... Aq 11 ) 
Next we examine the matrix of the kinetic coefficients with elements Mj^. Because 
of the relabeling, this matrix becomes partitioned into four boxes, namely 

Matrix (Mj^) 


Matrix(M ) = 


M ,1>k' 

M ,1’k" 

Vk’ 

M J"k" 

each 

other 

the off-dia 

Vv 

0 

0 

Vk" 


(185) 


A consequence of (185) is that the appropriate determinants satisfy the relation 


Det M 


Jk 


( Det Mj | ^ ? ) ( Det Mj ) 

le 


Then the quantity M - (Det M.. which occurs in (182) can be written as 

jK 

M = M» M" 


(186) 


(187) 


in rather obvious notation. We note also that M * and the matrix elements M^,^, 

t ' a 

depend only on the q , while M" and the matrix elements depend only on 

the q . The nature of these dependences is the same that we find in the single 


atoms. Again, because of the appearance of the zeros in (185), we have 




k _ *J 


X 'H jlX = X 


M 




j, k ,x m 3 , v ,x 


k” 


(188) 


k’ fL* l n 

But the x depend on both sets of coordinates q , q and the time t. And 

k" 

so do the x • 

To continue with our examination of the terms that occur in (182), we note 

that 

(x*^ + X J ^ U j = (x* 3 '»f» + X^'^)Uj, + (X* 3 + ^ V )Uj" (189) 

l' A" 

The quantity in (189) depends on both sets of coordinates q , q and the 


Jw * 

time t. However, the Uj, depend only on the set q , and the U^„ depend only 


W* n l. K 
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on the set of q . We see this when we examine equation (183) » which defines 
the Uj. In this equation, there occur the quantities Jj T lQW (q k ,r), w ^ich are 
of purely kinematical origin, in that the low-modes portion of the current 
density in our home-space is given hy 

W ftt> “ vJ ’ (t)7 J,T,low (<lk(t) * ;;) 

After the relabeling that we introduced earlier this equation becomes 


+ v^ (t)jj„^ low (q k (t),q k (t),r) (190a) 

Since the contribution to J_ . made by the first atom (the first term in 

x 9 low 

(190a) cannot depend on the configuration of the second atom, we conclude that 

_ jji - 

J y T can depend only on the set q . Similarly T low can depend only 

on the set q k . Equation (l90a)thus becomes 

^.lov^*’ ■ + ^’V.T.loX’ <*>>*> ; (190b) 

The equations that define the Uj , and „ (see (183)) then become 

U J ,(q k ’,t) =■ i///4Tj' j ,, T>lOT ( q k, ,f) • cA lQu (?,t) , (191a) 

and 

>b) = q* ,T,low^ * ^iow^’^ * (191b/ 

k' 

They show that the U^ f can depend only on the set q and that the Uj„ can 

k" 

depend only on the set q , as we already indicated by the way we wrote the 
left-hand sides. 

Furthermore, since the two atoms do not interact directly, the potential 
V(q k ) that occurs in (182) will be of the form 

V(q k ) * V'(q k ') + V"(q*") 


9 


(192) 
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where the two functions V’ and V !! are the potentials for the single atoms. 

Now we are ready to incorporate our previous considerations into the 
action functional. We shorten the notation somewhat hy the use of the volume 
element in configuration space. On combining (153) and (187) ♦ we have 
A<x • Aq J Aq 2 ... Aq n Af * (Aq^Aq 2 * ... Aq 11 M’ HAq^'Aq 2 " ... Aq 11 M" ) * Aa'Acr" 
Also we use two integration signs in front of Ao’Ao" in order to remind our- 
selves that we have to integrate over both configuration spaces. With these 
changes of notation equation (l82) becomes 


a,b 

appr, quant * a # 


/ b «{//*• W[[J$ -a. ♦ X* J " -4r - 


3q 


3q° 


Y y is Y r w > . //i *y M Y k» . i *r M Y k" 

-x ■ X m ^W) + 1(5 X M j»k* X + ^ x M 4 „ V „X + 


3q' 


9q c 




1 ' *1 » i» #. 1 / *1” 1" #, 

+ j (x J * + x J * )u. f + i(x J * * x J * )U,„ 


where we 
keep the 
And 
We shall 
are of a 


- - *V'»]] + 

+ "T W AT ^3ct c *lcv ] 


^ 3ct cS lowl 


- T ///»tW x cX lou ] • X =A lov ]) , (X93) 

have omitted the indication of the independent variables in order to 
notation as short as possible. 

now we come to the statement of the assumptions we are going to make, 
assume that the (n f + n" + l) functions (q k ,t), X^ (q k ,t), ^ ( q.^ , t ) 
rather restricted form, namely that 


t|)(q k ,t) = i|/’(q k ,t)^"(q k ,t)% (l9Ua) 

(q k ,t ) ='x^ (q k ,t)i|» f '(q k ,t) , (l9^b) 

X J (q\t) = "x J (q k ,t)i|»'(q k ,t) , (X?Uc) 
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(Similarly for the complex conjugates, the starred quantities.) Furthermore, 
the functions ^ and i|/" shall fulfill the normalization conditions 

* 1, jVy‘V “ 1 (195) 

(Actually we need to postulate these normalization conditions only for one 
particular time, say t = 0. Then they will be satisfied at all times, since 
the Schrddinger equations for 0’ « x° and those for i|>", x Jf? turn out to be 
satisfied. And this, in turn, guarantees the validity of the continuity 
equations of the type (159), once for the single-primed quantities and then 
again for the double-primed quantities. And when we integrate the continuity 
equations over the appropriate configuration spaces, we see that the normaliza- 
tion integrals of (195) do not depend on the time t.) 

As we shall see, the assumptions (194a) - (195) will yield a workable 
action functional, from which - by way of the action principle - equations of 
motion can be derived that are well defined. Of course, because of the restricted 
nature of the quantum-mechanical functions as embodied by (194a) - (195) » only 
a restricted class of system histories can be obtained from these equations of 
motion. The histories that are rejected by our assumptions are, though possible 
ones, not of technical importance at the moment. • 

When we use the equations (194a) - (194c), but not yet the equations (195), 
in the action functional (193), we obtain, after a reordering of the terms, the 
following expressions 


a,b 
appr, quant 


• / b d-b< / /Aa 1 f ’^T " ♦* + *X ^ 9 ^V _ * x -3 * _ii_. ) + 

t*a 9 a <3 ' 3qJ 

A8 ° 

+ ((| ’X #J? Mj, k , , X k + |<'X* J V + ’X J V * )Ujt - 'P •♦•V , )>]] ♦ 


i ft 

jl H 

+ //Ao'4o"**'*'[[i|(/" - r + "X 3 -^F - V 2±yr) + 

X" ^ ^ 

+ ( (§■ "x J Mj., k i."x l! " + f("x*' , V + "X J %*”)U,„ - *‘"tv'))]] + 

low ] • llfeew 

- •% M® t? * c5. 1 • [?> o5, 1} . 




2 j j t. ' ~ '- ri i ow J ' Lr " ''"low J J • (196) 

In the first configuration-space integral of ( 196 ) we first integrate over the 
double-primed space. Because of one of the normalization conditions (195), this 
portion of the integration yieldsunity. Thus there is only the integral over the 
single-primed space left. We proceed in a similar way in the second integral, 
where we first integrate over the single-primed space. Thus the expression ( 196 ) 
simplifies into 


a,b 

appr, quant 


■ / b « t /to.[tV' *Sr- Sfc* "x* J ’ -*$■- •/ *jf) * 

t=a Sq 3q 


+ ((| ’x^’mj^.’x^ + |(’x* J V + ’x J V')u jf - *%'V'))]3 + 
+ /Aa" [[%♦*" ^ ’ *" ^ + n X* r ^Ttr - V' &W) + 


9 9 


9 9 


1 #4>* V' ( 1 *4 It 4 It #tf #!! ft 

+ ((5 "X J M J(V ,"x k + ^("x J *" + ”X J ¥ )U„, - ^ * V"))]] + 


+ -f CA,1 • [r^r cL 3 - 


3ct low 
^low " 


- -| ///At [v x cA^J . [7 x cA w ]} . 


Set low ' 

low 


(19T) 


We see that the quantum-mechanical portion of the action functional has 

been decomposed into two completely separated portions, one for each of the 

two atoms. A consequence of this observation is that the quantum-mechanical 
1 ’ 

functions 'x of the first atom satisfy their own pair of SchrSdinger 
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equations, namely 


ii .'*Jr - 

3q J 


M 


J’k’ 


k 

X 


♦V 


0 
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(19B) 


and 



in analogy to (l6l) and (lf>5)» Equations (198) and (199) ensue from the action 

#4 t « 

principle in the usual manner, when we vary ’x and tp *. Similar statements 

pertain to the second atom. The foregoing remarks do not imply however that 

the two atoms are completely independent from each other. For they interact 

through the low-modes portion of the vector potential field cA, , which is 

low 

contained in the and „ according to(l91a) and(l9lb). The equation that 
governs cA^ qw is obtained from the action principle, when we vary °A^ 0W * Of- 
course, when we perform this variation, we must not forget that cA]_ ow is con- 
tained in the Uj , and Uj P . Then, with the abbreviations 

'K* 3 * * |('X #J V + ’X J V’) , (200) 

and 

*n n #4ti ii 4 » nil 

"K 3 = |(" X J Ip + V ") , (201) 

in analogy to (158), we obtain 


cAi ow ♦ V K (V K CA 10V ) = /AC V J* 3 , iT>low * 

+ T7 ! l °” (SOS) 

O 

This equation resembles (l8l). However, each atom contributes separately and 
independently to the exciting term, the right-hand side. 

We could have skipped the entire preceding discussion of this section. 
Instead we could have postulated (197) as a workable action functional. But 


32 

tj* 
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then we woul4 not know what kind of assumptions He behind this postulational 

* 

approach. But now that we have recognized the assumptions, namely the restricted 
product form for the quantum-mechanical functions, we can use the postulational 
approach for the case of more than two atoms. Let us tag the atoms with the 
index p. Then, when we write down the analog of(l97) we can afford to omit the 
index p for most of the symbols. Instead, we write it below the integration sign 
of the configuration-space integral. In this way, we imply that the p should be 
appended to each symbol behind the integration sign. This shortened notation 
should not cause any confusion. Then the generalization of (197) is: 


A u 

a,b 

appr, quant 


/ b l /Ao[[^~ (ip ||- + X #J “jr - X^ ^t) 


t*a atoms#p p 


9q 


3<1 


+ ((| x" J M JIt X k + §<X% + xV)U, - /*»>»]] 
+ °X lov l • t'KT 


act low J 


- -f ///dr [7 x cA 1rt J • [? x c 5 1<w ]} . 


low 


low' 


(203) 


where 

i 

U * — ///at J 4 _ . • c)L . . (204) 

p J c p J ,T,low low 

With each atom #p there is associated its own pair of Schrbdinger equations, 
which are built like (198) and (199)* And the generalization of (202) is 
*2 - - - - - 1 


cA, 


+ V X (v X CA. ) = 


l ! 4o K J J, 


® " low ' ■ V atoms #p p ~ “ “J- T * low 

This equation shows that each atom contributes separately to the exciting 


(205) 


term. Each of these contributions is the quantum-mechanical expectation value 

of the current density for the atom in question, except for the factor . 

e o° 
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The equations of motion for the whole system consist of the pairs of 
Schrbdinger equations, one for each atom (An example of such a pair is given 
by (198) » (199)), and equation (205) for the vector potential. These equations 

1 

are intertwined, because the SchrBdinger equations contain the vector potential 
(by way of the U.), and because equation (205) contains the quantum-mechanical 
functions (by way of the K^). It is impossible in practice to solve such an 
intertwined system exactly. Therefore, one has to use approximate methods for 
the solution. The favorite method is the perturbation calculus, which we shall 
discuss in the next section. 

10. The Perturbation Calculus . 

We start with the action functional (203), in which we regard the term with 
the U. as a perturbation. It is this term which causes the equations of motion 

J 

to become intertwined. As is the custom, we multiply this term by an "expansion 
parameter" X, which in the end we set equal to unity. We also expand the quantum- 
mechanical functions X‘^ for each atom #p and the vector potential cJ 1qw as 
power series in X. Accordingly we write 

(^)p * (ip Q + Xih + X 2 i|> 2 + ...) p , ( 206 ) 

(X j ) p “ (Xo + Xxf + X 2 Xg + ... ) p , (207) 

and similarly for the starred quantities. In the expansion for c A low , we delete 
the subscript "low" on the right-hand side, in order to keep the notation as 
concise as possible. Thus 

C ^low = C % + + + .«. (208) 

The terms with the index 0 constitute the zero-order approximation, those with 

index 1 constitute the first-order correction, and so on. 
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We insert the expansions (20 6) - (208) into the action functional (203) and 
multiply the terms with the (which wo write out according to (204J) by X. 

In doing so, we multiply the inserted series term-by-term, but break off after 
the terms in X H (for the N x -order perturbation calculus). We indicate this 
breaking-off by the prefix [8(X N )1. In our applications, we will be satisfied 
with the second-order perturbation calculus. Thus we break off after the terms 
in X 2 . The action functional (203) is then replaced by an approximation, which 
we denote by A . ("pert" for perturbation). We refrain from writing out the 

fty D 

inserted seriei e Snd the results of the multiplications in detail, and simply write 

A., b - [«**>!/ *< l 1 A«lt4 (/ f - * 4 + x* J 5 - x J 

p* r £ t*a atoms #p p 3a 

+ ((| X #J M X k - V* + \ (x% + xV£ ///Act 1 • ))]] 

+ T / / / At ["aSt" cI iow^ ' ^3ct c5 iow^ 

- J//4 t[V x 0 A low ] . [7 x cA lov l) (209) 

This is a workable action functional. 

When we apply the action principle, we can vary all the functions , 4*1 * 

IH 

if>2 ... • cA independently. We start by varying the x J for atom # p. Then, in 
the usual manner, we conclude that we must have 

tB(X 2 ) + X«X* J + X 2 «X* 3 + X + X 2 Itj.) 

3q 3q' 3q 

+ M. v (Xo + *X? + * 2 X,) + + ^l)~ ///Act J • (cA 0 + XcAj)]} = 0 

J ,T,low p 



( 210 ) 
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Here we have written out the series in detail. In the last term we could 

shorten the series, because of the factor X in front of it. Now we remember 

*1 *1 *1 

that 5 Xq t $Xj , Sx 2 can t>« chosen independently. We also remember that we 

9 * 1 

break off after the terms in XH Then, on equating to zero the factor of 6 X 2 » 


we obtain 


* 2[i<1 ^3 + V^p ■ 0 • 


( 211 ) 


The factor X 2 may be cancelled, of course. Next, we equate the factor of 6 X 1 
in (210) to zero. We obtain 

X[ih + M .X^L + X 2 [ ill + M ..Xi Ho ^ ///A t J • cA 0 ] * 0 . 

3q J P 3q J C j,T,low P 

Now, the first term vanishes, because of (211). Then, after we cancel the 
factor X 2 , we obtain 




[ii( -^4- + m .x^ + 'I'o r ///at 7 • cS 0 3 » 0 . 

3q J JK C t1,T,low P 

Finally we equate the factor of in (210) to zero. We obtain 

[Ml — T + M. k X k 0 L + ^ + M -Tf + * | ///At J • cA ] + 

3q J JK 0 P 3q J ° C J ,T,low P 

+ X 2 [i 6 + M.. X^ + 7 ///At J • cA + 

3q J JK 2 C J ,T,low 


( 212 ) 


+ 'I 1 ! 7 /// At 7 


0 c 
cA 0 ] = 0 . 


J ,T,low P 

Now, the first and second terms vanish, because of (211) and (212). Then, after 
we cancel the factor X 2 , we obtain 

+ M x * + *0 ^ J//A t 7 • cA + H ~ ///At 7 . cA ] = 0 

3q J 2 C j , T , low 1 C J ,T,low Up 

In the next step, we vary the for atoms # p in (209). In the normal 
manner, which involves some integrations by parts, we conclude that we must have 


[B(X 2 )]{(S^ + X6<h + X 2 6<f, 2 )[i*(2|j. + X % + x 2 %) 

ot 0*6 3t 

+ 4 iT~T (M(x 0 * x 4 + * 2 xi» - v(t p + + X 2 ^) 

3q 

* \ (Xo + x x{)^ ///At 7 (cA 0 + XcAi)]} = 0 . ({ 

J,T,low p 

We apply the same step-by-step procedure that we used before. That is, we 

equate to zero first the factor of 6 i|> 2 , then the factor of fit//*, finally the 

* 

factor of And at each step we use the results of the previous steps. 

Then we obtain 


(2lU) 


[itf + Ml _2_ (Mx J) - v* 0 ] p » o , (215) 

[!< 8|i + (Mxi ) - V*! * | Xq | ///ir j . cA 0 l = 0 , (216) 

3 <1 J ,T,low p 

% * 4aT 71 (M 4> - v *2 + I Xq I ///At -J • cA + 

3< 1 J ,T,low 

+ \ x 'l e / / / At . 7 * c5 oL= 0 . (217) 


J ,T,low 


J,T,low 


cAj + 


J ,T,low 


(217) 


In the last step we vary cA 1qv in (209). In the usual manner, which involves 
some integrations by part we conclude that we must have 

[ 8 U 2 )](«cA 0 + UcA, + > 2 «cA 2 ).{ l /Ao | [( X *J + A X JJ)(, + X , ) + 

atoms #p d 0 

+ (Xo + Xxi)(^5 + X X *)] 7 J 

J ,T,low 

gZ _ _ 

" e o Tact) 2 ^0 + XcAj + X 2 cA 2 ) - e Q V x (7 x (cA 0 + XcAi + X 2 cA 2 ))> » 0 . 

(218) 

Now we apply the same step-by-step procedure that we used twice before. It 


yields 


3 2 

' e o ° A ° ■ *o ? * (5 * cA 0> - 0 , 


(219) 
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3 ? 


op] f>. 

***%&« 


- e 0 TsHTT cS > - c o 5 x (? x c5 i> + 

♦ i /act i( X ; j * 0 + X V)| j 

atoms ip p o o e a^ 1ow 

■ c 0 Tsflr 5 ' c5 2 * E 0 5 * <? X c5 2) ♦ 


0 , 


(230) 


+ l /Aa |(xJ J <h + Xi J ^0 + Xo^l + Xi^o) 7 5 * o . 

atoms #p p J,T,low 

(221) 

We may rewrite these equations in a way that appeals to our physical 
intuition. According to equations (200) and (201), the probability current 
density components (K^) p for atom #p is given by 

(^)p * §(X% + xV) p • (222) 

On the right-hand side we use the expansions (206), (207). We use a similar 
expansion for the left-hand side, namely 

(K^p ■ (K J 0 + XKj + \ 2 K^ + ...) . (223) 

It will suffice for our purposes, if we break off the ensuing expressions after 
the term in X 1 . Then equation (222) becomes 

OcV XKj ) p = [B(X)]|[(xJ <j + + Wi) + (xl + *xJ)(¥o + X^)] p • 

Comparing the terms in X° and X 1 , we read out 

(Kj) p - |(xj^0 + xi*I) . (224) 

(kJ.) p = |(Xo J ^i + x[^o + xln + Xi’J'o ) • (225) 

Then, instead of (219) - (221), we can write after an obvious rearrangement 

9 2 

( 3 °t) 2 CA 0 + v x (v * cAfl) - 0 , 


(226) 
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♦ V» (V xcX,) -Jy I 

o atom 


c c 

o atoms /f'p p 


/AoK^ J , (227) 

J ,T,low 


T^TT cA 2 + ^ x (V x cA 2 ) = l /AoKj J . (228) 

o atoms #p p j,T,low 

We see that the exciting terms for the first and second order correction to the 
vector potential are the zero-order atomic current density and the first-order 

•j 

correction to it respectively (except for the factor ). 

E o c 

We summarize the results of this section by writing down the equations of 
motion we have derived, only this time in the order in which they are applied. 

We can delete equations (213) and (217), since they determine only the second 
order corrections to the atomic structures, in which we are not interested. ‘ 

The center of our interest is the electromagnetic field, because of its technical 
importance. And here it is mainly the second order correction that we need to 
know. For it constitutes the response to the excitation provided by the zero- 
order approximation to the electromagnetic field. To begin with, we have the 
equations of motion for the zero-order approximation, namely 

“J'/vV 0 ' (229) 

» ( Mx 4). ^o) . ° , (230) 

3q 

g2 

‘ (Tetyg °Ao + v * (V X cAq) = 0 (231) 

Usually we do not have to solve these equations, because the solutions are 
given to us as the starting point of a technical problem. Then we determine 
the first order corrections from the following set of equations. 
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'"I" ,1 


*<? 


J ,T,low 


cJ 0 V • 


+ . -|(x2 ^///At 3f • ca ) , 

2 M 3<1 «3 p * 00 J ,T,low °P 


32 


+ V MV *cAi ) * —r Z /Aa |(xj^6 + X^o) 7 


(232) 

(233) 
(23U) 


/ 2 cA i T v w «*i/ - a ; uu p'aq 

t3ct ' o atoms Ipp 2 u J,T,low 

The pair of Schrbdinger equations (229), (230) and the pair (232), (233) have 

similar structures. The only difference is that inhomogeneous exciting terms 

have been added in the latter pair. Equation (23^) is not of great importance 

in our technical applications, because cA^ does not depend on the "incident” 

field cA Q , Finally we have the important equation for the second order 


correction CA 2 , namely 


(Vct") ' ^ ' q ^2 + ? * (* * ci 2) = 


~ Z /Act |(xo J 'i'l + Xi J if»o + Xo'h + xj'l'o*) ? 

o atoms # p p 


(235) 


J ,T,low 

In order to use it we have to first determine the solutions of (232), (233) 
and their complex conjugates. 

Sometimes it may be helpful to rewrite the equations (232) - (235) in such 
a way that only the symbol appears in them, instead of the symbol S, ^ 

The reason for this is that analytical expressions for are more easily 
formulated than those for J, _ . , In order to explain the modification that 

J , i , 10W 

we are going to maJce, we review the definition of T low- * We had (see 
equation (132)) 

J = Z Ji P + Z E . 

J “ J,a a J J,a a 


( 236 ) 
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By definition, Jj T 1qv Is the low-modes portion of the transverse part of J^, 


so that 


J j,T,lov ’ a j ow J j,a E a 


(237) 


By assumption, cAo is purely transverse and its high-modes portion vanishes. 
This statement implies that cAq may be expressed in the series 


cA 0 = I Q, 
a, low 


a,o a 


(238) 


According to the orthonormality theorems (98), (100 ) of Box #5» all those terms 
in the series (236) that are not contained in the series (237) contribute zero 
to the integrals in (232) and (233). Thus we may delete the subscripts T, low 
in these integrals. 

However, this mere deletion is not permissible in equations (23*0 , (235). 
Here we must proceed in a different way if we wish to rewrite these equations 
in terms of Jj instead of Jj ^ ^ ow « We start with equation (237). This time 
we write in the independent variables. Thus 


Jj T iow (qk » ?) = E J 1 a (qK)S a (?) * 

^ j,r,iow a, low ,a a 

We express the expansion coefficients in terms of integrals taken over the 
home-space . We obtain 

J JtE (q k ) = ///4T'(E a (?') • 

Here we have affixed primes to k* and r in order to avoid confusion in the 
equations that follow. For the reasons that were explained in the paragraph 
following equation (238), we may delete the subscripts T and low. Thus our 
last equation becomes 

k 


kv= 


(239) 


J j,a (q > = ///At *E a (?» ) • J^q ,r») 


(240) 
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On combining (237) and (24o), we obtain 

J 1T1 J*'). I E (r)///4T'i (?•) • 
j,T,lQV a ,low a a 

or 




jj(9 k »r') , 


^J.T.low^*^ s fff* r 't I ’ < 2ltl ) 

Q» y low 

This equation shows that 5^ ^ lov is obtained from Jj by means of a linear 
operation, called a projection (into the sub-space spanned by the low-modes 
patterns E & ). We denote this projection by the symbol Then the short- 

hand notation for (24l) will be 


J J,T,low “ P T,lOw J J 

Now we come to the integrals in (234) and (235). 
integral in (234) becomes 


(242) 

Here we use (242). The 


fa ftxo^o + xi*.V Tfl jr 4 

P 

But the projection P^ 1qw , being an integration over the home-space (according 
to (24l)) may be interchanged with the integration over the configuration space, 
so that the integral becomes 

P T,lcvfa * XO ♦Stfj • 

P 

We may even interchange the projection and the summation over the atoms # p. 
Similar statements apply to the integral in (235). Therefore, equations (232) - 
(235) may be replaced by 
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original 
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(if! .-(*„! ///ArJj . olo > p , (2U3) 

3q 

<* + 4k 7) - v *V>p * - 3<*i 7 IH^i • cS 0 >p * > < 2W( > 

3q 

oI * + ? * (? * e5l) * r? P T,low t l „ / 40 l<x; j *0 ♦ xW)j r ' (2>*5) 

o atoms /rp p 

3 2 

( T ct) 7 °* 2 + 7 * ( 7 * cA 2) * 

= 77 P T lQW l /Act |(xo J ti + xVn + X^l + X^o >Jj • ^ ) 

o * atoms jhp p 


The main part of our report is based on these equations. Of course, they are to 
be supplemented by the equations ( 229 ) - (231) for the zero-order approximations. 

It may be helpful to rewrite equations (243) - (246) in a more abbreviated 
form, in which they may appeal to the physical intuition. To this end, we intro- 


duce the abbreviation 


(U. ) = (i HJtxJ. • cA ) . 

j,o'p c n) J o'p 


(247) 


The quantity (U. ) is the vector potential component in the configuration 

J»o p 

space of atom #p. The index 0 indicates that (U. ) is based on the zero-order 

J »o 

approximation cA q to the electromagnetic field. Furthermore, we abbreviate the 

two integrals in (245) and (246) by > and <Ji> . These two quantities are 

o p ■ p 

the zero-order approximation and first-order correction to the quantum-mechani- 
cal expectation values (as indicated by the carets < >) of the atomic electric 
current density (in our home-space) of atom #p. With these abbreviations, equa- 
tions (243) - (246) assume the form 
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In fact, someone who is not too concerned about mathematical proofs would be 
inclined to write down these equations at the outset of an investigation. He 
could well say to himself: Things Just have to be that way, they cannot be 
otherwise. Perhaps he may have forgotten to include the projection operator 
P T . But, after some reflection, he would have spotted the need for it, 
otherwise cAj and CA2 would not turn out to be purely transverse and devoid of 
the bothersome contributions of the high modes. 



11, The Special Case of Free Electrons. 


jjJjj 

or poor quality 


The atoms (or molecules) that we have considered so feu: are, in general, 
composite entities consisting of several particles (electrons and nuclei). 

Now we are going to examine an especially simple example, in which each of these 
entities consists of a single free electron; no nuclei are present. The electron 
will be treated non-relativist'ically and quantum-mechanically. In this simple 
example, the configuration space of each electron may be considered as being 
identical with our three-dimensional home-space. We may therefore identify the 
three configurational coordinates q 1 , q 2 , q 3 of a given electron with the three 
cartesian coordinates x, y, z of its position. And we may combine q 1 , q 2 , q 3 

into the position vector q * q^ej + q 2 e 2 + q 3 e 3 , where the e. are the cartesian 

J 

unit vectors. The configuration of the electron is then characterized by its 
position vector q. We may therefore choose a more appropriate notation for the 
function p (q^,r) of Section l) and the functions Jj(q*,r) of equation (5). The 
new notation will be p(q,r) and Jj(q,r). Because the charge (-e) of an electron 
(e = l.,6 x lCf 1 ^ coul) is concentrated in a single point, the functions p(q,r) 

and J (q,r) have a particularly simple form, namely 

J 

p(q,r) = (-e)6(q-r) , (252) 

and 

Jj (q,r) = (-e)6(q-r) e^, J = 1, 2, 3 , (253) 

where 6(q-r) is the three-dimensional Dirac delta-function (really a distribution, 
whose meaning is given in terms of integrals). It has the following properties: 
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for any arbitrary smooth function f. 

We can now proceed to adapt equations (24$)-(246) to the case of free 

electrons. First of all, since each electron #p is free, the function V is 

zero. Of course, we neglect the direct interaction of the electrons with each 

other, just as we did for the atoms. This is permissible if the electron density 

is not too great. In the second place, the coefficients in (243) have rather 

simple values. Here we go back to the dicussion that followed equations (125, 

repeated) and (l4l, repeated),, which may be found after equation (l45)- There 

1 da^ dqk 

we see that the expression^ is equal to the conventional kinetic 


.J 


energy of an electron. And since we are using cartesian coordinates as the q , 

m for J * k, where m is the conven- 
tional electronic mass. Then with the aid of (253) and (254)* equation (243) 
becomes 


we conclude that * 0 for J ^ k, and 


(ih — 7 + m X|)_ - - (-e) “ (t_ e • cA ) . 

aq J 1 P c o j o p 

We may regard -^4- and x 1 ^ as the cartesian components jfj of the vectors and 
3< 1 

Xi > The last equation we wrote stands for three equations, one for each of the 
three values 1, 2, 3 of the index j. We can combine these three equations into 
one vectorial equation, namely 

(ttivtyl + rnjci ) = -(-e)~(t|> cA ) ■ (255) 

p C O O. p ’ 

So far the independent variables on which the functions in (255) depend are 
<5> p ,t. (The ? that occurred In ( 253 ) disappeared hy virtue of the integration 
over the home-space, ) But since the configuration space of each electron #p is 
identical with our home-space, we may Just as well use r,t as the independent 


variables • 
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We use the delta-function expression (253) for the in the remaining 

equations (2Uh) - < 2U6 ) . And we also use vectorial notation. For convenience 

we add the equations for the zero-order approximations (ij» ) , (v ) , and cA . 

o p o p o 

We then arrive at the following system of equations, where, as before, the 
subscript p Indicates that the quantum-mechanical functions refer to electron #p. 


+ mx o ) » 0 , 

(256) 

. 3<j» .y 

•’‘o’p* 0 • 

(257) 

(Mi?*! + mx!) p > - (* Q ) p i(-e)oA o , 

(258) 

♦ 4 ' • * - R>p • ^- e)c5 o • 

(259) 

fact) 2 =*„ + ’“ (7 x eS o ) - 0 , 

(260) 

TSny =A, + 5 » (? » c Al ) - £ 1 c * T lOT I (• 

0 * electrons 

- e >K*o + Vo>p ’ (z6l) 



W CA 2 + ? X (7 K CA 2 ) - ^ P Tilw ei J tron (• 

-e)|(x**l + ♦ x/i + 

fo 

(262) 


Note that the projection operator p^ which selects the low-modes part of the 

transverse portion* appears in (261 ) and (262). Ag^in we observe that the excita- 
tion terms for the vector potentials cA\ and cA 2 are (apwt from the operation 
P _ . and the factor -^— ) the quantum-mechanical expectation values of the 

1 • XOW 6 c 

o 

current densities. Equation (256) - (262) are the basis cf the main part of 
this report. 

It may be helpful to adapt also the expression (203) fo the action func- 
tional to the case of free electrons, so that equations (256) - (262) may be 

* 

As we did it before, we have omitted the subscript "low” of the symbols 
cAo, cli, cl2. 
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derived also directly without the detour through the problem of composite 
atoms. In this derivation one uses the same rules of the calculus of variations 
and of the calculus IM’ perturbations that we have employed before. For the sake 
of convenience, we replace the triple integrals /// by a single integral /. On 
the other hand, for the sake of clearness, we will have to indicate the indepen- 
dent variables. In (203) ther'e occur the quantities U.(q,t), which are given by 

J 

(204). We repeat the latter equation, but omit the index p. 

Uj(q*k) * ^ /At cA low (r,t) • Jj ,t,1ow^*^ * (20U)» repeated 

Then, with (2l*l), we obtain 

Uj(q,t) « ~ //AxAT'cA low (r,t) • l E a (r)E a (r’ ) • Jj(q,r») . 

& y lOW 

Now we use the delta- function expression (253) for J.(q,r f ). It permits us to 

perform the integration over r’ in a simple way. The result is 

U.(q,t) = ^ (-e)/Ar cA. (r,t) • £ E (?)E (q) • 5 . 

a, low d 

Now, Uj (q,t ) is the cartesian component of the vector U(q,t). This vector 
is therefore given by 

U(q,t) = | (-e)/ArcA low (r,t) • J S (r)E (q) , 

a, low 

or, with a change of notation for the independent variables, 

U(r,t) * “ (-e)/At , cA 1 (?' ,t) • J E (r ? )E (?) . 

a, low 

We write this equation more succintly as 


0 = - (-e)/AT f cA* 


low 


Z 

a, low 


E'E 
a a 


(263) 


We use this expression in (203). At the same time we make the changes of 


notation that we made before. For instance, we use the vector % instead of its 
cartesian components > and we introduce the electronic mass m. Of course, we 


set V equal to zero, since we are dealing with free electrons. The final result is 




a » b t*a 
appr, quant 


/ b dt{ ° / AT fgct c5 low^ ’ C 3ct c5 low ] 
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Po °L>e 




--|/AX[VX cA low ] • [»* 0A lov ] . ♦ 

n 
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#P 


+ f X* 
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XL ♦ ^(-e)//AT»ATcA‘ ow • I fi;* a • «*♦ + *♦) JJ} 

^ a, low 


(26)4) 


When we apply the action principle to this action functional, we obtain the 
following set of equations in the usual manner: 

+ ”X ♦ * p I * a ®; • - o , 

y * * a, low 

3i|/ . 

& ll + — 5 • x„ + o7<- e >x„ '/ lT l KK • oS iov " 0 • 


'p 2c 


a, low 


a a 


a 2 * 

(act) 2 cA iow 


V X (7 X cA loV ) = ^ (-e) I /Ax' J E ft E 


e c 
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electrons 


a, low 


a a 


(?V +?<!'**) 

p P P P 


- ~ (-e) P m 1 _. I &}LK + XJ& J • 


& c 
o 


T > low electrons 2 ^ p p P 
#P 

To obtain the last equation we had interchanged the primed and unprimed 
coordinates in the double integral of (2 64). Since ci^ is a superposition of 
terms comprising only the low transverse modes, we have 


/At * T EE** cA* . = cA 

* u on 

a, low 


a a low '"'low * 


Thus, the set of three equations simplifies to 

= - * p | (-o)cA lov , 

ifi P + - ■ ■ V • y s — “ ( —g ) y_ • 
111 3t 2 X P 2c v e ' x p 


cA, 


low 


(265) 

(2 66) 

(267) 


7i5ET* c5 iow * ? * «*!*,> ■ To <- e > P T ,iow j„„ + V V (268) 


_ # 


electrons 


ORIGINAL PAGE IS 
OF POOR QUALITY 


A99 


These three equations are the exact equations of motion for the action functional 
(264), whereas the set (256) - (262) results from the perturbation calculus. 
Equation (268) shows that the excitation term for the vector potential is the 
quantum-mechanical expectation value of the current density, except for the 
factor . Of course, we must not forget to use the projection operator 

EC 

O 

P_ . Its presence ensures that cA, continues to be a linear combination of 
T,low low 

low transverse modes only, if it started out that way. 

If we combine equations (266) and (267) we arrive at the usual form of 
the Schrbdinger equation for electron #p, namely 

si *p = * ( -*^ - i ( - e)c! W*p • (269) 

In the entire development of this section we have treated each electron, as 
an independent entity, just as we did it for the atoms that were described by (203) 
In mathematical terms, this means that we regarded the quantum-mechanical func- 
tions of the many-electron problem as products of single-electron functions in 
the manner of equations (l?4a) - (194c). This is not quite correct because our 
procedure did not take account of the Pauli exclusion principle, which requires 
that, instead of products, we should have used determinants. But the pursuit of 
this matter is probably not worth the effort. 
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Note #1. Derivation of the Relation r J 4 - 1 = $ x (7 j. x j.) . 

3d 1 J 3qJ 1 P i J 

We start with the requirement that a tagged element of charge which happens to 
he at position r for the atomic configuration {q*} ought to with a velocity w 

given by 

pv = J , (Nl,l) 

where, of course, p and J are evaluated at (q d ,r). 

Next we consider a particular element which happens to he at the position 

a(q ) in our home-space when the atomic configuration j.s {q }. We chose the 

symbol a, not r, because r stood for the generic position, which ranges over the 

whole room that contains the atom and the radiation field, whereas a refers to 

Just one discrete and tagged point. Now we evaluate the velocity ^ a of this 

tagged point. Here we have to be careful with the notation for partial deriva- 

3 i - 

tives. We have used the symbol — rr when the q and r were the independent vari- 

i - i 

ables. But a depends only on the q , the r does not enter. So when only the q 

(not the q* and r) are the independent variables, we use the symbol — ("# slash") 

■ 

to denote partial derivatives. Then we have 


_d - ga dai ga i 

dt “ " - i dt " *_i v 

fq 




( N1 , 2 ) 


But according to (Nl,l) we should have 

4a = ■ — -■ - — s — J(q j ,a(q^)) * ; — J ,(q J ,a(q J ) )v i , (Nl,3) 

dt P(q j ,£(q j )) p(q J ,S(q J )) * 

where we indicated, how the p and depend on the n configurational coordinates q*^ . 

4 4 _ 4 

This dependence is not only explicit, as indicated by the first q J in (q ,a(q J )), 
but also implicit by way of the position r = a(q^), as indicated by the second q^ 
in (q^,a(q^)). On comparing (Nl,l) with (Nl,3) and remembering that the v d might 
be any arbitrary set of configurational velocity components we obtain 
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J& = 


J (q. J t l(q J )) 

Iq 1 p(q J ,S(q J )) 1 




J 


Now we take the partial derivative of this equation with respect to q , i.e. we 
take *77 , (not — -r). Here we must remember that the dependence on the q' J is not 




aq 1 


only explicit, but also implicit. In order to make the ensuing equation less 
cumbersome to write , we leave out the symbol ( q^ ,a( q~ ) ) for the independent 
variables. We obtain 

J_ Ji = -L- (I j ) + Ji. v (ij ) . 

2 q J Jq 1 Sq J P 1 2 q J P 1 

We express the - £ r > an the right-hand side with the aid of (N1,U) written for 

W 

the index j instead of i. We then get 

% tfk _ ■ 9 /I •s v l - 
p 

Similarly (interchange i and j). 




J- Ji= 9 . W Ij) . 

Jfq 1 0 q J aq 1 P 5 P 1 P J 


(Nl,5a) 


(Nl,5b) 


According to the principle of gene-identity the position a of the tagged 
particle is a function of the atomic configuration {q*} alone. A necessary and 
sufficient condition for this is that the two mixed partial derivatives on the 
left-hand sides of (Nl, 5a of b) have the same value. Equating the two right- 
hand sides yields 

-£) ( f V + f • ?( 7 V ■ i? ( iv + s- 5 i • ?( f V • (N1 ’ 6) 

This equation may be brought into a form that will be more useful to us 
later on, namely into the form shown in the title of this note. We multiply 
(Nl, 6 ) by p and differentiate out the products. At the same time we replace 

a- 

- — 7 and 7 by V • J. and 7 • J. respectively, according to equation (T) 

3q u 8q J 
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of the appendix. We then obtain 


7 J.V • 5, + T + I j • * 3 ; Cv ~)j, => 

Pi JgqJipJ i J .pi 


,ij.v . J 1 + -^J 4 •* jj *J t . . 


p J 

or after rearranging, 


-i- j - -2- j = (V i) x (J x J ) + i V x (J. x J ) . 

3q J 1 p i J P i J 

The two terms on the right may be combined. Thus finally 

j - -L. j = 7 x 4 j. x J.) . (Nl,7) 

3q 3q J P 1 J 

And this is the equation shown in the title. 

Equation (Nl,7) must hold wherever there is an element of charge that may 
be tagged, i.e. at all points r of the home-space where the charge density p does 
not vanish. But wherever the charge density p (and then also and J^) is zero, 
the equation is fulfilled by default. It simply states ”0 = 0". Thus (Nl,7) is 
fulfilled everywhere. 

The reader who is satisfied with this proof may stop right here. Others, 
however, might feel more comfortable if they saw a second and different proof. 

In this second proof, we use the same set of independent variables, namely the 
q 1 and r, that was used in the main text of the appendix. So there is no need 


for employing the slash-derivatives 


JL 


i ' 


Instead of tagging .fust one element 


of charge, as we did before, we tag all elements at once. Since the charge 
elements form a continuous distribution, we have to employ a triplet (C 1 ,? 2 ,? 3 ) 
of quantities to tag each particular element. For instance, for the £ a (Greek 


indices range from 1 to 3, whereas latin indices range from 1 to n.) we might 


use the cartesian coordinates of the element when it was placed according to 

some arbitrarily chosen fiducial configuration of the atom. We then have to 

deal with the three functions K (q »r). They tell us which tagged element is 

at some chosen location r for the atomic configuration {q 1- }. 

The principle of gene-identity allows us to evaluate the velocity v(q ,r) 

of an element of charge. We know that the 5 a -tags of an element never change 

dE° 

when the element is in motion. Thus the material time derivatives -4r- (taken 

dt 


as we move with the element) must vanish. But 

's,a 


i ; a = i r a + v ■ (VE U ) 
dt * 3t ■ W / , 


(Nl.ff) 


3 „a 


where tt* 5 is the local time derivative (taken at fixed position r), given by 


3t 

a o. _ aq 1 as 01 1 ae a 
* 4t a , 1 aq 1 ‘ 

d ot 

The requirement that ^ £ be zero yields 


(Nl,9) 


w • (VC°) = -v^-^t* . (Nl,10) 

3q 

Equations (Nl,10) are a system of 3 scalar equations (one for each a = 1, 2, 3) 
for the velocity vector w. Hence w is determined. If we were to solve these 
equations (but we are not going to solve them) we would find that w is a linear 


function of the n velocity components v . Thus we can write 


i- 

W = V w. 


»r) , 


(Nl,ll) 


where the n functions w^(q^,r) depend only on the atomic configuration {q^Jand 

the position r. On combining (Nl,10) and (Nl,ll) we obtain 

i- = c a _ i 3C a 
v w. • = - v — 1 —* . 

Sq 1 

But this equation must hold for any choice of the configurational velocity 
components v*. Thus we obtain 
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4--v t • H a . 

aq 1 1 


(HI, 12) 


We differentiate this equation with respect to q d and obtain 


3q J Sq 1 


dwi „ aP o 

— • Vt “ - w • V . 

1 -• W. V . , 

3q . 1 3q J 


— 3 3f u 

where we have interchanged V and — -r . We replace the — *r on the right-hand 

- - a 3q J 3q J 

side by - w^ • V£ , according to equation (Nl,12), written for the index J 
instead of i. Thus 

a a 3w. 

3 1L. L , VC® + w, • V|w • V? ) , 

3q d 3q 3q d 3 

or, on differentiating out the product in the last term, 

— ~ • ^5° + w. • (Vw.) • ^K a + wi • (775 a ) • w. . (Nl,13a] 

3q J Sq 1 3q d 1 d 3 

Similarly (interchange i and j), 

. 3w 

-rr-^r f • 1Z a + w • (7w.) • $5 a + w • (vv5 a ) • w . (HI, 13b! 

3q X 3q d Sq 1 3 3 1 

Now we subtract (Nl, 13b) from (Nl, 13a). The terms on the left-hand side will 

t 

cancel, because mixed partial derivatives do not depend on the order. Also the 
last terms on the right-hand side will cancel, because VV? a is a symmetric 
tensor. Thus, after rearranging and factoring out the common term • V5 a , 

{~r w. + w. • Vw. - — w. - w. • $w. } • 75° = 0 . 

3(l l J i J 3q j i J i 

This single printed equation really stands for three equations * one for each 
a = 1, 2, 3. But the three vectors 75° are linearly independent. Thus the 
expression in the parenthesis must vanish, and we obtain 


-~r v. + w, • Vw. = — w, + w. • Vw. . 

Bq J i J i 3q i J i J 


(Hl.lfc) 
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From the equations (Nl,l), (Nl,ll), and ( 5 ) of the appendix, repeated here 
for convenience, 


pw * J , 

(Nl.l) 

i- 

w* vw 1 , 

(»1 ,11) 

T ir 

j 1 v Jj , 

(5) 

pw^v 3- = JjV 1 . 



we conclude that 


But this equation must hold for any choice of the configurational velocity 


components v*, so that 


1 T 

w, = — J . . 
i P i 


(Nl,15) 


Inserting this into (Nl,l4) yields 


-J- (i ty ♦ i jj • 5 4 » - $ ty ♦> J,) <“.*> 

But this equation is the same as (Nl, 6), from which the desired relation (Nl,7) 


follows by means of purely algebraic manipulations. 
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Note #2. Vectorial Derivation of the Equations 

— A mm mm ~ mm mm mm 

6J * H" ^ plSr ) + v * ($? * J) and 6p * - V • (p6r) . 

One may always express the charge density p and the current density J in 
terms of three scalar functions a(r,t), 0(r,t), y(r,t) of the position r and 
the time t. These expressions are 

p » • 7M 7y , (N2,l) 

and 

(N2,2) 


j = . |2. x |i 7 Y x Va - Va x 70 . 


3t 


3t 


One can readily verify that the continuity equation 

f“ + V • J = 0 . (N2,3) 

is automatically satisfied. The triple {a,0,y} may be regarded as a tag for an 
individual element of charge. 

Sometimes one may wish to express p and J in terras of several triples 
{ V S n’ Y n } > as ln 


- I 


n 


70 x- Vy , 
n .'n * 


(N2,4) 


and 

3a 00 3y 

J ■ I { - K * Ft K * % - 17 K * 9 <y • (N2 -5> 

However, one such triple suffices. At any rate, our derivations will employ 
only linear operations. Thus what is true for one triple will also hold for 
the more general expressions (N2,U) and (N2,5). Therefore, we stay with the 
simpler forms (N2,l) and (N2,2). 

Now let us subject the three functions a, 0, y to the variations 6a, 60, 6y. 
The resulting variations 6p and 6J of p and J are then given by 


A107 
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6p ■ ($6a) • x + (V«0) • Vy x Va + (76y) • tfa * 70 
and 

« - (“ da)* V6 x 7y - || (fog) x 7y - |2- $0 x (?6y) - 

- (~ 60) * * Ya - || (^6y) * - f£ * (76a) - 

- (gf 6y) • Va x $0 - |£ (76a) x ?0 - |f 7a x (760) . 

By means of a straightforward, though tedious, calculation one shows that the 

preceding two equations may be written in the following form. 

« 

6p * 7 • [6a 70 x 7y + 60 7y x 7a + 6y tfa x 70] , (N2,6) 

and 


6? « - ~ [6a 70 x Vy + 60 ^y x 7a + 6y 7a x 70] + 

+ ? x [6a (|| 7y - |* ?0) + 60 (|j 7a-- |2. 7y) + 6y (|^ 70 - ||- 7a)] . (N2,7) 

So far we left the variations 6a, 60, 6y unspecified. But now we shall 
introduce specific variations related to shifts 6?, 6t in the event spac^, 
i.e. the four dimensional space of position r and time t. We shift the three 
function patterns a, 0, y by 6r and 6t. The variations 6r and 6t may still 
depend on the event r,t. Then, for each fixed (r,t), the variations 6a, 60, 

6y are given by 

6a = - • 6t - 5r • 7a , 

60 = - 6t |~ - 5r • 70 , (N2,8) 

6y = - 6t ^ - 5r • 7y . 

Note the minus signs, which come about in the following way. Since we shifted 
the functional pattern a, we know that the new a at the shifted event r + 6r, 
t + 6t is equal to the old a at the original event. Thus we have 
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(a + fia)(r + fir, t ♦ fit) * a(r,t) 

On expanding the left-hand side "by Taylor's theorem and retaining only the 
zero-order and first order terms, we obtain 

o(r,t) + fit + fir • + 5o * a(r,t) 

On cancelling a(i?,t), we obtain the first equation of (N2,8). The other two 
equations result in a similar way. 

We insert the expressions (N2,8) into (N2,6) and (N2,7) and use equations 
(N2,l) and (N2:,2). With the aid of the relation 

(Vo)(Ve X Yy) + (v$)(7y X 7a) + (Vy)(Va x 70) * (Va • $B x 7y) 1 , (N2,9) 

where I is the unit tensor (or idemfactor), we have 

6a VB x Vy + fiB 7y x Va + fiy Va x 7S “ - pfir + Jfit (N2,10) 

In the expression behind the curl sign ?xof (N2,7), the terms in fit will cancel, 
while the terms in fir combine to fir x 7. Altogether then, we obtain 

fip * - V • (pfir - Jfit) (N2,ll) 

fi J » — (pfir - Jfit) H x ({? x j) (N2,12) 

In the special*case fit » 0 (no shift in the time direction), we obtain the 
equations in the title of this note. 


